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FOREWORD 


This preliminary report on GEMINI Ascent Guidance Simulation is 
submitted to the McDonnell Aircraft Corp. (MAC), St. Louis, Missouri, in 
partial fulfillment of requirements under Task 2, Booster Guidance Studies, 
Statement of Work, Purchase Order Y20163. 



ABSTRACT 


This report presents the results of preliminary studies on the simu¬ 
lation of GEMINI Ascent Guidance. These studies have been made to simu¬ 
late the Stage 1 (open-loop) and Stage 2 (closed-loop) guidance equations to 
verify the compatibility of the computer equations with the vehicle and con¬ 
trol system and, to some degree, to determine the accuracy of the equations. 
The simulation is described in detail along with derivations of equations used 
in several areas of the simulation and operational math flow. 

The major results of the simulation studies discussed in Section III 

follow: 


• The ability of the guidance equations to control the ve¬ 
hicle has been demonstrated. 

• The navigation errors in the guidance equations are 
small. 

• The equations implemented to steer the vehicle during 
Stage 2 will bring the vehicle to nominal insertion con¬ 
ditions. 

In several areas, the simulation will deviate from the most recent 
math flow release. These deviations are indicated in the appendices. Plans 
for further study and refinement in the simulation are presented in Section IV. 
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Figure 1-1. GEMINI Ascent 




Section I 


INTRODUCTION 


This report presents the results of preliminary studies on the simula¬ 
tion of GEMINI Ascent Guidance at the IBM Space Guidance Center. These 
studies have been made to simulate the Stage 1 (open-loop) and Stage 2 (closed- 
loop) guidance equations to verify the compatibility of the computer equations 
with the vehicle and control system and, to some degree, to determine the ac¬ 
curacy of the equations. Figure 1-1 presents a typical ascent profile of a 
GEMINI mission. A detailed description of this profile is attached as Appen¬ 
dix E. 


The simulation is described in detail along with derivations of equations 
used in several areas of the simulation and operational math flow. 

In summary, the major results of the simulation studies discussed in 
Section III are as follows: 

• The ability of the guidance equations to control the vehi¬ 
cle has been demonstrated. 

• The navigation errors in the guidance equations are small. 

(See Figures III-14, -15, -17, and -18, and Table III-I in 
Section III.) The platform Z-axis velocity and position 
components differ from those in the environment; however, 
the addition of the deflection of the gravity vector will cor¬ 
rect this situation (Appendix A). 

• The equation implemented to steer the vehicle during 
Stage 2 will bring the vehicle to nominal insertion condi¬ 
tions. 

In several areas, the simulation will deviate from the most recent math 
flow release. These deviations are indicated in the Appendices. Plans for 
further study and refinement in the simulation are presented in Section IV. 
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Section II 


DESCRIPTION 


A. INTRODUCTION 


Initially, the GEMINI Ascent Simulation was to include the use of IBM's 
combined analog-digital simulation facilities with a Pace analog computer for 
instrumentation of autopilot and rigid body equations. Subsequently, however, 
IBM determined that an all-digital simulation was feasible, and study was 
defined and initiated accordingly. As a result, Stage 1 simulation was con¬ 
ducted on both the combined and all-digital programs; Stage 2 simulation was 
conducted only on an all-digital program. Due to the organization of the pro¬ 
gram, the only redundancy in the Stage 1 program occurs in the autopilot and 
rigid body dynamic equations. The remainder of the programs are identical. 

The entire simulation is based on the use of a main program, which de¬ 
fines the external world or environment (including time) into which a vehicle 
model is inserted along with a self-contained guidance computer. Each of 
these three subjects — main program,, vehicle model, and guidance computer — 
are discussed in this section. Figure II-1 is a functional drawing of the pro¬ 
gram. Dotted lines on the drawing indicate the relationship between Figure 
II-1 and these three subject areas. 

B. MAIN ENVIRONMENTAL PROGRAM 


Data constants are read into the computer to initiate the program on the 
ascent problem. The main program computes the total vehicle acceleration — 
gravitational and perturbative. The perturbative accelerations are provided as 
accelerometer inputs to the guidance equations which, in turn, provide inputs 
to the autopilot. The vehicle angular rates, as outputs of the autopilot, are 
then transmitted to the main program. These rates are used to update vehicle 
attitude, including the platform gimbal angles, which are also provided as an 
input to the guidance equations. 

The main program integrates total acceleration in the platform frame 
and updates position and velocity values. Aerodynamic forces are also com¬ 
puted in the main environment and are supplied to the vehicle model equations 
for inclusion in the solution of the rigid body equations. 
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Figure II-1. GEMINI Ascent Guidance Simulation 
































In order to provide the ability for changing parts of the program with¬ 
out affecting the rest of the program materially, as many functions as possi¬ 
ble were incorporated as subroutines to the primary program. Also, the 
program for the guidance equations and vehicle model were written as sub¬ 
routines to the main program. In Figure II- 1 , note that functions such as 
aerodynamic data, gravitational acceleration, atmospheric coefficients and 
vehicle attitude determination are logical choices for inclusion as subroutines. 

The coordinate systems employed in this simulation for presentation of 
results concur with GEMINI coordinate systems defined in the compatibility 
data sheets and various GEMINI Specification Control Drawings. The three 
basic coordinate systems include the platform, spacecraft and launch vehicle 
references. These systems as well as the various intermediate and compu¬ 
tational frames are described in Appendices A and D. 

Since computer subroutines previously written were available early in 
the program for the mathematical models which were used for atmospheric 
and gravitational terms, the subroutines were used without modification. 

Some variation exists between the models used and those intended for GEMINI 
simulation. However, specified models can be incorporated into the program 
by inserting the new subroutines, when available, without further modification 
to the program. 

The basic timing of the program is such that one computer cycle is the 
equivalent of 0.05 sec. This value was chosen on the basis of limiting round¬ 
off error while still providing updated platform gimbal angles to the fast-loop 
guidance equations at the required rate (20 times/sec). Note that the slow 
loop guidance equations are entered once every ten computer cycles (2 times/ 
sec). This closely estimates the real time computation cycle that will be ex¬ 
perienced in the GEMINI computer. 

Certain effects, assumed not to affect the simulation, were omitted: 
Thrust variation with change in altitude, body bending modes, fuel slosh, 
aerodynamic damping terms, gravitational oblateness beyond second harmonic, 
platform errors, misalignment errors, and errors in rate and angle terms 
resulting from displacement of instruments from the vehicle's center of gravity. 

A complete description of the simulation program can be found in Ap¬ 
pendix A, and the math flow and program listings are in Appendix D. 
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C. VEHICLE MODEL 


The vehicle model is conveniently divided into two parts. The first con¬ 
sists of the autopilot and engine servos; the second deals with the solution of 
the rigid body dynamics. In defining these parts, IBM had to rely on available 
data which, unfortunately, has been revised subsequent to writing the program. 
The data affected by these revisions is associated primarily with the aero¬ 
dynamic characteristics. Complete data sheets are available in Appendix G; 
however, the sources of this data are obscure at present. Although the data 
is not consistent with most recent compatibility data, the variation is small 
and represents a physical model differing only slightly from the currently de¬ 
fined GEMINI Launch Vehicle. In the case of the autopilot and engine servos, 
the model used agrees with the latest available compatibility data except that 
it has been simplified to a fourth order polynomial. 

The Stage 1 vehicle model was programmed in two mediums - analog 
and digital. These analog and digital programs are described in Appendices 
B and C, respectively. In each program, the model is reached by the main 
computer through a specific subroutine. The subroutine varies to provide 
the proper interface and real time control for the analog or digital model. 

While the rigid body dynamics is an integral part of the vehicle model, 
aerodynamic forces and location of the vehicle center-of-pressure are com¬ 
puted in the main environment and supplied as inputs to the vehicle model. 

The various vehicle constants (center-of-pressure location, aerodynamic 
coefficients) are stored as tabular data in the computer, and linear interpola¬ 
tion is used to evaluate specific values. 

D. GUIDANCE COMPUTER 

The model used for the guidance computer was formed by modifying the 
GEMINI Computer Ascent Guidance Math Flow as described in Appendix E. 
Additional logic was provided to permit entering specific parts of the program 
as a function of time. For example, in the case of fast-loop versus slow-loop 
computations, a simple counter is used to gate program flow around the slow 
loop nine out of ten times. An additional subroutine provides the discrete in¬ 
puts required by the guidance computer. Due to the sequential nature of 
writing and checking the simulation program, Stage 2 guidance equations are 
written in the form of a subroutine to the Stage 1 equations. 

In writing the programs for the guidance equations, some variation is 
expected in the actual performance of the guidance computer and the program. 
Some of these variations are as follows: 

1. The simulation program is written in floating-point arith¬ 
metic with eight significant places. 
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2. Time increment for the slow loop is an exact multiple 
of the fast loop and is constant. 

3. Time increment for both Stage 1 and Stage 2 slow-loop 
computations was assumed to be 0. 5 sec. 

4. Update logic was not exercised. 

5. Orbit Velocity Adjust Equations were not simulated. 

6. Quantization of the attitude error outputs by the D/A 
converters was not accounted for. 
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Section III 


RESULTS 


A. INTRODUCTION 

The primary purposes of the simulation are to prove (1) the adequacy 
of the GEMINI IGS guidance equations and (2) compatibility between the com¬ 
puter, launch vehicle autopilot, and control system. The work performed to 
date has been devoted primarily to obtaining a functioning model to satisfy 
these goals. Certain refinements to the current model are required (see 
Section IV); however, a working model has been defined and has produced pre¬ 
liminary data. These data are presented in this section. The specific data 
runs to be considered are defined as follows: 

• Run A — Combined analog-digital simulation, 0 to 150 sec., 

Stage 1 (Guid 1) equations. No aerodynamic forces con¬ 
sidered. See Figures Til-1 through III-4. 

• Run B — Digital simulation, 0 to 150 sec , Stage 1 equa¬ 
tions, aerodynamic forces included. See Figures III-5 
through HI-19. 

• Run C - Same as Run A except aerodynamic forces in¬ 
cluded. See Figures III-5 through III-19. 

• Run D — Digital simulation, 0 to 333 sec., Stage 1 and 2 
equations used. Aerodynamic forces included for Stage 
1 only. See Figures III-20 through UI-28. 

B. STAGE 1 

1. GUID 1 NAVIGATION 

A simplification was incorporated into the program which affected the 
calculation of the Z components of gravity (see Appendix A). In essence, 
this simplification neglects the fact that the actual platform is aligned to the 
local gravity vector rather than to the geocentric radius vector. Consequently, 
considerable error is experienced when computing the Z component of gravity 
in the platform frame (Figures III-16 and III-19). 
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Note that in Figures III-14, -15, -17, and -18, the errors experienced 
in velocity and position computations are within acceptable bounds. Specifically, 
the errors are as follows: 

1. The maximum error in V x occurs near 90 sec and attains 
a value of 0.07 ft/sec. This error then recedes to 0.025 
ft/sec at the termination of Stage 1 Guidance. 

2. The maximum error in Rx occurs near 120 sec, reaching 
a value of 2. 7 ft (average V x error of 0.02 ft/sec). This 
error then recedes to 1. 6 ft at the termination of Stage 1 
guidance. 

3. The error in V Y increases as a function of time, attaining 
a value of 0.14 ft/sec at the end of Stage 1 guidance. 

4. The error in Ry also grows as a function of time, reach- 
in a value of 326 ft at the end of Stage 1 guidance. This 
discrepancy is unrealistic since the error in velocity, 

Vy, is 0.07 ft/sec. or less. This apparent conflict can 
be attributed to computation round-off error in the en¬ 
vironment equations, which are being solved at 10 times 
the rate of the IGS equations. Solutions to this problem 
have been established and will be included in future runs. 


While the guidance equations do not result in appreciable navigation 
errors, remember, that the computations were programmed in 7090 word 
size with floating point arithmetic. The shorter word and fixed point opera¬ 
tion in the actual GEMINI computer may result in more error. This is dis¬ 
cussed in Section IV, "Plans for Continued Simulation Study. " 

2. SIMULATION NEGLECTING AERODYNAMIC FORCES 

The primary purpose of effecting simulation Run A was to provide vari¬ 
ous checks on the program. In addition to simplifying program debugging, 
the various autopilot gains were verified. Figure III-1 illustrates the rela¬ 
tionships between attitude errors, vehicle rates and engine deflection angles. 
Figures III-2, -3, and -4 are records of the angles of attack, velocity and 
position during ascent. Note the yaw angle of attack (ay) in Figure III-2. 
Ideally, this angle should start at zero, but due to a singularity in the com¬ 
putation equations at t = 0 (ay = tan-1 Vyb/V x b, V x b = 0 at t = 0) any spuri¬ 
ous velocity computed for Vyb produces a finite angle of attack. This condi¬ 
tion quickly corrects itself as V x b obtains a finite value. The environment 
body -platform matrix is not updated during holddown, and can result in 
an error approaching 0. 3 ft/sec. in Vyb at lift-off. At t = 2 sec. , this would 
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result in an indicated angle of attack of approximately 0. 3 deg. (V x b«Vy « 50 
ft/sec.). Further investigation and corrective action on this error will be 
performed even though its effect on the simulation is small. 

3. REPEATIBILITY 

The data runs presented in this report were compared to similar runs to 
demonstrate the repeatibility of results. The difference in data obtained from 
these comparisons was always two orders of magnitude less than the value of 
the function,or, in the case of zero magnitude, the difference approached the 
machine noise level (quantization). Since the largest variations were experi¬ 
enced in angles of attack, the actual comparison of this parameter is interest¬ 
ing to note. The maximum error occurs at 108 sec. (max a n) and has a 
value of 0. 08 deg. when the angle of attack a is 6. 3 deg. 

4. AUTOPILOT DELAY 

Provisions in the digital autopilot permit the effective time of data trans¬ 
fer to vary from the guidance equations to the autopilot model. A computer 
run was made with zero time delay and compared with Run B, which includes 
a 20-msec delay. This was done in an effort to identify the differences be¬ 
tween the combined and pure digital simulations. The results of this run have 
not been plotted, but the behavior of the difference terms between Run B and 
this run are similar to those between Run C and Run B. 

Due to the mechanization of the guidance equations, there is a 50 msec 
(one computation cycle) delay in sending the attitude error signals to the analog 
autopilot. Thus, the combined simulation results lag the digital results. This 
is shown in Figures m-8, -9, and -10. The maximum variation in pitch angle 
of attack, which occurs immediately after autopilot gain change, is approxi¬ 
mately 0.3 deg. 

5. GUID 1 CONTROL 

One of the primary objectives of the simulation program is to prove that 
the guidance equations, as programmed, are able to control the attitude of the 
launch vehicle. The various functions plotted in Figures III-7, -8, -9, and 
-11 demonstrate this controllability. Areas in the referenced Figures re¬ 
quiring additional comment are discussed below. 


a. Roll Program Transient 


The step seen in the roll attitude error signal (Figures III-1 and -11) at 
the conclusion of the roll program is due to the fact that the duration of the 
roll program is fundamentally controlled by the 0. 5-sec. slow-loop computation 
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cycle. Essentially, the start and stop of the desired roll program have been 
computed correctly; however, the equations are such that the duration of roll 
maneuver is controlled to some multiple of the slow-loop computation cycle 
(0.5 sec.). The step may be undesirable, and will be corrected by setting 
4>FL = <f >f in Block A162 of Figure E-1A in Appendix E. 

b. Pitch Angle of Attack Perturbation at Gain Change 
(a^ Excursion at 108 sec.) 


In Figure III-8, note the magnitude of angle of attack, a n, at 105 sec 
(just prior to gain change). Note also the perturbation on angle of attack fol¬ 
lowing gain change. A detailed check on these results as well as discussions 
with other Contractors doing similar work indicate that the perturbation is 
reasonable. Presumably, the three-step pitch profile, which is expected to 
be available for the launch vehicle in February 1963, will minimize some of 
these effects. The pitch maneuver used in the present simulation is a two- 
step profile which approximates a gravity turn. 

c. Pitch Rate 


In Figure III-11, the pitch rate is slightly erratic in the time period 
40-70 sec. This departure from the desired constant pitch rate occurs at a 
time when the angle of attack has an appreciable value and aerodynamic pres¬ 
sure is approaching a maximum. A detailed analysis of this variation will be 
performed in an attempt to explain this behavior. 


d. Yaw Attitude 


As mentioned above on the results of simulation which neglect aero¬ 
dynamic forces, a change will be made to the simulation to remove the ap¬ 
parent errors in yaw angle of attack. Further discussion on this subject is 
not warranted here except to repeat that this error is attributed to the simu¬ 
lation rather than guidance equations, and the results of the error are small. 

C. STAGE 2 


Vehicle staging is simulated by assuming that Stage 1 shutdown, Stage 
1 and Stage 2 separation, and Stage 2 ignition occur simultaneously after 148 
sec. following lift-off. A more accurate simulation of this operation is being 
considered; however, it is not important in terms of the over-all IBM ascent 
simulation objectives. 


A comparison of the environment and IGS navigation data at the time of 
Stage 2 ignition as well as at other discrete times during flight is contained 
in Table III-1. The data compares favorably. The major portion of the po¬ 
sition difference on the Y-axis is believed due to round-off errors in the en¬ 
vironment navigation equation and not in the IGS equations. The difference 
along the Z-axis is due to the difference in velocity obtained by the environ¬ 
ment and the IGS equation. Both of these subjects are discussed more com¬ 
pletely in the Stage 1 portion of this section and in Appendix A to this report. 

The various programmed constants used in the Stage 2 guidance equa¬ 
tions are contained in Appendices E and F. 

Some of the initial conditions on the vehicle at the time of first explicit 
guidance steering commands are outlined in Table in-2. Pj_ is -80,407. 52ft 
because Stage 1 roll azimuth offset (implicit steering) has not been used in 
the run presented. Note the magnitude of difference between the vehicle at¬ 
titude and commanded attitude at this time. The magnitude of the out-of-plane 
position accounts for the large commanded yaw angle. The relatively large 
pitch was commanded because the flight path angle is low and the vehicle is 
low in the trajectory. Figures III-21 and -22 present the first pitch and yaw 
maneuvers of the Stage 2 steering equations. A fade-in technique was used 
in this run for approximately 20 sec.at the start of Stage 2 guidance. This 
sophistication is not in the present math flow (Figures E-1A and E-2A), and 
it is questionable whether something similar would be desired. If a trajectory 
of this type were actually experienced and a fade-in technique were not em¬ 
ployed, a 6 deg. attitude error signal would be delivered to the autopilot. Con¬ 
sideration should be given to those effects as well as those due to malfunctions 
and switchover during this time if fade-in were employed. Note that the com¬ 
bination of IGS limiting and K-factors in the autopilot will limit these rates to 
9 deg./sec. 

Figures III-22 and -23 present additional steering data over the entire 
Stage 2 flight. The various computed quantities were rather stable, consider¬ 
ing the magnitude of the initial changes required in vehicle orientation. 

Vehicle attitude error signals are contained in Figure IE-24. These 
signals are stable except near shutdown when they tend to deviate considerably 
as the equations near the pole at insertion. Although the attitude error data 
presented is smooth, Figure IE-27 presents a detailed picture of the construc¬ 
tion of the signals being delivered to the autopEot in this simulation. In all 
the data presented, the resolution of the commands in the simulation exceeds 
the resolution to be obtained in reading the platform gimbals and also exceeds 
the resolution to be expected from the output ladder networks of the computer. 
In order to better represent the physical system, plans are being made to in¬ 
clude the effects of platform gimbal angle uncertainty and resolution as well 
as the output ladder network resolution in further simulation efforts. 
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Velocity in ft./sec. Position in ft. 


Table III-1 


COMPARISON OF ENVIRONMENT & IGS 
NAVIGATION DATA DURING STAGE 2 FLIGHT 



Time of Flight 


t = 148.0 

t = 156.5 

t = 331.5 

t = 333. 15 

X E 

469,940.36 

543,112.46 

3,102,063.3 

3, 143,765.6 

X G 

469,938.52 

543,110.61 

3,102,137.7 

Not Computed 

^ y e 

-21,081,311 

-21, 102,675 

-21, 185, 532 

-21, 179,483 

a y g 

-21,081,672 

-21,103,057 

-21,186,309 

Not Computed 

1 Z E 

7,248.3934 

7,657.5726 

86,528.575 

86, 522. 584 

1 Z G 

7,764.1767 

8,231.3245 

88,933.193 

•Not Computed 

Re 

21,086,349 

21,109,664 

21,411,610 

21,411,708 

Rg 

Not Computed 

21,110,046 

21,412,399 

Not Computed 

X E 

8,429.5037 

8,789.8511 

25,065.476 

25,486.752 

o X G 

8,429.4739 

8,789.8409 

25,066.128 

Not Computed 

QJ . W 

to v 

\ X E 

-2,586.4419 

-2443.0887 

3602.0795 

3780.4620 

« y g 

-2,586.3138 

-2442.9486 

3602.3227 

Not Computed 

't Z E 

48.406679 

47.870196 

21.662948 

-29.290453 

§ Z G 

55.049632 

54.865414 

35.436428 

Not Computed 

> V E 

8581.6850 

9123.182 

25,322.98 

25,765.622 

V G 

Not Computed 

9123.174 

25,323.68 

Not Computed 


Eryj Environment; G<\, Guidance. 


Last computed value of time to go, Tq, at 331. 5 sec. equaled 
1.650838 sec. This would define shutdown at 333.150838 sec. 

R f = 21,411,962 ft. 

V F - 25,765.9 ft./sec. 
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Table III-2 


-VEHICLE INITIAL CONDITIONS- 

AT INITIATION OF STAGE 2 STEERING COMMANDS 
(t = 158.0 sec.) 


Vehicle Velocity, Vrj, 

9,160.475 

f PS 

Vehicle Velocity Perpendicular to 

44.651 

f PS 

Orbit Plane, V^, 

Vehicle Velocity Along Radius 

2,655.541 

f PS 

Vector, Vp, 

Vehicle Position Perpendicular 

-80,407.52 

ft. 


to Orbit Plane, Pj^, 

Flight Inertial Path Angle, T , 

Vehicle Pitch Attitude, 

Vehicle Commanded Pitch 
Attitude, /3"jj 

Vehicle Yaw Attitude, AY 

Vehicle Commanded Yaw Attitude,)8y 


0.2939 Radians 
0.3287 Radians 
0.4843 Radians 

-0.0014 Radians 
-0.1953 Radians 
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The discontinuity in the attitude errors (Figure III-27) every 10 compu¬ 
tation cycles is attributable to the introduction of new information from slow- 
loop calculations. Yaw rate terms are not used in the fast-loop equations; 
however, the pitch rate terms are included. For this reason, the discontinuity 
in the yaw channel may be expected to be greater than in the pitch channel, 
especially if heavy yaw steering exists. 

Figure III-28 illustrates the vehicle angle of attack during Stage 2 flight. 
This information corresponds well with the data presented in Figures IH-22 
and -23. Some of the intermediate computer constants are shown in Figure 
III-25. Note that the loss terms at the end of flight are 770 ft/sec.for gravity 
and 190 ft/sec.for steering. These quantities may also be on the high side 
because of the magnitude of steering during Stage 2. Computation of V\ f on 
Figure III-25 is not started until approximately 200 sec. following lift-off. 

For this reason, the transient appears in 8Tu at this time. 

Plots of various Stage 2 steering parameters are included as Figure 
III-26. Unfortunately, the data does not extend to SECO and had to be ex¬ 
trapolated. Note that the operational math flow is set up to do a fast-loop 
countdown on time-to-go, Tq, beginning approximately 2 sec prior to shut¬ 
down. At the same time, the vehicle attitude error signals will be set to 
zero so as to zero out the vehicle rates. The GEMINI IGS System similarly 
will not be computing these quantities during flight. In the simulation, the 
environment is still keeping track of the vehicle over this period of time. 

Some of the results at shutdown are included in Table III-l. 

Note that Figure III-26 indicates that R is approximately 650 ft greater 
than Rjr, and the environment indicates Ru approximately 210 ft less. The 
difference in results can be explained by the Y-axis roundoff errors previ¬ 
ously discussed. However, the IGS equations should have done better in 
making insertion altitude. Several reasons may be postulated, the most sig¬ 
nificant of which may be the bypassing of the pitch and yaw rate, equations 
following 300 sec. Other sources such as differences betweenand P may 
contribute to the difference in results. This area requires more detailed in¬ 
vestigation. 

The quantities V v , V p , Vj_ in Figure m-26 are all nominally small at 
or near insertion. Table III-l would tend to indicate that vehicle velocity is 
within 0. 3 ft/sec of Vjr. 
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Figure III-8. Pitch Plane Angle of Attack (o^ Run B and Run C 
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Figure III-9. Yaw Plane Angle of 
Attack (ay) Run B and Run C 
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Figure III-12. Engine Angles and Aerodynamic Forces 
Results of Run C, Combined Simulation 
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Figure III-17. Down Range Position (Rx) 
and Difference Terms 
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Figure III-18. Radial Position (Ry) and Difference Te: 
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Figure III-21. Stage 2 Steering 
Commands (At Start of Stage 2 Guidance) 
















































































































































































































































































































































































































































































































































































































































































































Figure III-25. Stage 2 - Intermediate 
Computer Constants and Velocity 
Loss Terms 
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Figure m-26. Stage 2 Steering Parameters 
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Figure III-27. Stage 2 IGS Attitude 
Errors (Magnified To Emphasize 
The Need To Include Quantization In 
The Simulation 
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Figure III-28. Vehicle Angle of Attack, Velocity, 
and Radius to Earth Center During Stage 2 Flight 
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Section IV 


PLANS FOR CONTINUED SIMULATION STUDY 

A. OBJECTIVES 


The primary objective of future ascent guidance studies is to obtain an 
over-all error model for the IGS ascent guidance system. The secondary ob¬ 
jective is to study switch-over transients. In addition to these objectives, the 
various areas which required additional investigation and simulation will be 
considered. * Studies which are being planned are discussed in this section. 

B. FIXED POINT PROGRAM 


The simulation is programmed in 27 bit floating point binary. This im¬ 
plies that in many areas accuracy will exceed the capability of the 26-bit, 
fixed-point computations to be done in the actual GEMINI IGS. Errors due to 
this difference are difficult to assess; therefore, a fixed-point simulation of 
the IGS equations is now being completed. Plans are also being made to tie 
this system to the programs developed for this simulation to assess some of 
these errors. 

C. EFFECTS OF UPDATES ON STAGE 2 EQUATION 


The late updates during flight may tend to introduce slight discontinuities 
in the velocity data used in the steering equations. These effects along with 
the effects of poor updates on the solution to the steering equation will require 
evaluation. The total uncertainty in the update data will also be considered so 
that the accuracy of the required update data can be assessed. 

D. INITIALIZATION OF GIMBAL ANGLES PRIOR TO LAUNCH 


The reading of the platform gimbal angles prior to platform release will 
be considered along with the setting of these angles equal to the commanded 
angles following platform release, because the platform azimuth could be mis¬ 
aligned by as much as 45 min. This misalignment could result in a 45 min. 
attitude error as long as initial platform gimbal angles are obtained from 


♦Scheduling is not included in this discussion. Information on scheduling will 
be included in "Plans for Ascent Guidance Study," Feb. 1963, (to be published). 
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preset and computed values. Similarly, the initial commanded gimbal angles 
could be in error as much as 14 min. because of vehicle, spacecraft, and IMU 
misalignment as well as gimbal reading uncertainties in the platform. Since 
the resolution in gimbal reading and readout noise is much less than 14 min., 
accuracy can be gained by initializing as described. 

E. LOADING DATA OR CONSTANTS INTO THE COMPUTER 


From the operational standpoint, it may be desirable to load a number 
of constants into the computer within several weeks or possibly hours of a 
flight, e. g , Rf, Vf or capability to remove position steering computations. 

A thorough review of the equations and constants will be made to determine 
if this action is beneficial. Constants or variables requiring frequent or oc¬ 
casional changing will also be outlined. Procedures to be used to provide 
this update capability will have to be prepared should it become necessary. 

F. GROUND CHECK-OUT AND SIMULATION OF ASCENT PROGRAM 

During launch pad check-out, a capability must be provided to opera¬ 
tionally exercise the ascent program. This provision may be difficult with¬ 
out additional instructions in the program or some sophistication in the AGE 
equipment. The major difficulty is in the Stage 2 equations where the nomi¬ 
nal solution depends on an exponential thrust profile, and where the delivery 
of a SECO signal is required. A check-out procedure is being studied withr 
out additional instructions. However, additional instructions will be con¬ 
sidered if required. 

G. GIMBAL READING UNCERTAINTIES 

Various errors can be introduced through the platform gimbal reading 
due to misalignment of the IMU, spacecraft, and GLV, or through the uncer¬ 
tainties in the gyro reading. In fact, the gimbal reading uncertainty could be 
as much as 10 min. A model of the gimbal reading resolution, uncertainty, 
and possibly noise will be introduced in the simulation to estimate these ef¬ 
fects. 

H. OUTPUT LADDER NETWORK SIMULATION 

Figure III-27 indicates that the resolution of the IGS output ladder net¬ 
work is quite large compared to the accuracy to which the attitude errors 
have been computed and used in these programs. Consequently, these effects 
must be included in the simulation to assess their effect on vehicle stability 
and insertion error. 
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TIMING OF DISCRETES 


Timing of the discretes in the Ascent program is generally controlled 
in the slow-computation loop. The accuracy of these discretes can be con¬ 
trolled better if approximately one-half computation cycle was subtracted from 
each nominal discrete time. For example, the present accuracy of a discrete 
is * NOMINAL +°- 5 > -0. 0 sec. The above operation, however, would provide 
an accuracy of tfjOMINAL ±0. 25 sec. 

J. DIFFERENCE IN R AND Rf AT SHUTDOWN 

Reasons for the large discrepancy in R-Rf at insertion will require a 
more detailed study. The effects of extending position-steering computations 
closer to shutdown will be investigated in an attempt to improve these accur¬ 
acies. Judgment, however, must be exercised because the payload penalties 
to do position steering late in flight could exceed the advantages gained. 

K. FADE-IN OF INITIAL STAGE 2 ATTITUDE ERRORS 

As previously discussed, the fade-in of initial Stage 2 attitude errors 
should be considered to provide a smooth transition from Stage 1 to Stage 2 
flight. Advantages and disadvantages will be weighed before any decision is 
made. 

L. PROVISION FOR STUDYING SWITCH-OVER TRANSIENT 

Since the IGS is a back-up system to the primary Radio Guidance Sys¬ 
tem to be used on GEMINI, it will be necessary to include the investigation of 
switch-over transient in the simulation. The following methods are being con¬ 
sidered for this investigation: 

1. Include radio guidance system in the simulation which will 
allow switch-over at any time in the trajectory. 

2. Include capability in the present program which will allow 
entry into the ascent profile at any time along the trajec¬ 
tory with selected initial conditions. 

M. ATTITUDE ERROR QUANTIZING 


The attitude errors supplied to the autopilot from the guidance equations 
in this simulation are accurate to eight significant digits. In the actual 
GEMINI computer, these outputs will be supplied on the output ladder network 
where they are quantized to 0.12 deg. The effects of this quantization are dis¬ 
cussed under Stage 2 results in Section III-C. This quantization will be incor¬ 
porated in future computer runs to determine its effect on system stability. 
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N. ORBIT VELOCITY ADJUST 


Simulation efforts to date have not included the orbit velocity adjust 
equations which are a part of the operational math flow (Figure E-IB). Two 
plans are being considered for the simulation of these equations. 

The first plan is to use the present simulation and approximate the con¬ 
trol system with a perfect autopilot. This might be pictured as a third Stage 
in the GEMINI ascent. An approach of this type would certainly provide an 
indication of the adequacy of the equations; however, parametric analysis 
would be somewhat difficult. 

The second approach would use an entirely separate program. A forcing 
function would be provided and a perfect autopilot assumed. This would pro¬ 
vide maximum flexibility in the study of these equations. 

O. PROGRAM COMPATIBILITY PLANS 


Specific items will be modified to reach a state of compatibility with 
other GEMINI studies regarding data used in the program. These changes 
will facilitate comparisons of results which might be available through the 
Martin, Aerospace, and McDonnell simulation efforts. The items are as 
follows: 

1. Write and substitute into the program a subroutine de¬ 
scribing gravitational acceleration according to approved 
GEMINI models. 

2. Write and substitute into the program a subroutine de¬ 
scribing atmospheric coefficients, based on NASA TND- 
595 "A Reference Atmosphere for Patrick AFB. " At this 
time surface winds will also be introduced. 

3. Update vehicle characteristics to be compatible with doc¬ 
ument LV-140. Specific characteristics are aerodynamic 
coefficients and location of center of pressure. 

4. Provide for inclusion of engine pitch angle offset. 

5. Revise weight, thrust, and mass flow rates to latest com¬ 
patibility data. 

6. Provide correction to account for initial orientation of the 
platform with respect to the orbit plane as described in 
Appendix A. 
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P. MATH FLOW CHANGES 


The following list outlines areas in the math flow where changes have 
been incorporated or are currently being planned: 

1. Add Third Pitch Step 

2. Add Low-Level Sensor Discrete 

3. Modify Test onAQ (A205) 

4. Change Standard Update Increment to 40 sec. (A 128) 

5. Move platform release to engine ignition and delete 
ephemeris calculations. 

6. Set commanded platform gimbal angles equal to actual 
platform gimbal angles prior to platform release 

7. Modify Orbit Insertion Equations 

8. Modify SECO sequence. 

These changes will be added to the simulation as they are incorporated 
in the operational computer. Additional details concerning the above changes 
will be contained in the math flow change records as they appear. 
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Appendix A 


EXPLANATION AND DERIVATION 
OF GEMINI ASCENT GUIDANCE SIMULATION 

I INTRODUCTION 


This appendix defines the composite approach used initially in IBM 
studies on GEMINI ascent guidance simulation. The symbols and mathematics 
will not necessarily follow exactly in all programs and subroutines, but each 
interface will be compatible. 

The task of writing the main environmental program was simplified by 
using existing vehicle guidance simulation programs where possible. How¬ 
ever, coordinate systems and symbols were not always identical. Since the 
specific areas in which the differences exist are relatively unimportant, 
most of the differences will not be identified. One exception is the case of 
the coordinate systems chosen for programming. 

Certain simplifications have been made in preparing this simulation. 
These can be classified under one of two categories. One category covers 
the use of existing data rather than the rewriting of subroutines. The other 
covers simplifications which were made by neglecting elements of the prob¬ 
lem that were not under investigation or would not contribute to the final 
results. 

The first category, which covers gravitational computations, deter¬ 
mination of atmospheric coefficients, corrections for initial orientation of 
the platform, use of constant thrust, simplification of platform orientation 
and omission of zero offset of engine gimbal, are discussed in this report. 

In the second category the following simplifications were made: 

• The vehicle was considered to be a rigid body with no bending or 
fuel slosh modes. 

• All body rates and angles were assumed to be measured around 
the center of gravity. 

• Body rotational rates were considered constant during one 
computation cycle. 
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• Vehicle position and velocity were approximated by a quadratic 
for one complete cycle. 

• Only first-order aerodynamic forces and moments were con¬ 
sidered. 

• Instrumentation errors were not considered. 

II LIST OF SYMBOLS 

A. COORDINATE SYSTEMS 

X E> Y E> Z E - Earth referenced frame (inertial) 

*Z* " O r bit plane reference frame (inertial) 

I x , Iy, I 2 " Platform reference frame (inertial) 

*xb’ ^b’ l zb ~ s P acecraft referenced frame 
I x , V *z -Vehicle reference frame 

B. MATRICES 

a^ j - Earth to platform transformation matrix 
£ bjj j - Platform to spacecraft transformation matrix 

C. ALPHABETIC 

t - Time 

t pr - Time of platform release 

i - Inclination of desired orbit plane 

C Q - Magnitude of eastward inertial velocity at launch point due to 

earths rotation 

C 5 - Angle from east to the spacecraft I z ^ axis 

X s - Station number of engine gimbal 

Zj, Z 2 - Distance from vehicle centerline to engines A and B 
respectively 
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X C p - Station of the center of pressure of the vehicle 

X C g - Station of the center of gravity of the vehicle 

W - Vehicle weight 

M - Vehicle mass 

q - Aerodynamic pressure 

Sir - Vehicle frontal area 

Cg - Speed of sound 

M - Mach number 

I x . Iy, I z - Moments of inertia about I x , I y , and I z respectively 

(note that distinction in symbols must be made by reader 
where necessary) 

G 0 - Value of gravitational acceleration for conversion of fuel 
weight to mass 

- Axial force coefficient 

Cy - Side force coefficient 

Cn - Normal force coefficient 

F - Force 

subscripts X, Y, Z - Along platform axes 

x, y, z - Along vehicle axes 

xb, yb, zb - Along spacecraft axes 

A x , A , A z - Aerodynamic in vehicle frame 

Axb> A b; A 2)3 - Aerodynamic in spacecraft 
frame 


R - Position vector 


subscripts 

X, Y, Z - In platform frame 

X E> y e j Z E - In earth-referenced frame 
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V 


- Velocity 


D. 


subscripts 

X, Y, Z - Platform frame 

X E , Y e , Ze - Earth referenced frame 

w X» W Y> W Z " ° f wind in P latform f rame 
w x , w y , w z - Of wind in vehicle frame 

w xb> w yb , w zb - Of wind in spacecraft frame 


GREEK LETTERS 
Q - Platform pitch angle 

<J/ - Platform yaw angle 

cf) - Platform roll angle 

\ - Launch point latitude 

$ - Launch point longitude 

£l - Longitude of the ascending node of desired orbit plane 
8 - Angle between local geocentric vertical and orbit plane 

y - Angle from east to platform lx axis 

4> 0 - Initial roll angle of platform 

r - Integration time period 

€ 1) € 2 " Vehicle engine yaw angles 

8 i, 82 - Vehicle engine pitch angles 

oj - Earth's rotational rate 

e 

ou x> w yj t o z - Vehicle rotational rates 
ojxb, w y b > w zb " Spacecraft rotational rates 
aN - Pitch angle of attack 

q Y - Yaw angle of attack 

p - Air density 
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Ill COORDINATE SYSTEMS 


A. SYSTEM Xg, Y e , Z e 

This system (see Figure A-l)_is an inertial reference with the Z~ e axis 
through the north polar axis. The Xe axisjs determined by the position of 
the Greenwich Meridian at the time tp r . Y e forms a right-handed orthogonal 
set. 


Z 


E 



Figure A-l 

Data in this system is required to initially determine and continuously 
provide geographic latitude which is equal to inertial latitude. Geographic 
longitude does not have to be determined because the earth is considered to 
be symmetrical about the equatorial plane. However, it could be determined 
by rotating the system at earth's rate. 



The following angles are measured in this system: 

• ft, the longitude of the ascending node of the desired orbit plane 
at time t = tp r . For positive angles,ft is measured as a right- 
handed rotation about Z e . 

• i, the inclination angle of the desired orbit plane. For positive 
angles,_i is measured as a right handed rotation about the dis¬ 
placed X e axis. 

B. ORBIT PLANE SYSTEM T x *, Ty*, I z * 

This system is defined primarily by the normal to the orbU plane, 

Tz*. I x * is parallel to I x of the platform reference frame, and Iy* forms a 
right-handed orthogonal set. 

Figure A-l shows the relationship between system Xe, Ye, Ze and the 
orbit plane. 

C. PLATFORM REFERENCE SYSTEM T x , Ty, I z 

This system may not be the true platform reference. The frame is 
inertial after launch; but, due to platform drift, the actual platform which is 
also nominally inertial after launch may not coincide with this system. For 
the purposes of this simulation platform, drifts and initial misalignments 
will be ignored. 

This system is defined with theT x axis parallel to the desired orbit 
plane and horizontal with respect to earth at time of launch. TheTy axis 
is along the geocentric radius vector towards the center of the earth. Iz 
forms a right-handed orthogonal system. In practice, this system will be 
aligned to the local gravity vector. There is a basic simplification in this 
definition; see Section IV for explanation. 

D. SPACECRAFT REFERENCE SYSTEM I xb , l yb , T zb 

This system is_defined as the x b , y b , z b system used in the GEMINI 
capsule reference. I xb is aligned with the longitudinal axis of the space¬ 
craft; the positive direction is towards the nose of the capsule. I zb is aligned 
to the vehicle yaw axis; the positive direction is down when referenced to the 
capsule during normal orbital flight attitude. I yb is along the vehicle trans¬ 
verse axis and forms a right-handed orthogonal system. 
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E. LAUNCH VEHICLE SYSTEM I x , I y , I z 

The capsule is mounted on the launch vehicle_so there is_a 90 jieg. 
rotation_about the longitudinal axis of the capsule, I x b- Thus, I x b = I x , and 
Iyk and I z b are rotated 90 deg. fromTy andl z , respectively. This is a 
principal axis system located at the vehicle’s center of gravity. 

For the remainder of this Appendix, vector notation will be omitted 
although implied. 


IV COORDINATE TRANSFORMATIONS 


y Throughout the simulation, two coordinate transformation matrices, 
and [b..] , will be used extensively. The matrix [a^l defines the 
onship between the X e , Y e , Z^ system and the Ix, l Iy, *Z system. The 
matrix [b^l defines the relationship between the I x b, I y b, fzb system and the 
I X , Iy, Iz system. Thus, both earth and body angles are known with respect 
to the platform. As previously noted, both Ix, IY, Iz and x e> Y e , Z e are 
inertially fixed; therefore, [ ajj ] is a constant. Since the derivation of 
jayj leads to other subjects, it is convenient to discuss [by] first. 

A. DERIVATION OF MATRIX [bjj] 

The matrix [ bjj 1 is the transformation from capsule coordinates to 
platform coordinates. J 



Given: a) Two right-handed coor¬ 
dinate systems I x b, Iyb> 

I z b and Ix, Iy> aligned 
in a zeroed position (see 
Figure A-2). 

b) 0, the angle between plat¬ 
form and innermost gimbal 
(1st platform rotation). 

c) the angle between inner 
gimbal and outer gimbal 
(2nd platform rotation). 

d) 4>, the angle between outer 
gimbal and body (3rd plat¬ 
form rotation). 
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1 . 


First Platform Rotation - About Z Through 6 



'/'is positive with X" to the right of X'Y' plane. 
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Now, 


X" 

’ = l xb 

X = r x 

Y" 

' = Izb 

y = i y 

Z" 

' = - Iyb 

z = Iz 


Ixb 


cos'/'cos 8 -cos^sinfl -sinf 


ix 

Iyb 

- 

-cosf?cos^sini//+sin</>sin0 cosi^sint/rsind+sin^cos# -cos<£cos iff 


iy 

Izb 


sin<f>simf/cos8+sin8cos^> -sin8sin<f>sinf+cos<f>cos9 cos^sin <f> 


iz 


B. RELATIONSHIP BETWEEN I x b, Iyb, I Z b and I x , I y , I z 


Ix - Ixb 



! yb 


Figure A-3 shows the relation¬ 
ship of the capsule and launch vehicle 
systems. This configuration results 
from having the astronaut's head 
pointed towards the south with Iyb 
nominally eastward. During the 
pitching maneuver, vehicle pitch is 
negative around Iy. 


The following equations, based on Figure A-3, will be frequently 
employed throughout the report when relating launch vehicle and spacecraft 
quantities: 

*x = *xb 
Iy = _I zb 
*z = Iyb 
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C. DERIVATION OF MATRIX [ay] 

Recall that the initial platform orientation depends on the desired orbit 
plane. Some relationships will be established between the orbit plane, earth 
reference system and platform (Figure A-4). 

The following derivations are made whidh relate orbital, geocentric 
and platform quantities: 

1. Derivation of unit vector normal to orbit plane (Ig*) in terms of 
earth-referenced axes,XIand i. 

2. Determination of unit vector Iy in terms of earth-referenced 
axes, \and<£\ 

3. Determination of angle 8 between Iy and Iy* - - this angle is 
considered positive when orbit plane is above platform plane at 
launch site. 

4. Determination of unit vector lx = Ix*- 

5. Determination of angle / between east and I x axis - this angle is 
considered positive when I x is north of east. 

6 . Determination of matrix |ayj . 

The normal to the orbit plane can be derived by considering two 
rotations from the earth-referenced system: (1) A rotation about Zg through 
the angle SI, (2) about the displaced Xe axis through the angle i. The third 
rotation is not required since it would be about the displaced Ze axis. 



Figure A-4 
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1 st Rotation ft About Z, 



2 nd Rotation i About X A 



x a = x a 


From the above 

Z A ' = sin i sinft Xg - sin i cos ft Ye + cos i Zg 

and, 

Za = I Z * 

Thus, the normal to the orbit plane is given by 

I z * = Ixe sin 1 sin ^ - IyE sin * cos ft + I 2 E cos i (1) 

Since the platform Iy axis has been defined as along the launch point 
radius vector, it can be determined by two rotations. These rotations are 
about Ze through^), and then about the displaced Xg axis through X. 
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1. First Platform Rotation -- <ft' About 



From the above, 

X" = cos X cos <f> 1 Xg + cos X sin <£' Yg + sin X Zg 


By previous definition, 

X" = -I Y 
Therefore, 

-Iy = Ixe cos ^ cos<f>' + IyE cos X sin<£' + Izg sinX (2) 
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Now since lx and lx* have been defined as equal vectors, Iy and Ig* must be 
in a common plane. The cosine of the angle between Iy and Ig* is equal to 
the sine of the desired angle 8 and can be determined by: 



And, 

Sin 8 = sin i cos X sin (<f>' -ft) - cos i sinX (3) 

The vector lx can be evaluated in terms of the above quantities by 
evaluating the vector cross-product of Ig* and - lx as follows: 

I z * x - Iy - I-jj- sin (——n 8 ) = ^x 8 
2 

X XE l YE 

lx cos 8 = sin ft sin i -sin i cos ft 
cos X cos <t>' cos X sin<£' 

lx cos 8 = Ixe [ -sin i cos sin ^ - cos ^ sin<£' cos i j 

+ Iyg [cos i cos X cos <f>' - sinft sin i sinxj 

+ I Z e [ sinft sin i cos X sin<f> ' + sin i cosftcos Xcos<£* j (4) 

This vector lx is also defined by the third rotation which relates platform 
to earth values. This rotation is about X" through the angle y . 


J ZE 
cos i 
sin X 
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- x ,,r 


“cos X cos<£' 


cos Xsin<£' 

sin X ” 

Y’ M 

= 

-siny sin X cos <f>' 

- cos ysin<£ 

-siny sin Xsin<^ + cosycos<^ 

siny cos X 

Z”’ 


-cosysin Xcos<£' + 

sinysin<£’ 

-cosy sin Xsin^>'- sinycos<p 

cosy cos X 


By previous definition, 

X'” = -I Y 
Y"' = I X 


Z ,M = I Z (5) 

Therefore, 

l x = I XE [-sin / sin X cos<£' - cos y sin <f>' ] 

+I Y e [-sin y sin X sin<£’+ cos y cos<£’j 

+i ze [ sin y cos X 


( 6 ) 



From equations (4) and ( 6 ), 


lx - ^XE ^ + *YE ® + *ZE s ^ n ^ cos ^ 
lx cos 8 = Ixe C + IyE D + ^ZE |^sin&sin 
By equating Ize terms: 


i cos X sirup + sin i cosTi cos X cos< 


1 

sin y cos X =-Tsinfi sin i cos Xsin<£' + sin i cosftcosX cos^'l 

cos 8 ^ -1 

This expression may be reduced to determine the sine of the angle between 
Y M (east) and lx 

sin y = sin i cos (<f>' -fl) (cosS) - * ( 7 ) 

The matrix fa^l is defined by equations (4) and (5). 


-cosysin<£'-sinysin Xcos <f>' 

cosycos<£'-sinysinXsin<£' 

sin ycos X 

J XE 

-cos\cos<p 

-cosX sirup 

-sin X 

r YE 

- cosy sin Xcosc£’+sinysin</>' 

-cos ysinX sin<p-sin ycos<£' 

cosy cos X 

J ZE 


V 


X XE 

I Y 

■ M 

J YE 

l Z 


_ r ZE_ 


4. Difference Between Simulation and Actual Platform Alignment 

A basic relationship was neglected in the derivation of the equations 
above to define the relative orientations of the several coordinate systems. 
In the physical application of the guidance equations, the inertial platform 
will be aligned such that the Y-axis is along the local gravity vector. The 
error associated with this simplification is, as expected, primarily noticed 
in out-of-phase position and velocity computations. This discussion defines 
the changes that will be required in order to simulate the IGS platform more 
closely. 
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The figure shows the actual relationship that will exist between the 
referenced inertial frame and the actual platform. 



Thus, the actual platform frame can be described in terms of one 
rotation, 8 ', about X. 8 'is considered positive as shown, and is the angle 
between the local gravity vector and the geocentric radius vector. This 
relationship may be expressed by the following equations: 


Xp 


1 0 

0 


X 

Yp 

= 

0 cos 8* 

-sin 8 


Y 

z p 


0 sin 8 

cos 8* 


Z 

- - 


_ 

_ 


- 


Let sin 8' = * y, then 


~ Xp" 


1 0 o'" 


~ X ' 

Yp 

= 

01 -«y 


Y 

_ zp_ 


0 « y 1 


Z 


The following assumptions are implicit in this relationship: 

1. cos 8 ' = 1 -- This is valid since 8 ' < 15 min 

2. Xp = X -- Again, since S' is negligible, y may be assumed equal 
to the actual rotation of the platform in azimuth, (error is less 
than 30 sec) 
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The relationship between the orbit frame and the actual platform may 
be established by neglecting products € y Vy (less than 10 - 5) ; 


r x*~ 


1 

o 

o 


Xp' 

Y * 

= 

0 1 . ‘y-’y 


Yp 

z* 


_0 Vy - « y 1 


_Zp_ 


Consider the initial conditions inserted into the computer (Xc = 0, 
Yc = Rn, %c ~ 0)> which are required to satisfy the gravity equation 
(G = K/R2), then the following relationship are valid: 


X 

= 0, Y = 

Ro, 

Z = 0 


Xpo 

= 0 , Yp 0 = Rq, Z Po - 

€ y R o 

Xp 

= AXp + 

Xpo 

= AXp 


Yp 

=A Yp + 

Ypo 

= AY P + 

Ro 

Zp 

=A Zp 1 + 

Zp 0 

= AZp + 

«yRo 


Where AXp, AYp, AZp are the changes in Xp, Yp, Zp from time t 
to time t. Now the values in the computer can be identified: 1 

X C =AX P , Y c =AYp + R 0 , Z c = AZp 

Thus, 

X* = X C 

Y* = Y C + («y -Vy) Z C +fy («y -Vy) Rq 
Z* = Z C+ € yRo + (^y- e y) Y C 

This will be implemented in Stage 2 guidance equations in the following 
manner: 


R Z * = Z + * y R 0 + (v y - < y ) Y 
Vj_ =i + (V y -* y )Y 

Changes in the simulation will be (1) establish the matrix between the 
inertial reference and the platform frame, and (2) use this matrix to trans¬ 
form position, velocity and accelerations from the inertial reference to the 
platform frame. 
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V. INITIAL MATRIX y 

Figure A-6 shows the relationship between the platform, capsule and 
launch vehicle prior to launch. We assume that the platform is aligned to the 
orbit plane and earth reference system as previously noted. 


*zb 


!xb 



In Figure A-6, Iy 
points toward the earth's 
center; I x and I x b point 
away from the earth’s 
center (vertical). All other 
vectors are in a plane nor¬ 
mal to Iy. 

Y = angle from east 

to I X 

C5 Wangle from east 
to I z b 

4 >o = angle from Ix 
t0 ^b 


Figure A-6 

The following matrix may be written from Figure A-6. 


_ Ixb“ 


0 

-i 0 

I X ' 

r yb 

= 

sin</> 0 

0 -cos 4> 0 

iy 

Izb 


COS^q 

0 sin<#> 0 

_ l Z 


In matrix by j , 

bi3 = -sin<// = 0 
bn = cos\J/ cos 9 = cos# = 0 
bj2 - -cos \Jr sin 9 = -sin# = -1 
t>33 = cosvj/ sin<f> = sin<£ = sin<£ 0 
t»23 = -cos^ cos <f> = -cos<(> = -cos<t > 0 
From the above, 

\f/ = 0 
# = 90° 

4> = <t> Q = C 5 - Y + 90° 
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VI. DESCRIPTION OF LAUNCH VEHICLE MODEL (FIGURE A- 7 ) 


The six parameters described in this section have been included in the 
description of the vehicle. Other parameters not defined have been omitted 
because (1) we want a simplified model and (2) they do not contribute 
appreciably to an analysis of navigation equations. 

A. REFERENCE SYSTEM 

The capsule, vehicle and platform systems have been defined in 
Section III, Coordinate Systems. 

B. THRUST APPLICATION 

Thrust is assumed to be applied at the engine gimbal. During Stage 1 , 
two engines are in operation. They are symmetrically located about the 
vehicle longitudinal axis in the I x , Iz plane. The longitudinal station of the 
engines is denoted asX s ; the distance from the vehicle axis to the engines is 
Zj to engine A andZ2 to engine B. During Stage 2 a single engine is used. 

This engine is located on the vehicle longitudinal axis. All engines are 
gimballed in 2 -degrees-of-freedom, e and 8. Positive directions of these 
angles are shown in Figure A- 9 . 

C. AERODYNAMIC COEFFICIENTS 

The three coefficients, axial C^, side Cy and normal Cn, have been 
determined from preliminary GEMINI Launch Vehicle (GLV) wind tunnel data. 
Fixed data points, which show the coefficients versus mach number, are 
stored for various angles of attack, a n and ay. Linear interpolation is used 
to determine the required values. The angles of attack are defined by Fig¬ 
ure A- 7 . Aerodynamic forces act through the center of pressure X^. The 
location of X cp is also defined on the basis of mach number and angle of 
attack in a stored table. 

D. MOMENTS OF INERTIA 

The moments of inertia -- roll I x , pitch I y and yaw I z about the center 
of gravity -- have been determined from preliminary GLV data sheets. These 
values are stored as polynomials in the computer and determined as functions 
of time. Location of the center of gravity, X cg , is also stored as a polynomial 
in time. 
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E. VEHICLE WEIGHT 

The vehicle weight is a linear function of time determined by initial 
weight and fuel flow rate which were determined from preliminary data 
sheets. Mass is determined by converting weight through use of standard 
gravity value, G Q . 

F. THRUST 

Thrust is considered constant at the values specified in Appendix F. 



Figure A- 7 . Vehicle Model Used in the Simulation 


VII DETERMINATION OF AERODYNAMIC FORCES 
A. RELATIVE WIND 

Relative wind is defined asUie velocity of the vehicle with respect to 
the air mass. Vectorally, V w = V - V a . P'or initial studies, we assume 
that the air mass has zero velocity with respect to earth. Underjhis 
assumption, the inertial velocity of the air mass is given byw e xR, where 
w e is earth's rotation and R is the position vectorof the vehicle with respect to 
the earth's center. 
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Since w e can be defined as a vector in inertial space, this vector can be 
transformed into platform coordinates through matrix [a^ ]. The position 
vector, R, is updated continuously and is available in the platform frame. 

Thus, 


V e platform w e platform x R platform 
w e platform = [ a ij] w e = [ a ijj] | w e | Z E 

= a 13 "e l X + a 23 w e *Y + a 33 w e r Z 
= w eX r X + w eY I Y +w eZ *Z 
R"platform = Rx *X + R Y r Y + Rz *Z 

and 

Ve platform = V e x IX + V e y Iy + V e z Iz, 

where 


V e X = w eY Rz - w eZ Ry 
v eY = w eZ RX - w eX R Z 

v e z ^exRy-VRx 

and 

V w = V - V e 
V W X = V X - V eX 
V W Y = V Y - V eY 
V wZ = v z - v eZ 

Now, the relative wind velocity must be transformed into vehicle 
coordinates (Figure A-8): 

V w spacecraft = [ b ij] v w p i a tf 0 rm 
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Explicitly, 


v wxb = b H v wX+ b 12 v wY + b 13 v wZ 
v wyb = b 21 v wX+ b 22 V wY + b 2 3 V wZ 
v wzb = b 31 v wX+ b 32 V W Y + b 33 V w z 
Convert into vehicle coordinates: 
v wx = v wxb 
V W y = -V wz b 
^wz = V W yb 

The magnitude of wind velocity, V w , is foundry root-sum squaring the com¬ 
ponents in spacecraft coordinates. 

v w = v wxb^ + V W yb^ + v wzb^ 

B. ATMOSPHERIC COEFFICIENTS 

The air density, P , and speed of sound, C s , are determined by using a 
previously developed IBM 7090 subroutine, based on the 1959 ARDC Standard 
Atmosphere. Future plans provide for incorporation of data from NASA TND 
595 Reference Atmosphere for Patrick AFB. At this time, surface wind 
conditions may also be introduced. 

C. COMPUTATION OF AERODYNAMIC FORCES 

Three components of aerodynamic forces -- axial F^x> side F Ay and 
normal F^z -- are defined as follows: 

f Ax = "Qa 9 S it 

F Ay = " C Y 9 S 7T 

f Az = " C N q s ir, 
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Note that, 


and 


and 


a y = tan " 1 


_v wzb 


v wxb 


a 


N = tan 


-1 


v wyb 

v wxb 


F Axb = f Ax 
F A yb = f Az 
F A zb = " F Ay 


VIII. THRUST RESOLUTION 

As has been previously mentioned, the engines are gimballed in two 
degrees of freedom. This is shown in Figure A-9 along with the resolution 
of thrust into vehicle coordinates. 



(Note) forces are defined in a positive 
direction. Therefore, a 
positive 8 angle yields a negative 
force. 


Figure A-9 
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T Ax« T A COS&l COS€ l 
T Ay - t A sin «l 
t Az = -T A sin Si 


t Bx » t B cos 82 cos €2 
TBy = t B sin « 2 
t Bz = -Tb sin 82 


IX. RESOLUTION OF PERTURBATIVE FORCES 

Forces in vehicle coordinates may be summed as follows: 

F x = f Ax + t Ax + t Bx 
F y = F Ay + T Ay + T By 
F z = f Az + t Az + t Bz 

By transforming these forces through capsule coordinates and into the platform 
frame, the following equations are obtained: 

F platform = £ b ijj * F capsule 

F X = b ll F x + b 21 F z -b 31 F y 

F y = bi 2 F x + b 2 2 F z " b 32 F y 

F Z = b 13 F x + b 23 F z _b 33 F y 

X. CENTRAL FORCE FIELD 


Gravitational acceleration is determined by using a previously written 
subroutine which includes the first two oblateness teams. As previously 
specified, position of the vehicle may be transformed into an earth-referenced 
system, Xe, Ye, Z e , through [ay] . After computation of gravity in this 
system, G is transformed back into platform coordinates. 

R earth = [ a ij] 

G platform = [ a ij] 


R platform 
G earth 
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XI UPDATING POSITION AND VELOCITY 


The total acceleration acting on the vehicle may be determined from 
the equations in Sections EX and X: 

— F _ 

V =— + G 
M 


or, in terms of platform coordinates, 


F X 

v x =- 

M 

+ Gx 

Fy 

Vy=- 

M 

+ Gy 

F z 

v z =- 

M 

■ + G Z 


Velocity and position are updated by numerical integration by using 
Simpson's Rule. 




XII. VEHICLE ATTITUDE AND ANGULAR ROTATION 

Two cases are discussed in this section. The first case is a perfect 
autopilot, i.e., the attitude is governed solely by the guidance equations. 
The second is the tie-in of an autopilot model coupled with solving moment 
equations to determine angular rotation. 
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A. 


"PERFECT" AUTOPILOT (ATTITUDE GOVERNED BY GUIDANCE 
EQUATIONS) 


This simulation is designed to work in conjunction with simulated inputs 
determined by IBM guidance equations.* The primary outputs of these 
equations are vehicle errors e x , € y , « z about the three vehicle axes. Since 
the computation cycle is short and the expected values of the orientation 
errors are small, these errors are considered commutative, i.e., order 
of rotation about vehicle axes is not important. In this case, the following 
derivation demonstrates the method used to update the matrix JbjjJ . 

Given: 


Ixb 


b ll b 12 b 13 


ix 

!yb 

= 

b 21 b 22 b 23 


! y 

*zb 


b 31 b 32 b 33 


l z 


Differentiate with respect to time in X Y Z frame: 



From the theorem of Coriolis, 


dR 

dR 


— 

= — 

+ w 

dt 

dt 


XYZ 

xyz 


R 


(3) 


♦Appendix E 
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and, letting R = I x + Iy + I z , 
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Equate (6) and (2): 



The nine differential equations for bjj can now be written: 


b ll 

= w zb 

b 21 

w yb 

b 31 

b 12 

= "zb 

b 22 

w yb 

b 32 

b 13 

= w zb 

b 23 

w yb 

b 33 

b 21 

= w xb 

b 31 

w zb 

b ll 

b 22 

= w xb 

b 32 

w zb 

b 12 

b 23 

31 w xb 

b 33 

w zb 

b 13 

b 31 

- yb 

b ll 

w , 
xb 

b 21 

b 32 

= w yb 

b 12 * 

w xb 

b 22 

b 33 

= w yb 

b 13 ' 

w zb 

b 23 


Use the first two terms of a Taylor series expansion; then 

b ij t = b ij t -l + b ij At 

Let w At = e and simplify: 


b nt 


, + € zb b 2i 
-1 i-1 

' € yb b 31. 

b 12 i 

= b 12j 

_ 1 + € zb b 22 i l - 

’ e yb b 32j 

b 13 i 

73 b 13i 

+ 6 zb *>23^ ' 

’ € yb b 33j 

b 21t 

= b 21: 

i + € xb b 31. , ‘ 

• *zb b ll. 


( 7 ) 
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b 22. = b 22._ 1 + *xb b 32 i l ' *zb b 12 i _ 1 

b 23j = b 23 i _ 1 + *xb h 33 i _ 1 " *zb *>13^ 

*>3^ = *>31^ + € yb * € xb b 21i_i 

b 32i = b 32 i _ 1 + € yb b 12 i _ 1 " € xb b 22 i _ 1 

b 33i = *>33^ + *yb b 13i_i ' € xb 


where « xb , «yb, * z b are the attitude error outputs of the guidance equations in 
capsule coordinates. 

This detailed derivation need not be presented since the assumption was 
made previously that the rotations were commutative. However, equation (7) 
will be used in subsequent discussions. Note that the relationship between 
spacecraft and launch vehicle quantities has been previously stated. 

e xb = € x 

€ yb = € z 

€ zb = “ € y 

B. ANGULAR ROTATION (AUTOPILOT TIE-IN COUPLED WITH 
SOLVING OF MOMENT EQUATIONS) 

In the general simulation where the inclusion of an autopilot model is 
desirable, additional computations are made. This simulation uses com¬ 
bined analog and digital computers. In this discussion, we do not intend to 
define the interface or organization of computations between analog or digital 
equipment but to define the equations which must be solved. 

When we solve the moment equations, we must determine some 
physical properties of the vehicle specifically, moments of inertia, location 
of applied forces, and locations of the center of gravity. Recall from past 
discussion that functions of moments of inertia I x , Iy, I z , location of center 
of pressure X cp , and location of the center of gravity are stored in the com¬ 
puter. To determine moment arms, the distance from vehicle longitudinal 
axis to engine gimbal (Zj) and engine gimbal station (X s ) are given. 
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Thus, the moment arms of interest are 


Zl 

Z 2 

X 1 = x cg ' X s 
X 2 = X cg “ X cp 

The sign is changed because stations are measured from the vehicle 
nose aft while the x-axis is towards the vehicle nose. 


From the generalized diagram below, we can derive the moments 
generated by applied forces. 



Y f x 


/'y 
x v' 


Consider a positive 
moment which would rotate 
the vehicle in a positive 
direction about the axis in 
question according to the 
right-hand screw rule: 


M x = Zf z y - Zf y z 

My = 2f x Z - Ef z X 

M z = DfyX - 2f x y 

From previous discussions, we have expressions for the three force 
components due to engines A and B and aerodynamic forces. 
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Thus, 

M x = - T Ay Z 1 - T By Z 2 

My = T Ax Z, + t Bx Z 2 - T Az X 1 - T B z X 1 - F Az x 2 
M z = T^yXj + TgyX^ + F ay x 2 

From Eulers'dynamical equations, 

M x = Ix^x - w y « z (I z - Iy) 

My = Iy W y - "x ^X " ^ 

M z = h *z ~ w x w y (V - J x) 
or 



Vehicle rate, cj i? is generated by integrating oq. The results of this 
integration provide angular rotation rates about the body coordinate system 
axes. From the previous derivation [equation (3-7)] we can find the change 
in body-platform orientation defined by matrix pyj : 

0 

[*>!)] = - u *b 

“yb 


w zb " w yb 

0 "xb [ b ij] 

-“Xb 0 
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and 


N 


0 "zb ~"yb 
""zb 0 "xb [\] 

"yb -"xb 0 


To update matrix Tb^l we take the first three terms of the Taylor 
series expansion: L J 

At 2 

b ij t = bi 3 t -1 + h^A 1 + b ij ~ 

The differential equations for angular rotation are solved in vehicle 
coordinates, while matrix Qjjjl relates capsule to platform orientation. 
This requires that vehicle angular rotation be transformed to capsule 
rotations: 


w 


xb = "x 


w yb = "z 
"zb = '"y 

XIII DETERMINATION OF PLATFORM GIMBAL ANGLES 

We can assume that yaw angle iff is small -- nominally zero - during 
launch. Hence, values of gimbal angles may be obtained from matrix fb^jl 
and the terms defined in Section IV-A of this appendix. L J 

bj 3 = - sin ^ 

* = sin' 1 (-bi 3 ) -9001 f l 90° 
b^2 - sin 8 cos if/ 

bjj cos 8 cos iff 
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- sin 0 bj 2 

cos 9 bjj 

( 

b 03 sin cos tf/ 

b 23 - cos <j> cos \f/ 

sin <f> , b 33 

-cos <f> b23 

Since <j> and 9 will not be 180° during launch, 9 = 0 if bj 2 = 0 and 
t p= 0 if b 33 = 0 . 
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AUTOPILOT AND RIGID BODY SIMULATION -- ANALOG PROGRAM 
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Appendix B 

AUTOPILOT AND RIGID BODY SIMULATION -- ANALOG PROGRAM 


A. ASCENT GUIDANCE COMBINED SIMULATION 

In the hybrid simulation, the analog computer was used to simulate the 
autopilots, engine servos, and the rigid body model of the GEMINI launch 
vehicles in pitch, roll, and yaw channels. The operations performed are 
shown in the block diagram in Figure B-l, which indicates all the quantities 
used and transferred by each block. The contents of each block in Figure 
B-l are described in this appendix. 

The autopilots, engine servos, rate feedback filters, and associated 
gains have similar characteristics from channel to channel. A typical 
channel (roll, pitch, or yaw) is illustrated in Figure B-2; all gain factors, 
time constants, and transfer functions representing these items are given. 

In order to check the operation of each of the closed loops, each loop 
was independently subjected to a unit step input. For the purposes of the 
test, the rigid body model was assumed to be a constant gain factor. The 
response of the pitch and yaw channels is presented for times before and af¬ 
ter the gain change (at 105 sec. of flight time) in Figures B-3 and B-4. The 
gain of the roll channel is constant, and step response is presented as in 
Figure B-5. The areas on the analog computer diagram (Figure B-6), where 
the autopilot, engine servo, and rate filter transfer functions were simulated, 
are labeled and enclosed with dotted lines. The attitude error inputs A 8 , 

A*/', and A'/' in launch vehicle coordinates, are indicated in Figures B-l and 
B-6. 

The equations representing the rigid body model of the launch vehicle 
have been derived in Appendix A and are as follows: 

1. Moment arms for thrust and aerodynamic components of force, 
respectively: 

Xi = x cg - x s 

X2 = X cg - X C p 
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2. Thrust components along the x, y, and z axes for both engines:* 

T ax = T a cos « i cosS 

T ay = T a sin € 1 

T az = -T a sin 8 

T bx = T b cos € 2 cos 8 

T by = T b sin € 2 

T bz = _T b sin 8 

3. Euler's dynamical equations: 

. M x 
w x - 

*x 

w y = ly 1 [^y “ W z “x ( ! X “ !z)] 

"z^V 1 [Mz-"y"x(Iy-y] 

The effective thrust moment arm, Xi, is a function of time because the 
motion of the center of gravity is dependent upon mass flow rate which was 
assumed constant., Therefore, Xi was programmed on diode function 
generator F61. ; Time was received from the digital program on trunk 41. 
For similar reasons, the moment of inertia about the y-axis was set up on 
function generator F71. The actual functions programmed are illustrated in 
Figure B-7. These functions were derived from data in Appendix G. The 
aerodynamic moment-arm, x 2 , was formed by summing xj, x s , and x^ with 
the proper signs at the input to amplifier 33. x cp was received from the 
digital program on trunk 45. 

The thrust components are somewhat disguised in the analog computer 
diagram, but a check of the equations will provide the reason. First, in the 
equation for the moment about the y-axis, the terms Tax Z 1 and T Bx z 2 ap¬ 
pear. 

But, 


Tax Z 1 + T bx z 2 = T a cos 8 ( cos € 1 ’ cos € 2 >- 


* Bangles are assumed equal for engines A and B 
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Since the engine gimbal angles,and «2> are limited to 4.5deg.(Cos q- 
Cos € 2 will be very small. Typically, neglecting this term altogether would 
cause an error less than 0.01 percent according to actual runs, so the term 
was not included in the computation. By eliminating T^ x and Tg x , we also 
neglect all cosine terms. Since 8, «i, and 2 are less than 4.5 deg. the sine of 
the angle can be assumed equal to the angle in radians which could cause an 
error of no more than 0.1 percent. With these simplifications plus the 
constant thrust level, we can obtain scaled thrust components merely by 
setting a potentiometer following the engine servo output. Potentiometers 
P06, P26, and Q96 were used for this purpose. 

One further simplification was made in the simulation regarding the 
summing and differencing operations performed at the output of the autopilots 
as seen in Figure B-l. Since the engine servos are linear, the summing and 
differencing can be moved to the servo output. If we assign values to the 
quantities at the roll and yaw engine servo outputs and trace the flow through 
to the formation of the moment about the x-axis, the summing and differencing 
is unnecessary since a similar operation is performed at this point which 
cancels the first (zj = -Z 2 ). The net effect is that twice the quantity in the 
roll channel appears to be used in computing the x component of moment, and 
twice the quantity in the yaw channel appears to contribute to the moment 
about the z-axis. 

With these simplifications, it is relatively easy to find the total moments 
about each axis. The only additional quantities needed are the aerodynamic 
forces which are supplied on trunks 46 and 47. The scaled moments about 
the x, y, and z axes appear at the outputs of potentiometer 26 and amplifiers 
21 and M52, respectively, in Figure B-6. 

The moments and moments of inertia are next inserted in Euler's 
dynamical equations to yield angular rates about the launch vehicle axes. 

The moment of inertia about the x-axis is constant, and the moments of 
inertia about the y-and z-axes are assumed equal. Rates about x-, y-,and z- 
axes appear at the outputs of integrators 30, 00, and 01, respectively, and 
then are sent with engine gimbal angles to the digital program on the trunks 
indicated in Figure B-6. 

The gain change at 105 sec. was implemented using comparator relays, 
which have 105 V and time as inputs. The method in which information was 
exchanged between computers is described in the following section. 
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B. USE OF SIMULATION LABORATORY INTERCONNECTION EQUIPMENT 


The interface between the analog computer and the IBM 7090 DPS 
consists of four units: 

(1) D/A and A/D converters (digital-to-analog and analog-to-digital) 

(2) IBM 737 core memory 

(3) Data Controller (DACON), 

(4) Real-Time Control Unit (RTCU) 

Data is transferred between computers by feeding information from 
the 7090 through the RTCU to specified locations in the 737 memory 
(Figure B - 8). This information is taken from memory by DACON through 
the RTCU and sent to the D/A converters, one channel at a time. When the 
transfers to the analog computer are complete, DACON takes information 
from the A/D converters and relocates it in other memory locations, one 
channel at a time. This operation completes one cycle of the DACON pro¬ 
gram. The order in which the data transfers are made can be programmed 
on the DACON control board. 

The memory now contains information and is prepared to transmit to 
the 7090. When the 7090 requests information from the memory, data transfer 
between the 737 memory and the 7090 is initiated by an interrupt signal, which 
comes from the RTCU at fixed time intervals. This signal produces an 
unconditional GO TO statement in the 7090 program which causes data to be 
transferred. When the 7090 has received and transmitted information to the 
737 memory, it uses this information in the computation of new program 
results. The digital program must complete its computation cycle using the 
new information before another interrupt signal is transmitted by the RTCU. 
The analog computer was time-scaled by a factor of 4 to 1 because of the 
length of the digital program cycle (1 sec. analog time equals 4 sec. real 
time). The length of the digital program computation cycle was dictated by 
the large printout requirements. 


B-6 



B-7 


Pitch 



Figure B-l. Analog Portion of Ascent Guidance Combined Simulation 
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Figure B-2. Autopilot Transfer Function 
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Figure B-6. First Stage Ascent Guidance Analog Computer Diagram 



























Figure B-8. Simplified Diagram of Hybrid Simulation 
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Appendix C 

AUTOPILOT AND RIGID BODY SIMULATION--DIGITAL PROGRAM 


I. INTRODUCTION 

This appendix describes the digital simulation of the GEMINI launch 
vehicle autopilot and rigid body dynamics. The program consists of a con¬ 
trol program and eight subroutines for generating the required equations. 

The autopilot model used for simulation is shown in Figure C-l. This 
model also applies to the pitch, yaw and roll channels for both Stage 1 and 
Stage 2. Since the analysis of Stages 1 and 2 are similar, only Stage 1 is 
discussed in detail. 

The autopilot equations are developed in this appendix. The derivation 
of the equations used for rigid body dynamics is included in Appendix A. Cer¬ 
tain effects -- platform gimbal and rate gyro dynamics, and the gain amplifier 
lags in the autopilot — have not been included because they are believed in¬ 
significant and beyond the scope and purpose of the simulation. 

The symbols used in this program are consistent with those in Appendix 
A. However, a glossary has been attached which defines the "A" quantities 
not included elsewhere. 



Figure C-l. GEMINI Launch Vehicle Model 
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II. DETERMINATION OF DIFFERENCE EQUATIONS 


The servo dynamics in the pitch channel of Stage 1 (Figure C-2) have 
been reduced from compatibility data, supplied by Martin Aircraft Co. to 
fourth order by ignoring poles at (S + 220550.0) and (S + 511.67), and a 
x^zero at (S + 9918. 7). For the present time, the rigid body dynamics will be 
( represented by a first order integrator. The root locus plot of the continuous 
system (Figure C-2) is shown in Figure C-3 and is a function of (K r K b ). 1 
To obtain difference equations which accurately approximate the continuous 
system, a sampling rate was chosen in which the root loci of the sampled 
system and the continuous system agreed to within a reasonable tolerance. 
The z transforms of Figure C-4 are for a sampling period of 0.01 sec. Note 
that a one-period time lag has been introduced in the feedback loop. This 
lag is necessary to offset the inherent phase lead of the z transformations. 

If you compare the root loci of the sampled system (Figure C-5) and the 
continuous system (Figure C-3), you will observe that the two systems are 
essentially identical in phase and gain over the specified frequency range. 
Since the operating point (Kp Kg) is not expected to exceed about 5.0, the 
control poles of the z-transformed system will be accurate to better than 1 
percent. This accuracy is believed to be considerably higher than the data 
on which the equations are based. 

The next step is to program the difference equations resulting from the 
z-transformed system. A convenient and desirable programming check 
requires only the comparison of the response for a specific gain to that pre¬ 
dicted by the root loci. A good choice of gain from Figure C-5 is 21.05, 
which causes a conditionally stable system. The programmed system ’ 
response to a step input for this value of gain is shown in Figure C-6 and 
oscillates with no damping at 13. 3 radians/sec. as was indicated by the 
root locus plot. 



Figure C-2. Pitch Channel (Continuous) 


* All root loci and z transforms were obtained from an IBM 7090 program. 
See ’’Numerical Methods for the Synthesis of Linear Control Systems" by 
M.E. Fowler, IBM No. 62-907-417. 
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1 


The remaining pitch and yaw channels are treated in a similar manner. 
For the vehicle in question, the servo dynamics are identical for all three 
channels, and the filter constant (t ) is the same for the pitch and yaw 
channels. Thus, the yaw channel difference equations are the same as those 
for pitch. For the roll channel, the filter constant is 122.0 rather 17.0; 
therefore; different coefficients for the z transforms must be obtained. With 
the same sampling period of 0.01 sec. , the sampled system for roll is shown 
in Figure C-7. 

Now that the difference equations for the three channels have been 
obtained, programmed, and checked; the rigid body approximations are 
replaced by the true rigid body dynamics which include coupling. Figure 
C-8 shows the response of the final program to a step input at 100 sec. 



Figure C-4. Pitch Channel (Sampled at T = 0.01) 
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Figure C-6. Transient Response of Sampled System 
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Figure C-7. Sampled Roll Channel (T = 0.01) 
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Figure C-8. Transient Response in Final Program 
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III. COMPUTER PROGRAM 


The nine subroutines, which will be referred to collectively as the 
vehicle program, are listed below, and the primary functions of each are 
discussed. Because the vehicle program difference equations are for a 
sampling period of 0.01 sec. and the main program sampling period is 0.05 
sec., the vehicle program completes five cycles for every cycle of the main 
program. 

A. SUBROUTINE MC 

This subroutine provides communication between the main program and 
the remaining subroutines of the vehicle program. Since the main program 
and the vehicle program were not written by the same individual, this appeared 
to be the most convenient means of joining the two programs. The main pro¬ 
gram calls subroutine MC every 0.05 sec. which, in turn, calls the re¬ 
maining subroutines of the vehicle program. Received from the main pro¬ 
gram are the three error commands for pitch, yaw, and roll — the location 
of the center of pressure and the launch vehicle y and z components of the 
aerodynamic force. The three rates and the three engine gimbal angles are 
returned to the main program 0.05 sec. later in flight time. These quantities 
are the averaged values over the 0.05-sec. sampling period. This sub¬ 
routine also switches from Stage 1 to Stage 2 vehicle dynamics through 
sensing of main environment conditions. 

B. SUBROUTINE CONI 

Entry to this subroutine is made only once — the first time the vehicle 
program is called. It's function is to define the values of constants used in 
remaining subroutines for the Stage 1 equations. 

C. SUBROUTINE SD1 

SD1 computes the engine gimbal angles by solving the forward loop 
difference equations supplied with the error command signals. When nec¬ 
essary, it limits them to the maximum specified value of 4. 5 deg. The 
computation time lag (T 2 ), shown in Figure C-l, is also handled in this 
subroutine. According to whether the read-in value of the constant A(10) 
is 0, 1, or 2, the time delay will be 0, 10, or 20 msec, respectively. The 
final duty of SD1 is to change the values of autopilot gain constants Kr and 
Kj) at the end of 105 sec. of flight time. 
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D. SUBROUTINE RBI 


The rigid body moment equations are computed in RBI. The moment 
of inertia and the control moment arm (Distance from C. G. to point of thrust 
application) are both obtained from polynomials which were fitted by the 
method of least squares with weight as the independent variable. For the 
Stage 1 both polynomials are fifth degree; they are fourth degree for Stage 2. 
Plots of these polynomials as a function of time are shown in Figures C-9 
and C-lOfor Stages 1 and 2, respectively. These plots are obtained directly 
from data presented in Appendix G. 

The moment equations for Stage 1 are as follows: 

Mx =■ -TZ C ( sin « i - sin « 2 ) 

My = T [-2X C sin 8 +Z^. (cos « 1 - cos < 2 )] - F az x a 

Mz = T X c (sin « 1 + sin « 2 ) + F ay Xa 

Included in these equations is the approximation that (cos 2 « - sin 2 S )s 
(cos « cos 8) for the small angles involved. Fay and Faz are the aerodynamic 
forces and Xa is the aerodynamic moment arm; all of which are supplied by 
the main program. 

The moment equations for Stage 2 are as follows: 

Mx = Ti Y c sin r) 

My = -T X c sin 8 + Tj X c sin 77 
Mz = T X c sine 

Tj is the roll nozzle thrust which is located on the y-axis, Y c feet from the 
vehicle center line. The aerodynamic forces are considered to be zero for 
Stage 2. 

The rates are then obtained for both stages from the equations: 
u x = Mx/Ix 

«i» y = (My + UJ Z W X (Iz - lx) ) / Iy 
w z = (Mz + ui x w y (lx - Iy) ) / Iz 

E. SUBROUTINE RL1 

RLl simply solves the rate loop difference equations. 

F. SUBROUTINES CON2, SD2, RB2, RL2 

These subroutines perform in a similar manner to their Stage 1 
counterparts when the main environment indicates initiation of Stage 2. 
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Figure C-9. Stage 1 Vehicle Characteristics 
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Figure C-10. Stage 2 Vehicle Characteristics 







GLOSSARY OF "A" DEFINITIONS 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 


11 . 


12 . 

13. 

14. 

15. 

16. 


21 . 

22 . 

23. 

24. 

25. 

26. 

27. 

28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 

41. 


Time 

T (0.01 sec) 


K D (A0) 

k d (^) 

k d N>) 

Delay Control (0 = No Delay, 1 = 0.01 sec Delay, 2 = 0.02 sec 
Delay) 

fz c (Distance from thruster t to vehicle t; Stage 1) 
t Y c (Distance from vehicle t to Ti thruster; Stage 2) 

T (Engine Thrust) 


WT. 

X c (Control moment arm - Distance from STA 1274 to C. G. (Ft.) ) 
Tj (Roll nozzle thrust - Stage 2) 

Kr (Roll Channel) 

Kr (Pitch and Yaw Channels) 

Kd (Roll Channel) 

Kd (Pitch and Yaw Channels) 

Z0 num. coeff. (Pitch and Yaw forward loop) 

Z 1 num. coeff. (Pitch and Yaw forward loop) 

Z 2 num. coeff. (Pitch and Yaw forward loop) 

Z 3 num. coeff. (Pitch and Yaw forward loop) 

Z4 num. coeff. (Pitch and Yaw forward loop) 


A0 (Pitch error command) 
Av^ (Yaw error command) 

A<f> (Roll error command) 

ZO Denom. coeff. (Pitch and 
Zl Denom. coeff. (Pitch and 
Z 2 Denom. coeff. (Pitch and 
Z 3 Denom. coeff. (Pitch and 
Z 4 Denom. coeff. (Pitch and 


Yaw forward loop) 
Yaw forward loop) 
Yaw forward loop) 
Yaw forward loop) 
Yaw forward loop) 


Iyy and Izz 
Ixx 


8 F (Pitch feedback loop value) 
8 (T-5) (Pitch channel) 
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Glossary of "A" Definitions (cont) 


42. 8 (T-4) (Pitch channel) 

43. 8 (T-3) (Pitch channel) 

44. 8 (T-2) (Pitch channel) 

45. 8 (T-l) (Pitch channel) 

46. 8(T) (Pitch channel) 

47. «1 (Yaw angle of engine No. 1 Stage 1) 

48. 6 2 (Yaw angle of engine No. 2 Stage 2) 

49. w y (Pitph channel) 

50. u y or 9 or q 

51. f * (T-5) (Pitch channel) 

52. e ’ (T-4) (Pitch channel) 

53. € ' (T-3) (Pitch channel) 

54. € ' (T-2) (Pitch channel) 

55. € ' (T-l) (Pitch channel) 

56. « ’ (T) (Pitch channel) 

57. - 


59. 

60. 
61. 
62. 

63. 

64. 

65. 

66 . 

67. 

68 . 


74. 

75. 

76. 

77. 

78. 

79. 

80. 
81. 
82. 

83. 

84. 


V f (Yaw feedback loop value) 
8 (T-5) (Yaw channel) 

8 (T-4) (Yaw channel) 

8 (T-3) (Yaw channel) 

8 (T-2) (Yaw channel) 

8 (T-l) (Yaw channel) 

S (T) (Yaw channel) 


w z (Yaw channel) 
w z or '/'or r 
« ’ (T-5) (Yaw channel) 

« ’ (T-4) (Yaw channel) 

* ’ (T-3) (Yaw channel) 

< ’ (T-2) (Yaw channel) 
e ’ (T-l) (Yaw channel) 

6 ’ (T) (Yaw channel) 

Xa (Aerodynamic Moment arm (Ft.) ) 

Fay (Aerodynamic force in Y direction) 

Faz (Aerodynamic force in Z direction) 
Location of C.P. (Ft.) from Aft End 
Z° num. coeff. (Roll channel forward loop) 
Z 1 num. coeff. (Roll channel forward loop) 
Z2 num. coeff. (Roll channel forward loop) 
Z3 num. coeff. (Roll channel forward loop) 
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Glossary of "A” Definitions (cont) 


85. Z 4 num. coeff. (Roll channel forward loop - Stage 1) 

86 . - 

87. Mx (Moment about X axis) 

88. My (Moment about Y axis) 

89. Mz (Moment about Z axis) 

90. - 

91. Z° denom. coeff. (Roll channel forward loop) 

92. Zl denom. coeff. (Roll channel forward loop) 

93. Z 2 denom. coeff. (Roll channel forward loop) 

94. Z 3 denom. coeff. (Roll channel forward loop) 

95. Z 4 denom. coeff. (Roll channel forward loop - Stage 1) 


99. 

100 . 

101 . 

102 . 

103. 

104. 

105. 

106. 

107. 

108. 

109. 

110 . 
111 . 
112 . 

113. 

114. 

115. 

116. 

117. 

118. 

119. 

120 . 
121 . 
122 . 

123. 

124. 

125. 


4>y (Roll feedback loop value) 

8 (T-5) (Roll channel) 

8 (T-4) (Roll channel) 

8 (T-3) (Roil channel) \ Stage 1 
8 (T-2) (Roll channel) { 

8 (T-l) (Roll channel) 1 
8 (T) (Roll channel) J 


(T-l) 

(T) J 


Stage 2 


<1) x (Rol,l channel) 
u» x or 4> or P 
« ' (T-5) (Roll channel) 
e ' (T-4) (Roll channel) 

6 ' (T-3) (Roll channel) \ Stage 1 
e ' (T-2) (Roll channel) f 
« ’ (T-l) (Roll channel) 

€ ' (T) (Roll channel) J 


€ ' (T-4) 

* ’ (T-3) 

« ' (T-2) Stage 2 

* • (T-l) 

€ ’ (T) 


Sum of w y (A(50) ) over one sample period (0.05 sec) 

Sum ofwz (A(70) ) over one sample period (0.05 sec) 

Sum ofu>x (A(H0) ) over one sample period (0.05 sec) 

Sum of 8 (A(46) ) over one sample period (0.05 sec) 

Sum of e i (A(47) ) over one sample period (0.05 sec) - Stage 1 

« (A(66) ) - Stage 2 # v „ 

Sum of e 2 (M 48 ) ) over one sam P le period (0.05 sec) - Stage 1 

A(105) ) - Stage 2 
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PARK BEDFORD 


* LABEL 
CMC 

SUBROUTINE MC (NOT USED) 

DIMENSION OV1(24),DV2(20),DV3(131),FA(3),SF1(3),SF(2),A(125) 

C »DV5(100) 

COMMON DV1»TI ME,DV2,THETD,YAWD,ROLLD,DV3,EXX,EYY ,EZZ,EPSIL1»EPSIL2 
C ,DELTA2,XCP,FA,SF1,SF,SFT,DV5,A 

IF (TIME-.05) 1,1,2 

1 CALL CONI 

2 A(28) =-EYY 
A(29) =—EZZ 
A{30) =—EXX 

A(78) * FA(2) 

A(79) = FA(3) 

A{80) « ( 1274.-XCPJ/12. 

CALL SD1 

YAWD = AM 18J/5. 

ROLLD = A(119J/5- 
THETO = A(1201/5. 

DELT A2= -A(121)/5. 

EPS IL1= -A(122)/5. 

EPS IL2= -A(1231/5. 

RETURN 

END 

• LABEL 
CSD1 

SUBROUTINE SOI 

DIMENSION DV1(24),DV2(20),DV3(131),FA(3),SF1(3),SF(2),AC125) 

C ,DV5(100) 

COMMON DVl,TIME,0V2,THETD,YAWD.ROLLD,DV3,EXX,EYY,EZZ,EPSIL1,EPSIL2 
C ,DELTA2,XCP,FA,SF1,SF,SFT,DV5,A 

DO 20 1=1,6 
20 A(1+117) = 0. 

N = A(10) +1. 

DO 13 K= 1,5 
A(1) = A(1)+A{2) 

DO 1 1=1,5 
A(1+40) = A(1+41) 

A(1+50) = A(I+51) 

A(I+60) = A(I+61 ) 

A(1+70) = AlI+71) 


C-18 



1 


A(1 + 100) •= A( 1 + 101) 

A(1 + 1 10) = A(1 + 111) 

GO TO 12,3,4,5),N 

2 A(7) = A(20)*A(28) 

A(8) = A(20)*A(29) 

A(9) = A(19)*A(30) 

N = 4 

GO TO 5 

3 N = 1 
GO TO 5 

4 N = N-1 

5 A(56) = A(7)—A(40) 

A(76) = A(8)-A(60) 

A(116) = A(9)-A(100) 

A(46) = 0. 

A(66) = 0. 

At 106) = 0. 

DO 6 1=1,5 
J = 6-1 

A(46) = A146)+A tJ+20)*A(J + 50)-AtJ + 30)*A(J+40) 

A(66) = A(66)+A(J+20)*A(J+70)—A(J+30)»A(J+60) 

6 A(106) = A{106)+A(J+80)*A(J+110)—A CJ+90)*A(J+100) 

A147) = A(66)+A{106) 

A{48) = A(66)-A(106) 

IF (A8SF(A(46))-.078539816) 8,8,7 

7 A14 6) = SIGNFl.078539816,A(46)) 

8 IF (ABSF(A147))-.078539816) 10,10,9 

9 A ( 47) = SIGNFl. 078539816, A(47)) 

10 IF (ABSF(A(48))-.078539816) 12,12,11 

11 A(48) = SIGNFl.078539816,A(48)) 

12 CALL RBI 

13 CALL RL1 

IF (All)—105.) 15,14,14 

14 Al20) = .21 

A 118) = .165 

15 RETURN 
END 

* LABEL 
CRB 1 

SUBROUTINE RBI 

DIMENSION DV1(24),DV2(20),DV3(131),FA(3),SF1{3),SF(2),A(125) 
C ,DV5(100) 


C-19 



COMMON DVl,TIME?DV2tTHETDtYAWD»ROLLD»DV3,EXXfEYY»EZZtEPSILl,EPSIL2 
C ,DELTA2,XCP,FA,SFl,SF,SFT,DV5,A 

A( 14) = 329000.-1662. 1622*A(1) 

WT2 = A( 14)*A(14) 

WT3 = WT2*A(14) 

WT4 = WT3*A(14) 

WT5 = WT4*A(14) „„„„ 

A(15) = 154.86-.0017201*A(14)+.10734E-7*WT2-.3329E-13*WT3+. 493255 
C E-19*WT4-.26337E-25*WT5 

A{37) * -999960.+11.074*A(14)+.0003861*WT2-.23705E-8*WT3+.50698E- 
C 14*wT4-.36014c-20*WT5 

A{77) « A(80)-A(15) 

A {87) = A{12)*A(11)*(SINF(A(47))—SINF(A(48))) 

A(88) - At12)*(2.*A{15)*SINF(A(46))+(COSF(A(47))-C0SF(A(48)))*A{11 
C ))-A{79)*A(77) 

A(89) = A(12)*A(15)*C SINF CAC47))+SINF(A(48)))+A(78)*A(77) 

AC 109)= A(87)/A(38) 

A £ 49) = {A(88)+A(70)*A(110)*(A{37)-A(38)))/A(37) 

A(69) = (A(89)+A(110)*A(50)MA{38)-A(37)))/A(37) 

A(50) = A(2)*A(49)+A(50) 

A(70) = A(2)*A(69)+A(70) 

A(1 10)= A(2)*A(109)+A(110) 

A( 1 18) = A(118)+A(50) 

A( 119) = A(119)+A(70) 

A(120) = A(120)+A(110) 

A(121) = A(121)+A(46) 

A{122) = A(122)+A(47) 

A(123) = A(123)+A{48) 

RETURN 

END 

LABEL 


CRL1 

SUBROUTINE RL1 

DIMENSION DV1(24),DV2(20),DV3(131),FA(3),SF1(3),SF(2),A(125) 

C .DV5(100) 

COMMON DVl,TIME,DV2,THETDtYAWD,R0LLD,DV3,EXX,EYY.EZZ,EPSIL1,EPSIL2 
C ,DELTA2,XCP,FA,SF1,SF,SFT,DV5,A 

A(40) = .15633518*A(50)*A(18)+.84366482*A(40) 

A{60) = .15633518*A(70)*A(18)+.84366482*A(60) 

A( 100)= .7047 6983*A(110)•A(17)+.2952301*A(100) 


RETURN 


END 
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PARK BEOFORD 


LABEL 


CKC 


1 

2 


SUUKUUlINt: 

DIMENSION DVl(24),DV2(20),DV3(I 31) f FA(3),SF1(3)»SF(2),A<12b) 
C ,DV5(100) 

COMMON DV1,TIME,DV2,THETD,YAWD,ROLLD,DV3,EXX,EYY,tZZ,EPSILl, 
C ,DELTA2,XCP,FA,SF1,SF,SFT.DV5.A 

IF(LCI) 1,1,2 
CALL C0N2 
LC1 = 2 
A (28) =-EYY 
A (29) =—EZZ 
A(30) =—EXX 
CALL SD2 

YAWD = A(118)/5. 

ROLLD = A(119)/5. 

THETD = A(1201/5. 

DELTA2= At 1211/5. 

EPSIL1= At 1221/5. 

EPSIL2= At 1231/5. 

RETURN 

END 

LABEL 


CC0N2 

DIMENSION^Vl1241,DV2 1201,DV3(1311,FA(31,SF1(31,SF(21,A(1251 
C ,DV5(100) 

COMMON DV1,TI ME,DV2,THETD,YAWD, ROLLD ,DV3,EXX,EYY,EZZ , EPS I LI , 
C ,DELTA2,XCP,FA,SF1,SF,SFT,DV5,A 


DO 1 1=1,9 

1 A(I 1 = 0. 

DO 2 1=11,125 

2 All) = 0. 


A ( 2 1 = 

.01 

Alii) = 

3.67 

A C12) = 

100000. 

A t 1 61 = 

565. 

At 171 = 

5. 

At 181 = 

.144 

At 191 = 

10. 

A120)= 

.216 


EPSIL2 


EPSIL2 
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A(21)=- 

.16298145E-8 

AC 22) = 

.00022354559 

A(23)= 

.0048435752 

A(24)= 

.010580247 

A (25) = 

.002042806 

A(31)=- 

.026487565 

A(32)* 

.11564639 

A(33)=- 

.91904442 

A(34)* 

2.4805997 

A(35)=- 

2.6330241 

A(38) * 

5159. 

A < 81) = 

.13299286E-5 

A182) = 

.055928213 

A(83)* 

.10959722 

A(84)= 

.20813586 

A (9 1) = 

.066116963 

A{92)=- 

.12566 

A(93)=- 

.35719627 

A(94)=- 

.21303089 

A(50) = 

-THETD 

A(70) - 

YAWD 

A(110) 

* ROLLO 

A(118) 

=5.*A(50) 

A C119) 

=5.*A{70) 

A(120) 

=5.*A(110) 

RETURN 


END 


LABEL 


SUBROUTINE SD2 

DIMENSION DV1(24),DV2(20), 


tOV5(100) 

COMMON 

DV1,TIME,DV2,THETD, 


,DELTA2,XCP,FA,SF1, 

DO 20 I 

-1,6 

A(1+117) = 0. 

N » A(10)+1. 

00 15 K 

-1,5 

A(l) = 

A(1)+A(2) 

DO 1 1 = 

1,5 

A( 1+40) 

= AC 1+411 

A CI+50) 

* Ad +51) 


1DV3(131),FA(3) 

YAWD,R0LLD,0V3, 
SF» SFT»DV5*A 


SF1 

EXX 


3),SF(2),A(125) 
EYY,EZZ,EPSIL1,EPSIL2 


C-22 



A{1+60) - A(I+61) 

1 AU+70) = A(I+71) 

DO 2 1=1,4 

AU + 100) = A ( 1 + 101) 

2 Al1+110) = A(l+111l 
GO TO (3,4,5,6),N 

3 A(7) » A(20)*A(28) 

A(8) ■ A(20)*A(29) 

A(9) - A(19)*A(30) 

N * 4 

GO TO 6 

4 N = 1 
GO TO 6 

5 N = N-1 

6 A(56) = A{7)—A(40) 

A(76) = A(8)-A(60) 

A{115) = A(9) - A(100) 

A(4 6) = 0. 

A(66) = 0. 

A(105) = 0. 

DO 7 1=1,5 
J = 6-1 

A(46) = A(46)+A(J+20)*A(J+50)-A(J+30)*A{J+40) 

7 A(66) = A(66)+A(J+20)*A{J+70)-A(J+30)*A(J+60) 

DO 8 1=1,4 

J = 5-1 

8 AC 105) = A(105)+ACJ + 80)*A(J+l10) — A CJ+90)*A{J+100) 

IF (ABSF(A(46))-.038920842) 10,10,9 

9 A (46) = SIGNFU 038920842, A(46) ) 

10 IF (ABSF(A(66))-.038920842) 12,12,11 

11 A(66) = SIGNFU038920842,A{66) ) 

12 IF (ABSF(A(105))-.57595865) 14,14,13 

13 A(105 ) = SIGNFl.57595865,A(105)) 

14 CALL RB2 

15 CALL RL2 
RETURN 
END 

* LABEL 

CRB2 

DIMENSION DVl(24),DV2(20)»DV3(131),FA(3),SF1(3),SF(2),A(125) 
C ,DV5(100) 
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2 

3 


COMMON DV1,TIME,0V2,THETD,YAWD,R0LLD,DV3,EXX,EYY,EZZ,EPSIL1,EPSIL2 
C ,DELTA2,XCP,FA,SFl,SF,SFT,DV5,A 

At 14) » 71500.-344.118*A(1) 

WT2 * A( 14)*A(14) 

WT3 * WT2*A(14) 

WT4 = WT3*A( 14) 

A(15) * 29.555-.0011695*A(14)31521E-7*WT2-.37686E-12*WT3+.16829 
C E-17*WT4 

IF lA(14)—35000.) 2,1,1 
A(37) * 103000. + .9041 1 *CAC 14)-35000.) 

M37? = -86553.+21.493*A(14)-.0010424*WT2+.23863E-7*WT3-.20586E-12 


A(87) * 
A(88) * 
A(89) - 
A(109)= 
A(49) » 
A(69) = 
A(50) « 
A(70) » 
A (110 ) = 
A (118) < 
A (119) > 
A(120 ) - 
A(121) 

A ( 122) 

A(123) 
RETURN 
END 
LABEL 


*WT4 

A(16)*A(11)*SINF(Al105)) 

A(15)•(A(12)*SINF(A(46))+A(16)*SINF(A(105) 
A(12)*A(15 )*SINF(A(66)) 

A(87)/A(38) 

(A(88)+A(70)*A(110)*(A(37)-A(38)))/A(37) 

(A(89)+A(110)*A(50)*(A(38)-A(37)))/A(37) 

A(2)*A(49)+A(50) 

A(2)*A(69)+A(70) 

A(2)«A(109)+A(110) 

* A(118)+A(50) 

. A(119)+A(70) 

* A(120)+A(110) 

* A(121)+A(46) 

« A(1221+A(47) 

- A(123)+A(48) 


)) 


CRL2 

DIMENSION DV1(24),DV2(20),DV3(131),FA(3),SF1(3),SF(2),A(125) 

^COMMON DV1|tiME,DV2,THETD,YAWD,R0LLD,0V3,EXX,EYY,EZZ,EPSIL1,EPSIL2 
c ,DELTA2,XCP,FA,SF1,SF,SFT,DV5,A 

A140) = . 15633518*A(50)*A{ 18)+.84366482*A(40) 

A(60) = .15633518*A(70)*A(18)♦.84366482*A(60) 

A(100)= .70476983*A(110)*A(17)+.2952301*A(100) 

RETURN 

END 
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C JAN. ,1962 

CMC 

SUBROUTINE MC 

DIMENSION DV1(24),DV2(20),OV3(106),FA(3),SF1(3)» SF(2),A(125) 
C ,DV5(100),DV6(24) 

COMMON DV1 ,TIME,DV2,THETD, YAW0,R0LLD,D\/3, JUP,DV6,EXX,EYY,EZZ 
C,EPSIL1,EPSIL2,DELTA2,XCP,FA,SF1,SF,SFT,DV5,A 
GO TO(1,4,4),JUP 

1 IF(LCl) 2,2,3 

2 CALL CONI 
LC1 = 2. 

GO TO 3 

4 IF(LC2) 5,5,3 

5 CALL C0N2 
LC2 * 2. 

LC3 - 2. 

3 A(28) = -EYY 
A(29) * -EZZ 
A(30) = -EXX 
A(78) * FA(2) 

A(79) * FA(3) 

A(80) ■ (1274.-XCPJ/12. 

IF(LC3) 6,6,7 

6 CALL SD1 
GO TO 8 

7 CALL S02 

8 YAWD » A(118)/5. 

ROLLD = A(119)/5. 

THETD = A(120)/5. 

06LTA2 *—A(121)/5« 

EPSIL1 =-A(122)/5. 

EPSIL2 =-A(123)/5. 

RETURN 

END 

* LABEL 

CC0N1 

SUBROUTINE CONI 

DIMENSION DV1(24),DV2(20),DV3(131),FA(3),SF1(3),SF(2),A(125) 
C ,DV5(100) 
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COMMON DV1,TIME,DV2,THETD,YAWD,ROLLD,DV3,EXX,EYY,EZZ,GPS 111,EPSIL2 
C ,DELTA2,XCP,FA,SF1,SF,SFT,DV5,A 

DO 1 1=1,9 

1 A(I) = 0. 

DO 2 1=11, 125 

2 A(I) = 0. 

A { 2) = .01 

A( 10)= 2. 

At 1 1) = 2.4225 
At 12) = 215000. 

At 17)= .07 
At 18)= .4 
At 19)= .31 
A(20)= .65 
A121)=-.23283064E-9 
A t 22)= .00088665786 
A{23)= .016191196 
A t 24)= .027056967 
A{25)= .0042071799 
A(31)=-.074140752 
A{32)= .60805228 
A(33)=-2.0132761 
A(34)= 3.1890013 
At35)=-2.6612953 
A t 38) = 9730. 

A(8 1)= . 18626451E—8 
A(82)= .002755933 
A t 83)= .062395737 
A(84)= .12029903 
A(85)= .024504986 
A(91)=-.025944648 
A t 92)= .26990752 
A(93)=-l.105326 
A t 94)= 2. 1921498 
A(95)=-2.1128606 
RETURN 
END 
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Appendix D 


GEMINI ASCENT GUIDANCE SIMULATION - 
COMPUTER PROGRAM 


I. INTRODUCTION 


This appendix presents the math flow and the FORTRAN listings used 
for the simulation, which has been previously described and derived in Sec¬ 
tion II and Appendix A of this report. In addition, the specific coordinate 
systems and symbols are defined in this appendix. Figure D-l illustrates 
the complete program flow and will be referenced frequently in this discus¬ 
sion. Also, the various printout formats, used at various points in the pro¬ 
gram, are included as Tables D-I through D-IV, symbols as Table D-V and 
FORTRAN listings for the main environment and subroutines as Table D-VI. 

H. VARIATIONS BETWEEN THE PROGRAM AND APPENDIX A 


Several variations exist between the actual program and the description 
in Appendix A: 

• Coordinate Systems — Since the actual program was an 
adaptation of a more general ascent guidance simulation 
program, the coordinate systems had to be changed. 

These coordinate systems are described in Section IV and 
Figure D-2 of this appendix. 

• Symbols — As a result of other changes, the symbols 
vary significantly. Symbols used in the computer program 
are listed in Table D-V. 

• Matrices — Transformation matrices between coordinate 
systems are identified in Section IV. 

• Position and Velocity Update - Basically, Simpson's in¬ 
tegration is employed, but present gravity components 
are estimated, thus providing for more accurate computa¬ 
tion than described in Appendix A. 


D-3 
































North Polar Axis 



Greenwich Meridian 
at Time of Platform Release 


Figure D-2. Relationships of I, J and X Coordinate Frames 


• Vehicle Attitude — The program uses a method of updating 
the vehicle platform matrix which is more accurate than 
the method described in Appendix A. Runge-Kutta inte¬ 
gration with double precision is employed for this purpose. 

• Hold Down — Provisions are incorporated into the math 
flow to permit updating of the position and velocity compo¬ 
nents during hold-down conditions. This provision is an 
addition to the program described in Appendix A. 

IE. PROGRAM SUBROUTINES 


Basically, the simulation consists of the environmental or main pro¬ 
gram, which controls the program and calls the required subroutine as shown 
in Figure D-l. These subroutines, listed below, are accompanied by a 
description of their respective functions. 

1. TLU (Table Look-Up).accepts inputs of Mach number 

and angles-of -attack (on an( f ay) an< * P rov ides outputs of 
aerodynamic coefficients (Ca» Cn an d Cy) and location of 
the vehicle center of pressure (X^p). 

2. OBLATE.accepts inputs of position components in an 

earth-reference coordinate system, RX(1), RX(2), 

RX(3), and provides outputs of gravitational acceleration 
components in the same coordinate system GXE(l), 

GXE(2), GXE(3) . 

3. ATMOD.accepts an input of altitude above sea level 

^z) and provides outputs of air density (DE) and the speed 
of sound (SP). Pressure (PE) is also provided but is not 
currently being used in thrust computations. 

4. RUN-KUT.accepts inputs of body rates (Thet^ ^ 

Roll D, PsiD) and vehicle-guidance frame matrix, |D4(i,j)j 

and provides output of updated matrix,p4(i,j)l . L J 

5 COM.provides interface between the main program 

andthe autopilot subroutines.. Two programs were used; 
one for the digital autopilot simulation, and a second for 
the analog autopilot simulation. 
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Other subroutines used are the FIXER-FLOATER, DATA, PLOT, STOP 
and EXIT. 

IV. DEFINITION OF COORDINATE SYSTEMS (Figure D-2) 

A. lORJ FRAME 

An intermediate north-east oriented inertial reference frame estab¬ 
lished at time of platform release: 

1(1) = east 1(2) = north 1(3) = vertical (+up) 

B. 2 OR I FRAME 

The primary guidance frame, inertial reference defined at time of 
platform release. Basically, this frame is related to the J frame by a single 
rotation about 1(3). Relationship of this frame to the platform frame is given 
below: 

2(1) = 1 X 2(2) = l z 2(3) = -1 Y 

C. 3 OR F FRAME 

An earth-fixed geometric north-east frame with one axis through the 
launch point. Prior to platform release this frame is identical to the J frame. 

D. M FRAME 
Spacecraft-oriented frame: 

M(l) = l x b M(2) = l yb M(3) = l zb 

E. X FRAME 

An inertial frame described by two rotations from the I frame such that 
X(3) is through the polar axis. This frame is used for gravity computations 
and is established at platform release. X(l) and X(2) are in the equatorial 
plane such that X(l) has the same longitude as 2(3) prior to platform release. 

Launch vehicle coordinates are not directly defined in this simulation 
although implicit in several equations. The relationship between launch ve¬ 
hicle and spacecraft coordinates is described in Appendix A. 
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V. 


TRANSFORMATION MATRICES 


A. D3, J-to-I transformation matrix — Defined by angle AZ: 

T = |p3(i, j)J j 

B. D4, Vehicle-to-I transformation matrix — Defined by platform gimbal 
angles 0 b, ^ b, <f> b: 

I = [D4 (l,j)] IS 

C. AP, X-to-I transformation matrix — Defined by anglesX and AZ: 

T = [AP(i,j)]x 

D. AJF, F-to-J transformation matrix — Defined by X and angle cor¬ 
responding to angular rotation of the earth since platform release: 

J = [AJF(i,j)] "F 
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Table D-I 


INPUT DATA FORMAT 


■ 

Position Number j 








1 

2 

3 

4 

5 

1 

SLAT 

SLONG 




2 

AZ 

DEL = 0 

_ 

XO 

ZO 

3 

SF1(1) 

SF 1 (2) 

SF1 (3) 

w A 

C(l) 

SF (1) 

4 

C (2) 

— 

— 

Debug 

5 

SF (2) 

CON 4 

H (1) 

H (2) 

— 

6 

— 

— 

— 

— 

— 

7 

— 

— 

— 

— 

— 

8 

— 

— 

— 

_ 

_ 

9 

— 

E 

C3 

C4 

D5 

10 

SWT 

GC 

EZ 

EXP P 

QPT 

11 

TIM I 

Y0 

— 

— 

— 

12 

— 

— 

— 

_ 

_ 

13 

— 

Y1{1) 

Yl (2) 

— 

— 

14 

— 

— 

— 

— 

X INCL 

15 

16 

«I 

cup 

_ 

— 


17 

18 

19 

SH 

Stage 1 

— 

— 

— 

20 

DFL 

DT 

cue 

C5 

TU 

21 

TBECO 

TRF 

TR 

TPl 

TP2 

22 

TGC 

GK 

CPI 

CP2 

Cl 

23 

C2 

CT (1) 

CT (2) 

CT (3) 

D Time 

24 

GIP 

FP 

AR (1) 

AR (2) 

SFT 

25 

A3 

A4 

A5 

A6 

A7 

26 

A8 

A9 

u>N 

— 

— 

27 

_ — 

Qn 

v F 

vx 

Vg 

28 

<IF 

k! 

kq 

V A 

t e 

29 

v x F 

Kl 

*cc 

Rp 

30 

*K 

*81 

G epf 

MS2D 

— 
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Table D-II 


INITIAL ENVIRONMENT PRINTOUT 


Line 

Position Number 

Number 

1 

2 

3 

4 

5 

6 

1 

DL 

DDV 

Sin AZ 

XPHIO 

V(l) 

V(2) 

2 

V(3) 

Ome Z 

D4(l, 1) 

D4(2,1) 

D4(3,1) 

D4(l, 2) 

3 

D4(2, 2) 

D4(3,2) 

D4(l, 3) 

D4(2,3) 

D4(3,3) 

TFSPR 

4 

TFPR 

R(l) 

R(2) 

R(3) 




Table D-III 

MAIN ENVIRONMENT PRINTOUT 


Line 

Position Number 

Number 

1 

2 

3 

4 

5 

6 

1 

Time 

W 

QX 

Thet B 

Psi B 

Phi B 

2 

Epsil 1 

Epsil 2 

Delta 2 

Thet D 

Yaw D 

Roll D 

3 

TI{1) 

TI(2) 

TI(3) 

VR 

ANP 

AYP 

4 

A(l) 

A (2) 

A (3) 

DE 

SP 

Q 

5 

F(l) 

F(2) 

F(3) 

FA(1) 

FA (2) 

FA (3) 

6 

RA(1) 

RA(2) 

RA(3) 

VA(1) 

VA(2) 

VA(3) 

7 

RP(1) 

RP(2) 

RP(3) 

VP(1) 

VP(2) 

VP (3) 

8 

SR 

CA 

CN 

CY 

XCP 

- 

* 

« R 

«Y 

€ P 

Psi B 

Phi B 

Thet B 


♦Fast-loop printout not shown in Figure D-l. 
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Table D-IV 


DETAILED MAIN PRINTOUT 


Line 

Position Number 

Number 

1 

2 

3 

4 

5 

6 

1 

AT(1) 

AT(2) 

AT(3) 

AG(1) 

AG(2) 

AG (3) 

2 

TV(1) 

TV (2) 

TV(3) 

B(l,l) 

B(2,1) 

B(3,1) 

3 

B(l,2) 

B(2, 2) 

B(3,2) 

B(l, 3) 

B(2,3) 

B(3, 3) 

4 

BT(1, 1) 

BT(2,1) 

BT(3, 1) 

BT(1, 2) 

BT(2, 2) 

BT(3,2) 

5 

BT(1,3) 

BT(2,3) 

BT(3,3) 

RE(1) 

RE (2) 

RE(3) 

6 

«L 

D4(l, 1) 

D4(2,1) 

D4(3,1) 

D4(l, 2) 

D4(2,2) 

7 

D4(2,3) 

D4(3,3) 

D4(2, 3) 

D4(3, 3) 

D4I(1,1) 

D4I(2,1) 

8 

D4I(3, i) 

D4I(1,2) 

D4I(2,2) 

D4I(3,2) 

D4I(1, 3) 

D4I(2,3) 

9 

D4I(3,3) 

VEl(l) 

VEI(2) 

VEI(3) 

VEM(l) 

VEM(2) 

10 

VEM(3) 

TA 

F(l) 

F(2) 

F(3) 

D6(l, 1) 

11 

D6(2,l) 

D6(3,1) 

D6(l, 2) 

D6(2,2) 

D6(3,2) 

D6(l, 3) 

12 

D6(2,3) 

D6(3,3) 

D3(l, 1) 

D3(2,1) 

D3(3,1) 

D3(l, 2) 

13 

D3(2,2) 

D3(2,3) 

D3(l,3) 

D3(2,3) 

D3(3,3) 

D2(l, 1) 

14 

D2(2,1) 

D2(3,1) 

D2(l,2) 

D2(2,2) 

D2(3, 2) 

D2(l,3) 

15 

D2(2,3) 

D2(3,3) 

Thet D 

Yawd 

Rolld 

Alpha 

16 

Beta 

Theta 

PE 

EP 

DELV(l) 

DELV(l) 

17 

DELV(3) 

vi(i) 

Vl/(2) 

VI (3) 

R 1(1) 

Rl(2) 

18 

Rl(3) 

v(i) 

V(2) 

V(3) 

AP(1,1) 

AP(2,1) 

19 

AP(3,1) 

AP(1, 2) 

AP(2,2) 

AP(3,2) 

AP(1,3) 

AP(2,3) 

20 

AP(3, 3) 

RX(1) 

RX(2) 

RX(3) 

GXE(l) 

GXE(2) 

21 

GXE(3) 

— 

— 
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Table D-V 


SYMBOLS FOR GEMINI ASCENT GUIDANCE 
ENVIRONMENT PROGRAM 

(Parentheses indicate maximum dimensions) 


Symbol 

Nearest 
Reference 
Number 
(Figure D-2) 

Definition 

A(3) 

94 (194) 

Total acceleration in I 

AADP(800) 


Angle of attack for digital plot 

AG(3) 

0 

Gravity acceleration in I 

AJF(3, 3) 

45 

Transformation matrix from F to J 

AP(3, 3) 

51 

Transformation matrix from Equatorial 
plane to I 

AR(4) 

3194 

Vehicle area constants for drag 

AT(3) | 

9 

Thrust acceleration in I 

ATT 

3195 

Angle of attack 

A TTY 

3195 

Angle of attack (side) 

AZ 

52 

J to I azimuth angle (equal - Y) 

B(3, 3) 

32 

Intermediate terms in Simpson's Integration 

BT(3, 3) 

32 

Intermediate terms in Simpson's Integration 

C(3) 

17 

Fuel flow rate 

CC2 

17 

Predicted weight for stage testing 

CON4 

2 

Time increment 

CPI 

GUIDl 

Pitch rate from 23. 04 to 89. 6 sec 
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Table D-V. Symbols for GEMINI Ascent Guidance 
Environment Program (cont) 


Symbol 

Nearest 
Reference 
Number 
(Figure D-2) 

Definition 

CP2 

GUID1 

Pitch rate from 89. 6 sec to end of 1st stage 

CPT 

111 

Limit on initial small interval printouts 

CY,CY,CA 

3197 

Aerodynamic coefficients (side, normal, 
axial) 

C2C 

21 

No. of increments near cutoff 

C5 

43 

Angle from east to Izb 

C7 

41 

Constant in geodetic to geocentric conversion 

C8 

41 

Constant in geodetic to geocentric conversion 

DDV 

43 

Longitude difference between launch point 
and ascending node of orbit plane 

DE 

81 

Air density 

DEL 

52 

J to I elevation angle 

DELC 

20 

Smaller increments initiated (near cutoff) 

DELCTl 

18 

Time increment to Stage 1 cutoff 

DELCT2 

19 

Time increment to Stage 2 cutoff 

DELO 

20 

Smaller increment test parameter 

DELT 

22 

One half of the basic time increment 

DELTA 

221 

Basic env. time increment 

DELV(3) 

33 

Intermediate quantity for Simpson's Integra¬ 
tion 

DL 

8 

Intermediate guidance angle 


D-14 




Table D-V. Symbols for GEMINI Ascent Guidance 
Environment Program (cont) 


Symbol 

Nearest 
Reference 
Number 
(Figure D-2) 

Definition 

DTR 

2 

Degrees to radians constant 

DW 

17 

Weight increment 

D3(3, 3) 

52 

J to I transformation matrix 

D4(3,3) 

7 

Spacecraft to I transformation matrix 

E 

111 

Determines interval between printouts 

EXPP 

12 

For detailed printouts 

EZ 

23 

For early termination of program 

FA(3) 

3197 

Aerodynamic Force components in S/C 

F(3) 

515 & 33 

Thrust acceleration in I 

FCY 

111 

Counter for initial printouts 

G(7) 

10 

Constants in gravity series (not used in 
simulation) 

GC 

3197 

Gravity constant 

H(3) 

2 

Beginning of stage weights 

Ip 

112 

Index for graph outputs 

J 

24 

Integer for proper cycling in integration 
loop 

JUP 

16 

Present stage 

KDELTA 

22 

Time Updating integer (to prevent roundoff) 

LDELTA 

22 

Summation of time updating integers 
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Table D-V. Symbols for GEMINI Ascent Guidance 
Environment Program (cont) 


Symbol 

Nearest 
Reference 
Number 
(Figure D-2) 

Definition 

LC1 

133 

Used in staging 

LC2 

136 

For initializing Stage 3 quantities 

MEND 

132 

End program parameter 

MJF 

41 

F to J transformation matrix 

NC 

212 

Cutoff counter 

NR 

214 

Counter for new data read in 

NSD 

214 

No. of new sets of data 

NC2C 

212 

Integer value of C2C 

PE 

8 

Atmospheric pressure 

PfflB 

336 

Roll gimbal angle 

PSIB 

336 

Yaw gimbal angle 

Q 

3194 

Dynamic pressure 

QN 

1 

Counter for early stop for initial debugging 

QPT 

132 

Constant for small interval printouts 

R(3) 

33 

Position vector in I 

RA(3) 

33 

Position vector in GEMINI inertial co¬ 
ordinates 

RL 

51 

Launch site radius magnitude 

RLJ(3) 

51 

Launch site vector in J 

RM 

3191 

Magnitude of position vector 
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Table D-V. Symbols for GEMINI Ascent Guidance 
Environment Program (cont) 


Symbol 

Nearest 
Reference 
Number 
(Figure D-2) 

Definition 

RTD 

11 

Radians to degrees conversion 

RX(3) 

291 

Position components in equatorial frame 

Rl(3) 

331 

Launch position values in I 

SH 

1 

Launch site height above sea level 

SLA 

3194 

Sea level altitude 

SLAT 


Launch latitude 

SLONG 

43 

Launch longitude 

SP 

81 

Velocity of sound 

SR 

3194 

Velocity ratio (mach no.) 

Stage 1 

17 

Stage 1 cutoff time 

Stage 2 

181 

Stage 2 cutoff weight 

t 

221 

Time from launch 

T 

81 

Thrust 

TAS(3) 

3197 

Thrust acceleration in S/C 

TA 

111 

Main printout counter 

TB 

3194 

Second engine thrust (Stage 1) 

THEB 

336 

Pitch gimbal angle 

TIMI 

221 

Initial holddown time 

TTT 

22 

Floating point value of time parameter 
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Table D-V. Symbols for GEMINI Ascent Guidance 
Environment Program (cont) 


Symbol 

Nearest 
Reference 
Number 
(Figure D-2) 

Definition 

TV (3) 

33 

Thrust velocity in I 

Tl(3) 

3194 

1st engine thrust components in S/C 

T2(3) 

3194 

2nd engine thrust components in S/C 

V(3) 

33 

Velocity vector in I 

VA(3) 

33 

Velocity vector in GEMINI coordinates 

VEI(3) 

11 

Velocity of air mass in I 

VEM(3) 

11 

Velocity of air mass in spacecraft frame 

VI 

515 

Previous value of velocity 

VLFU 

51 

Initial earth rate velocity in F 

VM 

11 

Velocity magnitude 

VR 

113 

Relative velocity magnitude 

VI 

33 

Previous value of velocity 

W 

3194 

Total present weight of missile 

XBP(I) 


Pitch attitude error for digital plot 

XO,YO,ZO 

10 

Components of launch site vector 

XINC 

213 

Cutoff time increment 

XINCL 

43 

Inclination angle of desired orbit plane 

XPHIO 

43 

Initial platform roll angle 

YI(8) 

81 

Thrust values for each stage 

Z 

3194 

Height above sea level 

XL, <#>L 

41 

Launch latitude and longitude 

u) 

51 

Earth's rate 

a 

43 

Longitude of ascending node of desired 
orbit plane 
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Table D-VI 


A. 

Main Program 

B. 

Analog Com. 

C. 

Digital Com. 

D. 

Table Look-Up 

E. 

ATMOD 

F. 

OBLATE 

G. 

RUNKUT 
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Table D-VI 


A. Main Program 


C 

C 


R. E. COFER 
JAN. *1963 
RTIM * CLOCKF(QQQ) 

DIMENSION V(3),R(3).D3( 3, 3), D4(3,3),D2(3,3),OME1(3),0ME2(3)* 
AT(3)»D4I(3,3)*VEI(3),VEM(3),BT(3»3), 
B(3»3)*TV(3)»D6t3,3),TVA(3)»VBI(3)»RBI(3)*DVl(8)» 
T2(3)*C(4)* G(7),T1(3)»A(3) ,RE(3), 

TS(3*3)»AG(3)»DV4(12),6ET(3,3),AJF(3,3)» 

VLJl3)»RLJ(3),0ME3(3)»TI{3* 3)* R1(3),DELV(3)* 
AD(2),YI(8),FF{8),SL(8),TTC(3)*H(3)*AR(4), 
VEl3),RA(3),VA(3)»ZP(33),BC(34),DV5(1)»TC(6) 
DIMENSION QZQ(800),XBP(800),AADP(800),V1(3) ,AVEM(3) 

*DV6(3)* DV7(15),F(3),VG(3),CTI 3),VI 1(3),DELTVt3),SF(2)*SF1 
(3),FA(3),TAS(3) 

,RX(3),AP(3,3),GX(3),GXE(3) 

,DV8(3),S2K(11),S2C(17),VP(3),RP(3),DV9(14) 

COMMON SLAT * SLONG* RL» T GLAT * TGLONG,RT »T8,VBI,RBI,AZ,DEL» SWT »DV1 * 

TIME*ALPHAI *BETTAI»THETAI,TVA,THETE,YAWE*ROLLE»R* AT*A,XO* ZO, 
THETD,YAWD,ROLLD,UT,WT,S2K,MSC,MVC,S2C,DELT,D4,G,DV4,AG,YO 
*VT, DL»V*CTfTFSPR»TFPR»DTIME»DI1,DI2,DI3» 

XINCL»OFL»OT.OMEGAE,OMEGAPtVLFU,CS,YO,TU*SS,TT,GK,TDECO 
C»TRF»TR»TP1fTP2»TGC*CPlfCP2fClfC2fG1P»FP f DV6,JUP,F»VG» 

C THEB,PSIB,PHIB,DV7,EXX,EYY,EZZ, 

C EPSILl.EPSIL2.DELTA2,XCP,FA,SFl,SF,SFT,DELTA,DV8,VP t RP 

C *DV9*MS20,DTSC0 

READ INPUT TAPE 5,3,SLAT,SLONG.TGLAT,TGLONG,RT,AZ,DEL,AO ,X0,Z0, 
SF1,OMEGA,C,DEBUG ,SF,C0N4,H,G,TI,E,C3,C4,D5, 
SWT,GC,EZ,EXPP.QPT,TIM I,Y0,BET 
,YI,XINCL,OMEGl,OMEGAP,DELTL,C9,XG,XJ,XK 

,SL( 1),CD,TTC,EO»AD»SH,STAGE 1,NSD,EOK, 
TC,DFL,DT,OMEGAE,C5,TU,TB£CO,TRF,TR,TPl, 
TP2,TGC,GK,CP1,CP2,C1,C2, CT,DTIME,G1P,FP 
,AR(1),AR(21,SFT,S2K»S2C,MS2D 
WRITE OUTPUT TAPE 9,4,SLAT,SLONG,TGLAT,TGLONG,RT,AZ,DEL,AO ,XO, 

ZO,SF1,OMEGA,C,DEBUG ,SF,CON4,H,G,TI,E,C3,C4 
,D5,SWT,GC,EZ»EXPP,QPT,TIM I,YO,BET 
,Y I,XINCL,OMEGl,OMEGAP,DELTL,C9,XG,XJ,XK 

,SL(1),CD»TTC,EO,AD»SH,STAGEl»NSD,EOK, 
TC,DFL,DT,OMEGAE,C5,TU,TBECO,TRF,TR,TPl, 
TP2,TGC,GK,CP1,CP2,C1,C2, CT,DTIME,G1P,FP 
*ARl1),AR(2),SFT,S2K,S2C,MS2D 


3 FORMAT (5E14.8) 

4 FORMAT (5E18.8) 
SENSE LIGHT 1 
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CALL TLU (SR,ATTY,ATT,CY,CN,CA,XCP) 

LC 1= 1 
LC2= 1 
C2C =NC2C 
CPT = QPT 
DO 2 I = 3,8 
2 SL( I ) = 0. 

TB = YII1) 

W = HU) 

DELTA= CON4 
KDELTA = 1 
MEND = 0 
JUP = 1 

DTR = .0174532925 
RTD = 57.295780 
NN = 3 
K = 4 

C7 « .0033841003 
EOC = 1. 

IP = 1 
XBP( 11 = 0. 

AADP(1) « 0. 

RAD = DTR*SLAT 

SLONG = DTR*SLONG 

TGLAT = DTR*TGLAT 

TGLONG = DTR*TGLONG 

AZ = DTR*AZ 

DEL =DTR*DCL 

AO = DTR*AO 

XINCL = DTR * XINCL 

C8 = (1.1731*DTR)/3600. 

DFL = DTR * DEL 
DDV = DFL 
C5 = DTR*C5 
CPI = DTR*CP1 
CP2 = DTR*CP2 

SLAT = RAD - C7 * SINF12. « RAD) + C8 * SINF(4. * RAO) 
SLAT = RAD 
TT = COSF(SLAT) 

SS = SINF(SLAT ) 

RL * SQRTF(X0*X0+Y0*Y0+Z0*Z0) 

SLA= RL - SH 
VLFU = RL*OMEGA*TT 
ON = QN +1. 
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IF (QN - DEBUG) 1,23,23 
1 IF (D13) 43,42,42 

42 CALL GUID1 
GO TO 1 

43 ODV * ODV ♦ (OMEGAE- OMFGAP)• TFSPR 
AA1 » SINF(XINCL) 

SINDL —SS*COSF(XINCL) ♦ AA 1*TT*SINF(DDV) 
OL * ARSINF(SINDLI 
COSDL - COSFIOL) 

SINAZ - -AAl»COSF(ODVI/COSDL 
AZ * ARSINFISINAZ) 

XPHIO * C5 *AZ 
VVV - SINF(XPHIO) 

WWW * COSF(XPHIO) 

OME3I1) - 0. 

0ME3I2) - OMEGA*TT 
0ME3I3) * OMEGA*SS 
Z 1= COSF(DEL) 

Z2* SINF(OEL) 

X = COSFIAZ) 

COSAZ « X 
Y = SINAZ 
03( 1,1)* Z!*X 
03(2,1) * Y 
03(3,1) *-Z2*X 
03(1,2) —Z1*Y 
03(2,2) * X 
03(3,2) * Z2»Y 
03( 1,3) * 12 
03(2,3) * 0. 

D3(3,3) * Z1 
AP( 1,1) ■ SINAZ-SS 
AP(1,2) * COSAZ 

API 1,3) * -SINAZ•TT 
AP(3,1) * *TT 
AP(3,3) * +SS 
AP(2,1) - -COSAZ*SS 
AP(2,2) * SINAZ 

API 2,3) = COSAZ*TT 
45 VV * COSF(OMEGA»TIMI) 

WW * SINF(OMEGA*TIMI) 

AA * SS*TT 
BB * SS»WW 
DO * TT*WW 
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GG = SS*SS 
HH = TT*TT 
A JF( 1,1) = VV 
AJF < 2, 1 ) = BB 
AJF(3,1) =-DD 
AJF( 1,2) = -BB 
AJF(2,2) = HH+VV*GG 
AJF(3,2) * AA*(1. - VV) 

AJF(1,3) - DO 

AJF(2,3) * AA*( 1. - VV) 

AJF(3,3) * GG+HH*VV 
DO 51 1=1,3 

VLJU ) = A JF ( I , I ) *VLFU 

51 RLJU) = AJFU,3)*RL 
DO 52 1=1,3 

OWE HI) = 0. 

DO 52 J=1» 3 

52 OME HI) = OMtKI) ♦ A JF (I, J ) *0ME3( J) 

DO 6 1=1,3 

V11 (I) = VII) 

R(1 ) = 0- 
V CI)~0. 

0ME21I)=0. 

DO 6 J= 1,3 

R(I ) = R(I) ♦ D3(I,J)*RLJ(J) 

V(I) = V(I) * D3(I,J)*VLJ(J) 

6 0ME21I) = 0ME2(I)+D3( I,J)*OMEHJ) 

D0291 I = 1,3 
RXII) = 0. 

D0291 J = 1,3 

291 RXII) = RXII) ♦ AP(J,I)*R(J) 

CALL OBLATE (RX,GXC,1, 1 ) 

D0292 I * 1,3 

AG(I) = 0. 

D0292 J = 1,3 

292 AG1I) = AG(I) + AP(1,J)*GXE(J) 

44 IF (DI2) 516,39,39 

39 IF (LC12) 515,511,511 
511 LC12 = -2 

DO 518 I = 1,3 
518 FII) =-AG(I)*DT 
GO TO 513 
516 LC13 = - 2 
515 DO 512 I = 1,3 
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512 F(I) = F ( I ) + (V(I) - Villi)) -AG(I) * DEL T A 

513 CALL GUID1 
TIKI = TFPR 

IF (LC13) 514,45,45 

514 T I M l = 0 

DO 61 I =1,3 
Rill) = R(I) 

6i vim = vm 

D4( 1, 1) = 0 

D4(2, 1) = 0 
D4(i,l) = 1 
D4(1,2) = + VVV 
04(2,2) =-WWW 
D4(3,2) = 0 

D4( 1,3) = WWW 
D4(2,3) = VVV 
D4(3,3) = 0 

WRITE OUTPUT TAPF 9,13,DL ,DDV , S INA Z , XPHI 0, V, OME 7 ., D4 , TF SPR, 
C TFPR.R 

CALL ATMOO (SH,DE,PC,SP) 

PEI = PE 

81 T = GC •(YI(JUP) *2. )/W 

DO 9 I =1,3 
9 AT( I ) = T *D4(I, 1 ) 

10 CALL COM 

DO 91 I = 1,3 
RX(I) = 0. 

DO 91 J = 1,3 

91 RX( I ) = RX(I) + AP(J,I ) *R(J) 

CALL OBLATE (RX,GXE, 1,1) 

DO 92 I = 1,3 
AG(I) = 0. 

DO 92 J = 1,3 

92 AG(I) = AG(I) + AP(I,J)*GXE(J) 

DO 94 1=1,3 

94 A ( I ) = AT ( I ) -*• AG ( I ) 

VE I ( 1 ) = V(1) - 0MF2(2 )»R(3) + OME 2(3)*R(?) 

VEI(2) = V(2) - 0ME2(3)*R(1) +OME2( 1 >*R(3) 

VEI(3) = V(3) - 0ME2( 1 )*R(?)+0MC2(2)*R( 1) 

DO 11 I =1,3 
VEM(I)= 0. 

DO 11 J =1,3 

11 VEM(I) = VEM( I ) + D4(J,I)•VEI(J) 

DO 113 I =1,3 
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113 AVEM(I) = ABSF ( VFMU) ) 

ATT * ATN1F(AVEM(2),AVFM( 1)1 
AAD = RTD* ATT 
lF(VFM(2)) l|lJ,IIS,HS 

114 ANP = -AAD 
GO TO 1 151 

115 ANP = AAD 

1151 ATTY = ATN1F(AVEM(3),AVFM( 11) 

AADY » R TO * A T T Y 
IF (VEM(3) ) 116,117,117 

116 AYP = AADY 
GO TO 1171 

117 AYP =-AADY 

1171 VM = SORTFIVl 1 )*V( 1)+V(2)«V(2>+V( 31 * V( 3)) 

VR = SQRTF(VEM(1)**2 ♦ VEM(2)**2 ♦ VEM(3)**2) 

FCY = FCY + 1. 

IF (FCY - CPT) 12, 112,111 

111 TA= TA+1. 

IF (TA - E) 133, 112, 112 

112 IP = IP ♦ 1 
AADP(IP) = AAD* 9. 

XBP(IP) = ABSF(EYY*KTD*90.1 
OZQ(IP) = ABSF(THETD*RT0*90.1 

12 WRITE OUTPUT TAPE 9, 1 3 , T I ME , W, OX , TFIFB , P S I B , PH I B , F PS I L 1 , EPS I L2 , 
C DELTA2,THETD.YAWO, ROLLD,T1,VR,ANP,AYP ,A,OF,SP,W,F,FA, 

C t<A,VA,RP,VP,SR,CA,CN,CY,XCP 

13 FORMAT (1H /(1H ,6E 19.81) 

IF (EXPP -2.) 132, 131,131 

131 WRITE OUTPUT TAPE 9,37 , AT , AG , T V ,B , BT,KF,RL,U4,04 I,VEI,VEM,TA 

C F,D6,D3,D2, I HE TD,YAWD,KOLLO,ALPHA,PE IlA 

C 1HETA,PE,FP,DELV,Vl,Kl,V,AP,RX,GXE 

132 TA = 0. 

IF (MEND - 1 ) 153,23,133 

133 IF (LCD 134, 135, 135 

134 JUP = 2 

TB = 0. 

CPT = OPT - 1. 

FCY = 0. 

W = H ( 2 ) 

DELCT3 = DELTA - DELCT1 
DELT = DELCT3/2. 

LC1 = 2. 

FO = 0 
IOC = I. 
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EAK = 0 
GO TO 221 

135 IF (LC2) 136,16,16 

136 JUP * 3 

FCY = 0. 

W = H ( 3} 

0ELCT4 * OELTA - DELCT2 
DELT * DELCT4/2. 

LC2 = 2. 

GO TO 221 

16 GO TO ( 17,181,20,212),JUP 

17 DW *= C( 1)*DELTA 
CC2 = TIME ♦ OELTA 

IF (CC2 - STAGED 22,18,18 

18 LC1 = -2. 

CPT = 1.E6 
GO TO 22 

181 DW = C(2)‘DELTA 
CC2 = W - DW 
IF(CC2 -STAGE2) 19,19,22 

19 LC2 = -2. 

DELCT2 * DELTA*(W -STAGE2l/DW 
DELT = DELCT2/2. 

GO TO 221 

20 DELO * DELTA - DELC 
IF (DELO) 22,21,21 

21 XINC = DELC/C2C 
JUP = 4 

GO TO 12 

212 IF (NC2C- NC1214,214,213 

213 TIME * TIME ♦ XINC 
DELT * XINC/2. 

CPT * 1.E6 
NC = NC ♦ 1 
GO TO 222 

214 IF (NSD -NR 1 23, 23,215 

215 READ 3, YI,SL(1) 

NR * NR*1 

EOC = 1. 

GO TO 81 

22 LDELTA = LDELTA ♦ KDELTA 
TTT = FLOATF(LOELT A) 

DELT = DELTA/2. 

221 TIME=TIMI ♦ DELTA*TTT + DELCT1 +DELCT2 + DELCT3 *DELC T 4 
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227 E Z =EZ + 1. 

IFIEZ-8000.) 24,24,23 

23 CALL STOP 

RTIM = CLOCKF(UQQ) 

CALL PLOT (90.,0.,3.,600 ,3, XBP,AADP,OZQ) 
CALL EXIT 

24 JL= JL+ 1 

GO TO (31,26,26), JL 

26 GO TO (27,27,29), JL 

27 DO 28 1= 1,3 

R(I) = R(I)♦ DELT*V(I) + DELT*B(),I)/4. 

28 V(I) * V(I) + A(I)*DEL T 
GO TO 303 

29 DO 30 1= 1,3 

R ( I ) = RID* DELT * V (I ) + DELT* B (2 , I )/4 . 

30 VII) = V(I) ♦ A(I)*DEL T 
303 RZZ = R(3) 

VEIdJ = VII) - 0ME2( 2 )*RZZ +OME2(3)*R(21 
VEI(2 ) * V(2) - OME2(3)*R(1) ♦0ME2())*RZZ 
VEI(3) = V(3) - 0ME2(1)*R(2)*OME2(2)*R()) 

D030) I =1,3 
VEM(11 = 0. 

D030) J =1,3 

301 VEM(I) = VEM(I) ♦ D4(J,I)*VEI(J) 

VR = SQRTF(VEM(1)«*2 + VEM(2)**2 ♦ VEM(3)**2) 
DO 3195 1=1,3 

3195 AVEM(I) = ABSF(VEM(II) 

ATT = ATN1F(AVEM(2),AVEM(1)1 
ATTY = ATN1F(AVEM(3),AVEM(11) 

QT = R( 1) 

ZQ = R(2) 

ET = Rl31 

D3191 RM = SQRTF(QT**2 +Z0**2 +ET**2) 

DW = C(JUP) * DELT 
3194 W = W - DW 

Z = RM - SLA 
Q X = Z 

CALL ATMOD ( Z,DE,PE,SP) 

VV1 = COSF(EPSIL11 
WWl = COSF(EPSIL2) 

DELTA) = DELTA2 

TUI) = YI(JUP)*VV1*C0SF(DELTA1) 

T1(3) =-YI(JUP)*SINF(EPSILI) 

T1(2) =—YI(JUP) *SINF(DELTA 11 
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T2(1) * TB*WW1*C0SF(DELTA2) 

T2(3) * - TB*SINF(EPSIL2) 

T2(2) » - TB *SINF(DELTA2) 

SR - VR/SP 

CALL TLU (SR,ATTY,ATT,CY,CN,CA,XCP) 

Q »(.5*DE*VR**2) 

FA(1) » -Q*AR(JUP)* CA 
FA(2) - —Q»AR(JUP)• CN 
FA(3) ■ +Q*AR(JUP)* CY 
IF (VEM(2)1 3192,3196,3196 
3192 FA(21 * -FA(21 

3196 IF (VEM(3)1 3199,3198,3198 

3198 FA(31 - —FA 131 

3199 00 3197 I - 1,3 

3197 TAS(I)-(FAlI) ♦ T1(I1 ♦ T2(I))*GC/W 

CALL RUNKUT 
00302 I* 1,3 
AT(11 * 0. 

DO 302 J- 1,3 

302 AT(11 = AT(I) + 04(1,J)* TAS1J1 
00191 I - 1,3 
RX(11 - 0. 

DO 191 J » 1,3 

191 RX(I) * RX(I) ♦ AP(J,I)*R(J1 
CALL 08LATE (RX,GXE,1,1) 

00192 I * 1,3 

AG(11 - 0. 

00192 J « 1,3 

192 AG(11 » AG Cl) ♦ API I,J1*GXE(J) 

D0194 1-1,3 

194 A (II * ATI I) + AG(I1 

31 DO 32 I- 1,3 

BT(JL,I) * AT(I)»DELT*2. 

32 B(JL,I) » A(I1*0ELT*2. 

IFIJL-3) 24,321,345 
321 JL- 0 

00 33 I- 1,3 

DELVIIl = (B(1,1) + 4.»B(2,I 1 ♦ B(3,1)1/6. 

R(I) - R1(I) +(Vl(I) ♦ (B(1,I) +B(2, l) ♦ DELV(I))/6.)*DELT*2. 
V(I) -VKI) ♦ OELV(I) 

DELTV(I) - (BT(1,I) ♦ 4.*BT(2,I1 + BT(3,1)1/6. 

TV(I1- TV(I)+OELTV(I) 

33 F(I) » F(I) + DELTV(I) 

RAM) - R( 1) 
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RA(2) = -R(3) 

RA ( 3) * R (2 ) 

VA ( 1) * V( 1) 

VA ( 21 = -V ( 3 ) 

VAC 3) * V (2) 

DO 331 I =1,3 
VIII) = VCI) 

PSIB * ARSINFC-D4C 2,1)) 

IF (D4(2,2)) 336,335,336 

335 PHIB = 1.5707963 
GO TO 337 

336 PHIB = ATN1F (04(2,31,-04(2,2)) 

337 IF (D4(1,1))333,332,333 

332 THEB » 1.5707963 
GO TO 334 

333 THEB = ATN1F(+04(3, 1 ),04(1, 1) ) 

334 CALL GUI01 

IF (DTSCO) 339,339,338 

338 DELCT2 = DTSCO 
STAGE2 = 500000. 

E = 0. 

EXPP = 2. 

339 GO TO 10 

37 FORMAT (6E18.8) 

345 PAUSE 
END 
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Table D-VI 


B. Analog Com. 


C BOB COFER 

C ANALOG SIMULATION 

SUBROUTINE COM 

DIMENSION DV1(24),DV2(20),DV3( 131),FA(3),SF1(3),SF(2) 

COMMON DV1»TIME,DV2,THETD,YAWD »R0LLD»DV3»EXX»EYY»EZZ» 

C EPSILl,EPSIL2,DELTA2,XCP,FA,SF1,SF f SFT 

IF (LC1) 2,1,1 

1 LC1 = -2 
CALL START 

2 CTM * TIME* SFT 
CEX =EXX * SFltl) 

CEY =EYY * SF1(1) 

CEZ =EZZ * SF1(1) 

CXCP * XCP * SF l (2) 

CFA2 =-FA(3)*SF1(3) 

CFA3 = FA(2)*SF1(3) 

CALL FIXER (CTM,0,CTM) 

CALL FIXER(CEXfO.CEX) 

CALL FIXER(CEY»0»CEY) 

CALL FIXER(CEZ»0»CEZ) 

CALL FIXER(CXCP,0»CXCP) 

CALL FIXER(CFA2»0»CFA2) 

CALL FIXER(CFA3»0»CFA3) 

CALL DATA (7, 0,CTM,CEX,CEZ,CEY,CXCP,CFA2,CFA3,0) 

CALL WAIT 

CALL DATA (0,7,-7,CEPSIL,CDELT1,CDELT2,CEXR,CEYR,CEZR.CEDR ) 
CALL FLTR(CEPSIL,0,CEPSIL) 

CALL FLTR(CDELT1,O.CDELT1) 

CALL FLTR(COELT2,0,CDELT2) 

CALL FLTR(CEXR,0,CEXR) 

CALL FLTRICEYR.O.CEYR) 

CALL FLTR(CEZR»0,CEZR) 

CALL FLTR ( CEDR,0, CEDR ) 

EPSIL1 =CEPSIL * SF(1) 

EPSIL2 = CDELT1 * SF( 1 ) 

DELTA2 =CDELT2 * SFC 1) 

ROLLD = CEXR* SF C 2) 

YAWD = CEZR* SF t 1 ) 

THETD =-CtYR* SF ( l) 

CHECK = -CEDR * SF( 1) 

WRITE OUTPUT TAPE 9,37, CHECK 
37 FORMAT ( 6E18.8 ) 

RETURN 

END 
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Table D-VI 


C. Digital Com. 


C DIGITAL SIMULATION 

SUBROUTINE COM 

DIMENSION DVl(24) t DV2(20)»DV3(131),FA(3),SF1(3i,SF(2) 
COMMON DV1,TIME f DV2,THETD,YAWD,ROLLD,DV3,EXX,EYY,EZZ, 
C EPSIL1,EPSIL2,DELTA2,XCP,FA,SF1,SF,SFT,OELTA 

FAST » FA(2) 

FA (2) = -FA(3) 

FA(3) = FAST 
CALL MC 

YSTORE * ROLLD 
ZSTORE = YAWD 
ROLLD = THETO 
YAWD = YSTORE 
THETD = -ZSTORE 
RETURN 
END 
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Table D-VI 


D. Table Look-Up 

* DICK WADDING 6059 TABLE LOOK-UP AND INTERPOLATION 48 

C 6059 TABLE LOOK-UP AND INTERPOLATION DICK WADDING 

SUBROUTINE TLU (X,AY,AN,Z1,Z2,Z3,Z4) 

DIMENSION TX(50),TA1(10),TA2(10),TZ1(50,10),TZ2(50,10),TZ3(50),TZ4 
1(50,10) 

IF (SENSE LIGHT 1) 1,4 

1 READ INPUT TAPE 5,2,NX,NA1,NA2 

2 FORMAT (2413) 

READ INPUT TAPE 5,3,(TX(I),I*1,NX) 

READ INPUT TAPE 5,3,(TA1(J),J=2,NA1) 

READ INPUT TAPE 5,3,(TA2(J),J=2,NA2) 

READ INPUT TAPE 5,3,((TZ1(I,J),J=2,NA1),1=1,NX> 

READ INPUT TAPE 5,3,((TZ2(I,J),J=2,NA2),I-1,NX) 

READ INPUT TAPE 5,3,(TZ3(I),1=1,NX) 

READ INPUT TAPE 5,3,((TZ4(I,J),J*1,NA2),I * 1,NX) 

3 FORMAT (10E7.3) 

TA1(1) = 0. 

TA2( 1 )=0. 

DO 30 1=1,NX 
TZ 1 (I,11=0. 

30 TZ2(1,11=0. 

4 DO 5 1=1,NX 

IF (X-TXim 6,6,5 

5 CONTINUE 

RETURN 

6 DO 7 J=1»NA1 

IF (AY-TA1(J)) 70,70,7 

7 CONTINUE 

RETURN 

70 DO 71 J2=1,NA2 

IF (AN-TA2U2)) 8,8,71 

71 CONTINUE 

RETURN 

8 F*(X-TX(1-1))/(TX(I) —TX(I— 1 ) ) 

G=(AY-TA 1(J-1))/(TA1(J)-TA1(J- 1 ) ) 

G2*I AN—T A2(J2—1))/(TA2(J2)-TA2(J2-1)) 

ZB1=TZ1(I—1,J—1)+F*(TZ1(I,J-1)—TZllI-1.J-1)) 

ZT1 = TZ1(I—1,J)+F*(TZ1(I»J) —TZl(I-1, J)) 

Z1 = ZBHG*(ZT1-ZB1) 

J=J2 

ZB2=TZ2(I-1,J-1)+F*(TZ2(I,J-1)-TZ2(I-1,J-1)) 

ZT2=TZ2(I—1,J)+F*(TZ2(I,J)-TZ2(1-1,J)) 

Z2=ZB2*G2*(ZT2-ZB2) 

Z3=TZ3(I-1)+F*(TZ3(I)—TZ3(I— 1)) 

ZB4=TZ4(I—1,J-1)+F*(TZ4(I,J-1)-TZ4(I-1,J-1)) 

ZT4=TZ4(1-1,J)+F*(TZ4(I,J)-TZ4(1-1,J)) 

Z4=ZB4*G2*(ZT4-ZB4) 

RETURN 

END 
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E. ATMOD 


C BOB COFER 

C JAN. ,1963 

SUBROUTINE ATMOD(A,KUN,0,P,SS) 
DIMENSION HB( 11),TB( 1I),PB< 11),C(11) 
IF (LC1) 2,1,2 
1 LC1 = 2 

R = 6356.766 
U = .034 1647942 
AA= . I4503776E-3 
B8 = . 1940331 75E-2 
CC= .014837582 
EE= 3.4838395E-3 
C3 = 2.3025851 


HB ( 1 ) 

= 

0. 

HB (2 ) 

= 

11. 

HB (3) 

= 

25. 

HB (4 ) 

as 

47. 

HB (5 ) 

= 

53. 

HB (6) 

= 

79. 

HB (7) 

= 

90. 

HB (8 1 

= 

105. 

HB (9) 

= 

160. 

HB{10) 

= 

170. 

HB(ll) 

= 

200. 

TB ( 1 1 

= 

288. 16 

TB (2) 

= 

216.66 

T B ( 3 ) 

= 

216.66 

TB(4 ) 

= 

282.66 

TB (5 ) 

= 

282.66 

TB (6) 

= 

165.66 

T B ( 7 ) 

= 

165.66 

TB (8 ) 

a= 

225.66 

TB (9) 

= 

1325.66 

TB C101 

= 

1425.66 

TB( 1 l) 

= 

1575.66 

PB ( 1 1 

= 

. 50057165E l 

PB {2) 

= 

.435471 94E1 

PB(3) 

= 

.33959537C1 
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PB(4) = . 20807728EI 
P8<5) = . 1 7658166E1 
PB 1 6) = .40772669E—2 
PB(7) =-.98115386 
PB(8) =-.?1276979EI 
PB{9 ) =-.34412957E1 
PB(10)=-.35492009E1 
PB(11)=-.38460691E1 
C(l) =-.0065 
C(2 ) =0. 

C (3) = .003 
C ( 4 ) =0. 

C C 5) =-.0045 
C C 6) =0. 

C(7) = .004 
C (8) = .02 
C(9) = .01 
C(10)= .005 
C( 11 )= .0035 

2 GO TO (41,42),KUN 

41 AM = .30480061*A 
GO TO 43 

42 AM = A 

43 Z = .001*AM 

H = (R*Z)/CR+Z) 

IF (H—11.) 3,3,4 

3 TP= 288.16 - 6.5*H 
LEV = 1 

GO TO 26 

4 IF IH-25.) 5,5,6 

5 TP=216.66 
LEV = 2 
GO TO 26 

6 IF (H-47. ) 7,7,8 

7 TP = 216.66 + 3.*(H-25.) 
LEV = 3 

GO TO 26 

8 IF (H -53.) 9, 9,11 

9 TP = 282.66 
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LEV = 4 
GO TO 26 

11 IF IH - 79.) 12,12,13 

12 TP - 282.66 - 4.5* 1H-53.) 

LEV * 5 

GO TO 26 

13 IF (H - 90.) 14,14,15 

14 TP « 165.66 
LEV » 6 

GO TO 26 

15 IF (H - 105.» 16, 16, 17 

16 TP - 165.66 ♦ 4.*CH -90.) 

LEV - 7 

GO TO 26 

17 IF (H - 160.) 18,18,19 

18 TP « 225.66 ♦ 20.*(H - 105.) 
LEV - 8 

GO TO 26 

19 IF (H - 170.) 21,21,22 

21 TP - 1325.66 ♦ 10.*(H - 160.) 
LEV * 9 

GO TO 26 

22 IF IH - 200.) 23,23,25 

23 TP - 1425.66 ♦ 5.*CH - 170.) 
LEV * 10 

GO TO 26 

25 TP - 1575.66 ♦ 3.5MH - 200.) 
LEV * 11 

26 V * 1000.*(H-H8(LEV))/TBlLEV) 
IF(C(LEV))27,28,27 

27 V3 - Q/CILEV) 

UU - LOGF(1.+C(LEV)*V)/C3 
ALOGP - PB(LEV) - V3*UU 
GO TO 29 

28 VV - V*CC 

ALOGP - PBCLEV) - VV 

29 PM « 10.**ALOGP 
P =» PM/9.80665 
G -IFF* PM/TP 

SS a 20. Oh;. •’.3 3 * SO" TF ( IP) 

GO- TO (3I,12),KU.M 
’ I P = A A » P M 
0 - fh;*0 

SS a 3.220.133 i* SS 
KLItMN 
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Table D-VI 


F. OBLATE 


SUBROUTINE OBLATE (A,B,N t M) 
DIMENSION A(3)«B(3) 

X=A(1) 

Y=A(2) 

Z=A (3) 

IF (N-1J 4,3,1 

1 IF (N-2) 4,2,4 

2 X=3.2808333*X 
Y=3.2808333*Y 
Z=3.2808333*Z 

3 IF (M) 4,10,11 

10 AJ2=0. 

GO TO 12 

11 AJ2 = . 108228E-02 

12 IE (M-l) 13,13,14 

13 AJ3=0• 

GO TO 15 

14 AJ3=-.23E-05 

15 IF (M-2) 16,16,17 

16 AJ4=0. 

GO TO 18 

17 AJ4=-.212E-05 

18 IF IM-3) 19,19,20 

19 AJ5=0. 

GO TO 21 

20 AJ5=-.2E-06 

21 IF (M—4) 22,22,23 

22 AJ6=0. 

GO TO 24 

23 AJ6=.IE-05 

24 IF IM-5) 9,4,4 
9 CK=. 1407639E+ 17 

A=20925640. 

ZJ2-3./2.*AJ2 
ZJ3=5./2.*AJ3 
ZJ4=5./8.*AJ4 
ZJ5=3./8.*AJ5 
ZJ6=1./16.*AJ6 
ZJ10=3./2.*AJ3 
ZJ 1 1= 15./8.*AJ5 
RH02=X**2+Y**2+Z**2 
RH0=SQRTF(RH021 
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RH03=RH02*RH0 

ZRl=Z/RHO 

ZR2=ZR1**2 

ZR4=ZR2**2 

ZR6~ZR4*ZR2 

AR1=A/RH0 

AR2=AR1**2 

AR3=AR1*AR2 

AR4=AR2**2 

AR5=AR4*AR1 

AR6=AR3**2 

H1 = ZJ2*AR2 

H2=5.*ZR2 

H3=ZJ3*ZR1*AR3 

H4=7.*ZR2 

H5=ZJ4*AR4 

H6=63.*ZR4 

H7=ZJ5*ZR1*AR5 

H8=231.*ZR4 

H9=ZJ6*AR6 

H 10=3003.*ZR6 

RHOX=HI *(l.-H2)*H3*(3.-H4)-H5*(3.-42.*ZR2+H6)-H7«(35.-210.*ZR2*H8) 
*+H9*(35.-945.-ZR2+3465.*ZR4-H10)+l. 

RHOZ=H1*(3.-H21+H3*(6.-H4)-H5*(15.-70.*ZR2+H6)-H7*<105.-315.«ZR2+H 
•81+H9M245.-2205.*ZR2+4851.*ZR4-H101+1. 

RHOZ2=CK/RHO2*(ZJ10*AR3-ZJ1l*AR5) 

H 11=—CK/RH03 
FX = H11*X*RHOX 
FY=H11*Y*RHOX 
FZ=H11*Z*RHOZ+RHOZ2 
B ( 1 ) =FX 
B(2 )*FY 
B(3)“FZ 

IF IN-1> 4,7,5 

5 IF IN—21 4,6,4 

6 FX=FX/3.2808333 
FY*FY/3.2808333 
FZ*FZ/3.2808333 
X=X/3.2808333 
Y=Y/3.2808333 
Z=Z/3.2808333 
GO TO 7 

4 WRITE OUTPUT TAPE 9,8 

8 FORMAT (1H1,40H ERROR IN ARGUMENTS OF SUBROUTINE OBLATE) 

7 RETURN 
END 
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Table D-VI 


G. RUNKUT 


C BOB COFER 

SUBROUTINE RUNKUT 

DIMENSION D4(3,3),DV1(45),DV2(32) 

COMMON DV1,THETD,YAWD,ROLLD,DV2,DELT,D4 
DO 4 I * 1,3 

Cl - D4 (I»2)*THETD - 04(1,31 * YAWD 

D1 * D4(I,3)*R0LLD - 04(1,11 * THETD 

El * D4(I,11*YAWD - 04{1,21 * ROLLO 

A = 04(1,1)+C1 *.5 

B - D4(I» 2 ) +01 *. 5 

C * 04(I» 3)+E 1 *.5 

C2 » B*THETD - C*YAWD 

D2 * C*ROLLD - A*THETD 

E2 * A*YAWD - B*ROLLO 

0 * 04(I,11+C2*.5 

E = D4(I,21+02*.5 

F = 04(I,3)+E2*.5 

C6 * E«THET0 - F*YAW0 

06 * F*ROLLD - D*THETO 

E6 * D*YAWO - E*ROLLD 

R * 04(1,1) ♦ C6 

S = D4(1,2) ♦ 06 

T * D4(1,3) ♦ E6 

C7 * S*THETD - T*YAWO 

07 = T*ROLLO - R*THETD 

E7 * R*YAWD - S*ROLLO 

04(1,1) = 04(1,1) + (Cl+C2*2. +C6*2. +C7)*DELT/6. 

D4(I,2) ■ 04(1,2) ♦ (Dl+D2*2. +06*2. +D7)*0ELT/6. 

4 04(1,3) - 04(1,3) ♦ (El+E2*2. +E6*2. + E7)*DELT/A. 

RETURN 
ENO 
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GUIDANCE EQUATIONS 
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Appendix E 

GUIDANCE EQUATIONS 


I. INTRODUCTION 

This appendix describes the operational IGS guidance equations. 

In general, these equations have been derived in Appendix A and in ref¬ 
erences 1 through 7 listed in Appendix H. 

Certain portions of the guidance equations are presented which have 
specific application to the GEMINI IGS and have not been derived in the refer¬ 
enced documents. A discussion of the IGS equations as implemented in the 
simulation has been included in this appendix to present the simulation oper¬ 
ation in more detail. 

n. DESCRIPTION OF IGS ASCENT GUIDANCE OPERATION 

This section defines the operation of the Ascent Guidance equations dur¬ 
ing launch. Numbers in parenthe s#s refer to areas or blocks in the Ascent 
Math Flow Diagrams (Figures E-1A and E-IB). See Table E-I for the definitions 
of symbols. The platform is aligned as described in Appendix A. The word 
"continuously" when used in the discussion of computer operations means that 
the information is updated at approximately 0. 5-sec-intervals. 


A. FROM START TO PLATFORM RELEASE 

The computer will first be turned to the Ascent mode. All quantities 
which require initialization are contained in block A102. Agena ephemeris 
data will be inserted into the computer approximately 60 sec. prior to launch. 
This information will be continuously updated (A112) and will be used to de¬ 
fine the azimuth orientation of the platform X-axis with respect to East (y) and 
the angle between the orbit plane and the launch site (S). The initial conditions 
on the spacecraft velocity along platform axis (A 117) and the initial platform 
roll gimbal angle (<£ N - A116) will also be computed continuously. 

At approximately 30 sec. prior to launch, the platform will be released 
from torquing. * The guidance system is now in the inertial mode. 


♦As of 28 January 1963, this time - 30 sec. prior to launch - is not definite 
and may eventually be "time of engine ignition." 
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Figure E-IB. GEMINI Computer Math Flow Diagram Ascent 
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Table E-l 


SYMBOLS FOR MATH FLOW 


A. ALPHABETIC (CAPITAL) 


a 3 


Stored constant used to inhibit spacecraft yaw 
angle computation. 


A 4 

A 5 =d 2 r a /dV 2 P 


Stored constant used to limit spacecraft allowable 
yaw angle. 


a 6 = dr a /av 
A7 = d r a /d R 
A 8 = (dra^vj^ 1 

Ag = f drp/dvl -1 


Sensitivity coefficients used for Spacecraft 
Insertion Velocity Adjustment Evaluated at 
perigee. 


Same as above only evaluated at apogee. 


C* Effective exhaust velocity. 

C o Initial eastward velocity at launch point due to 

earth rate. 

C 5 Constant used to define orientation of spacecraft 

y b axis)with respect to east, while vehicle is on 
the lauhch pad. 


C P1> C P2 Constants defining step 1 and step 2 pitch rate of 

the launch vehicle. 

F X» f Y> f Z Integrated acceleration components along platform 

axes. 
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Table E-l. Symbols For Math Flow 


Kj’, K 2 ', K 3 ' Constants used in update of commanded platform 

gimbal angles after the platform is released but 
prior to liftoff. 

K 4 * Gravity component along vehicle velocity vector. 

Kg*, K 7 ', Kg', Kg’ Constants used to compute Stage 1 azimuth angle 

offset to compensate for initial velocity and posi¬ 
tion perpendicular to the orbit plane. 

Kg' Intermediate computer quantity used in calculating 

pitch rate in steering. 

Kg Intermediate quantity used in gravity computation. 

P Commanded vehicle pitch rate. 


Q 5 , Qg, Qn Intermediate coefficients used in steering 

equations. 

R Vehicle distance from earth center. 

Rp. Desired insertion altitude. 

T e Time read from TRS. Will remain zero until 

lift-off. 


t g 


Time to go. Used in Stage 2 steering equation. 
Defines time to effective thrust termination. 


T u 

V 


XL 


Time of update. 

Total vehicle velocity. 

Vehicle velocity perpendicular to the orbit plane. 

Desired insertion velocity. 

Velocity loss term approximating expected loss 
due to gravity. 


Velocity along the vehicle position vector 
(Radial velocity). 
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Table E- 

■1. Symbols For Math Flow (cont) 

V X 

Velocity loss term approximating expected loss 
due to steering (angle of attack). 

V gf 

Final value of gravitational velocity loss term 
remaining at shutdown. 

V\f 

Final value of steering velocity loss term re¬ 
maining at shutdown. 

V XG> V YG> V ZG 

Velocity update components transmitted from 
ground. 

v Xi> V Yi> V Zi 

Measured platform velocity components. 

V XUP. v YUPj 

VZUP 

Platform velocity components interpolated to 
update time. 

Vga 

Horizontal velocity increment required at perigee 
to reach apogee. 

V gP 

Horizontal velocity increment required at apogee 
to reach perigee. 

V’z 

Ground velocity update components corrected for 
azimuth orientation of platform. 

W N 

Intermediate computer quantity used in vehicle 
kinematic computations (W = e AV/C ) 

X, Y, Z 

Vehicle position components. 

Xo> Z G 

Eastward velocity of earth resolved into platform 
frame. 

B. ALPHABETIC (SMALL) 

a l 

Sine X, where X is latitude of launch point. 

a 2 

Cosine X, where X is latitude of launch point. 

a f 

Final value of thrust acceleration. 

a T 

Thrust acceleration now. 
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Table E-l. Symbols For Math Flow (cont) 


*u 

ar T 

g 

gep 

gepf 

gX> g Z 

t 

*B 

*cc 

kx 

tR 

l RF 

l Pl 

tp2 


Value of thrust acceleration at nominal shutdown 
time when W = 1. 

Average value of thrust acceleration between 
nominal and actual shutdown time. 

Average thrust acceleration over previous computation 
cycle. 

Gravity at vehicle altitude. 

Effective gravity - gravity minus centripetal 
acceleration - along geocentric vertical. 

Effective gravity at nominal shutdown condition. 

Gravity components resolved along platform axis. 

Inclination of Agena orbit plane. 

Time after lift-off 

Constant - time bias used to correct for delays in 
receipt and detection of lift-off signal. 

Constant approximately equivalent to the duration 
of computer Stage 2 slow loop. Used in connection 
with the initiation of SECO countdown. 

Time of autopilot gain change following lift-joff. 

Constant used to compensate for any thrust 
acceleration imparted to the vehicle following 
the issuance of the shutdown discrete (mostly 
cut-off impulse). 

Time to start the roll program. 

Time to stop the roll program. 

Time to begin the step 1 pitch program. 

Time to begin the step 2 pitch program. 
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Table E-l. Symbols For Math Flow (cont) 


Computer quantity used to delay entry into 
insertion velocity adjust program and to delay 
turn-off of the SECO discrete. 

Time to start Stage 2 guidance. 

Intermediate quantity used in vehicle kinematic 
computations. 

GREEK (CAPITAL) 

Thrust attitude required to compensate for vehicle 
radial velocity. 

Final value of commanded thrust attitude. 

Computed value of vehicle pitch attitude. Actual 
vehicle pitch attitude with respect to local 
horizontal. 

Commanded vehicle yaw angle (with respect to 
orbit plane). Includes explicit yaw steering. 

Vehicle yaw angle (with respect to orbit plane) 
required to kill velocity perpendicular to the orbit 
plane. 

Thrust attitude required to compensate for 
effective gravity. 

Final value of thrust attitude required for gravity 
compensation. 

Commanded vehicle thrust attitude with respect 
to local horizontal. 

Vehicle flight path angle with respect to local 
horizontal. 

Length of slow loop computation cycle. 

Accumulated time following entry into "SECO 
countdown" loops. Also used to correct initial 
computation cycle time following detection of 
lift-off. 



Table E-l. Symbols For Math Flow (cont) 


at r 

At* 

AV 


AV E 


AY 


^lv^lv^lv 

^°hVO’^LVO> 

LVO 

A0 sc» A< ^sc» A ^sc 


A 0 A</> 

SCO’ 

SCO 


SCO’ 


a^a 


*’o 

*b>*b> *b 
®N> *N> *N 


Time vehicle will roll at a constant rate to reach 
the proper azimuth orientation. 

Fast-loop cycle time (50 msec). 

Total velocity to be gained prior to nominal 
shutdown including approximated velocity loss 
due to gravity and steering. 

Velocity to be gained corrected for actual 
shutdown time cutoff 

( V e = / t a T ) 

Computed value of vehicle yaw attitude. Angle 
between vehicle X-axis and orbit plane. 

Refer to pitch, roll, and yaw, respectively. 

Computed vehicle attitude errors. 

Limited vehicle attitude errors delivered to 
autopilot. 


Computed spacecraft attitude errors during 
insertion velocity adjust. 

Limited vehicle attitude errors displayed to the 
astronaut during insertion velocity adjust. 


Vehicle roll offset required to compensate for 
vehicle position and velocity perpendicular to the 
orbit plane. 

Vehicle pitch attitude error quantity. 

Measured gimbal angles. 

Commanded platform gimbal angles during Stage 
During Stage 2, % and'/'N are equated to actual 
platform gimbal angles once per slow loop. 
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Table E-l. Symbols For Math Flow (cont) 

®FL> ^FL’ ^FL Fast-loop commanded gimbal angles. < This 

includes the effects of the pitch rate (P) term. 

Commanded gimbal rates. 

Coefficient used in computation of steering loss. 

Average value of steering loss coefficient between 
now and shutdown. 

Final value of steering loss coefficient at time of 
shutdown. 

Longitude of vehicle with respect to Greenwich. 
Vehicle yaw attitude error quantity. 

Longitude of ascending node of Agena orbit. 

Rate of earth rotation. 

D. GREEK (SMALL) 

y Prior to platform release, angle between east and 

the platform X-axis. The platform is torqued so 
that its X-axis is parallel to the orbit plane. 
y is positive when X is displaced north. 

y Q Value of y at the time of platform release. 

8 Represents the angle between the launch site and 

the orbit plane; positive when vehicle is below 
orbit plane. 

8y Spacecraft yaw angle required to kill velocity 

perpendicular to the orbit plane. 

Sj Quantity used to allow convergence of certain 

Stage 1 computations upon initiation of Stage 2 
guidance. 

Position increment above or below nominal insertion 
altitude. 


X 

X 

X f 

* 

%' 

n 

fl E 
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Table E-l. Symbols For Math Flow (cont) 


gt Computed quantity which is used in the linear 

interpolation of velocity data. Velocity data is 
corrected to update time. 

8T U Computed quantity representing the adjustment 

to nominal Stage 1 engine shutdown time. 

8V Velocity increment above or below nominal 

insertion velocity. 

gr a Total computed position increment above or 

below apogee. 

grp Total computed position increment above or below 

perigee. 

80 N> Fast-loop gimbal angle increments used to produce 

" desired pitch rate. 

1 Y) % Matrix coefficients used to obtain platform com- 

^ ponents perpendicular to the orbit plane. 

w * Nodal precession rate of the orbit plane. 

w * Pitch rate term used to keep thrust attitude con¬ 

stant with respect to local vertical. 

^p r Pitch rate term used to satisfy vehicle altitude 

constraint. 

w Pitch rate term used to compensate for apparent 

® rotation of gravity vector. 
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B. PLATFORM RELEASE TO LAUNCH 


Following platform release, the computer will commence with the navi¬ 
gation function (A133, A134, A135 and A136). In addition, the computer will 
update platform gimbal angles, which are changing due to earth rotation prior 
to launch (A118). The launch vehicle attitude errors (A256) will be computed 
during this time as well as prior to platform release. These quantities, if 
monitored, could serve to provide some information on the operational readi¬ 
ness of the IGS system. 

In block A106 the computer will be continuously reading the output of the 
spacecraft time reference system (Te) . Any change in this value from zero 
will indicate that lift-off has occurred. The computer will then go into a fast- 
loop (A120, A121 and A122) to obtain the time of lift-off within approximately 
10 msec. 

C. STAGE 1 - OPEN-LOOP STEERING 

The magnitude of the roll maneuver (including offset for vehicle position 
and velocity perpendicular to the orbit plane) is computed (A 125) following lift¬ 
off. This angle in combination with tRF will be used to compute the time to 
start a constant rate (1. 25 deg./sec) roll program (tp),which will bring the ve¬ 
hicle to the required azimuth. 

The vehicle will rise vertically and the computer will test (approximately 
every 0. 5 sec.) for time to start the roll program (A156). At the proper time, 
gimbal angles and rates will be defined in block A161 and attitude errors for 
the launch vehicle will be generated in block A256. Following the completion 
of the constant roll program, the commanded roll gimbal angle will be 
set equal to the value computed for the final roll gimbal angle (<f>i - A162). 

The start of the first- and second-step pitch maneuvers will be con¬ 
trolled by blocks A158 and A159. The time to provide the output for the gain 
change discrete will be controlled by A164. The pitch profile produced is such 
as to approximate a gravity turn, thus minimizing the angle of attack and, 
therefore, the normal forces on the vehicle. During Stage 1 as well as Stage 2 
operation, the platform gimbal rates will be computed in block A168 so as 
to produce the pitch rate desired. The open-loop pitch profile will be continued 
through the staging interval and will be concluded upon computation of the first 
Stage 2 steering commands (A223 and A229). 

D. STAGE 2 - CLOSED-LOOP STEERING 

Block A155 will control the time to initiate the Stage 2 steering computa¬ 
tions. Two passes are used to initialize Wn (A208) in Stage 2. During this 
time, the open-loop pitch maneuver is continued for approximately 1 sec. On 
the third entry into the Stage 2 equations, the commanded attitudes and rates 
obtained from the explicit steering equations are computed. 
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The equations as programmed will steer the vehicle into a plane, de¬ 
fined by the ephemeris data, to the height desired for insertion and will orient 
the velocity vector so as to achieve the desired orbit. When the desired mag¬ 
nitude of velocity is reached, an engine shutdown command (A234) will be 
given. This function is discussed in detail in the following paragraphs. 

The vehicle velocity perpendicular to the orbit plane (Vj_), the velocity 
along the vehicle radius vector from the center of the earth (Vp), and the ve¬ 
hicle position perpendicular to the orbit plane (PjJ will be computed in block 
A203. These quantities will be used in A221 and A222 to compute the vehicle 
commanded angles and rates. Actual pitch attitude of the vehicle with respect 
to its local horizontal and actual yaw attitude with respect to the orbit 
plane (AY) will be computed in A227. The quantities will be used in A229 to 
determine the pitch ( 9 ' 0 ) and yaw ($' Q ) attitude errors of the vehicle. 9' Q and 
\J/q are then inserted into A256 where the fast-loop attitude errors are gen¬ 
erated. The desired vehicle pitch rate (P-Block A223) is used in A168 to ob¬ 
tain the desired gimbal rates which, in turn, are also used in the fast-loop at¬ 
titude error equations. (P, as computed above, does not include any excess 
rate which might be required to bring the vehicle to the commanded pitch at¬ 
titude. ) 

Time-to-go (Tq) is continuously tested in A232. When this quantity is 
reduced to approximately 2 sec.,the attitude errors (A236) will be set to zero 
(thus allowing the vehicle rates to go to zero) and a fast countdown on SECO 
will begin. At the proper time (A232), a SECO signal (A234) will be delivered. 

E. ORBIT VELOCITY ADJUST 

The ascent equations provide a capability to refine the spacecraft veloc¬ 
ity to meet the insertion conditions. The equations do this by using space¬ 
craft energy. The capability is required for the following reasons: 

1. The guidance system may not satisfy these insertion con¬ 
ditions accurately. 

2. Uncertainties associated with residual thrust of the ve¬ 
hicle may exist following insertion. 

3. The payload capability of the booster may fall short of 
the energy required to meet insertion. 

When the payload capability of the booster falls short of the insertion 
conditions, a test is provided on AQ (A205). This test allows entry into the 
orbit velocity adjust equations even when a SECO signal is not delivered. 
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The perturbations from nominal insertion conditions (Vj_, Vp, S y, and 
SR) are computed continuously and form the basic inputs for the orbit velocity 
adjust equations. Approximately 20 sec. (A249, A281) is allowed to elapse be¬ 
fore any commands are generated for spacecraft thrusting. This elapse allows 
the zero attitude error signals (A236) to remain available for the launch ve¬ 
hicle during thrust decay; it also allows the astronaut time to separate the 
spacecraft from the launch vehicle. 

Following this time, the commanded platform roll gimbal angle will be 
set to approximately zero, and the spacecraft attitude errors will be computed. 
The astronaut will first roll the vehicle approximately 90 deg. to null the roll 
error. This will be done in response to the roll attitude error display in the 
capsule. He will then null the yaw and pitch error appearing on the same dis¬ 
play. 


The horizontal velocity will then be computed. This velocity will be 
either added to or subtracted from the spacecraft at perigee to reach apogee 
and at apogee to reach perigee. These quantities will appear on the Av me¬ 
ters in the spacecraft. While nulling out the attitude errors, the astronaut 
will then thrust the vehicle to either add or subtract the velocity appearing on 
the AV indicators. As the spacecraft approaches the desired apogee (A246), 
the velocity to be added or subtracted at perigee will go to zero. When this 
condition is reached, thrusting is discontinued. 

The velocity to be added or subtracted at apogee will be recorded by the 
astronaut. At this point, ascent guidance is essentially concluded. The 
astronaut will use the "Catch-Up" mode of the computer to obtain the velocity 
increment desired at apogee. 

F. PLATFORM UPDATES 

The computer has the capability of accepting velocity data from the 
ground to correct for platform misalignment, and integrated errors in the 
platform (A127, A128, A129 and A138 through A147). The use of velocity 
data in correcting the azimuth orientation of the platform is discussed in Sec¬ 
tion II-B of this appendix. 

The early updates in flight (t < 240 sea) will be used to correct the azi¬ 
muth alignment of the platform; the updates following this time will be used to 
correct the measured platform velocities. 

G. SWITCHOVER FADE-IN 

One additional feature,which is not shown on the math flow but is 
scheduled to be inserted into the computer, is the set of equations used to 
fade in the attitude error signals when switchover occurs from the primary 
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system to the IGS. These equations will have the effect of fading in the atti¬ 
tude errors when a slow drift malfunction occurs. They will also allow the 
major percentage of the attitude error signal to be delivered to the autopilot 
when a rapid malfunction occurs. The result is to allow maximum control 
and response during rapid malfunctions and limited response during slow 
drift malfunctions. 


H. DERIVATIONS 


A. RELATIONSHIP BETWEEN ANGULAR RATES OF GIMBALS, 
SPACECRAFT AND LAUNCH VEHICLE 


Vehicle 



This diagram shows the relation¬ 
ship between the launch vehicle and 
spacecraft reference axes. 


follows: 


Target 


Pb = rotation about body Xaxis - spacecraft roll, positive 
as indicated 

q\j = rotation about body yaxis - spacecraft pitch 
rb = rotation about body Zaxis - spacecraft yaw 


Pm = vehicle roll 


q m = vehicle pitch 
r m = vehicle yaw 
The following relationships exist: 


Pb 

- Pm 

w xb _ 

“xm 

% 

= r m 

Wyb = 

w zm 

r b 

= _( lm 

w zb = 

- w ym 
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The autopilot command signals are documented as follows: To cause 
the launch vehicle to rotate in a positive direction ( + roll, + pitch, and + yaw 
as defined above), the d-c voltage signal to the autopilot must be negative. 

The derivation of the relationship between gimbal rates and body rates is 
equivalent to the development of Euler's kinematical equations with specific 
application to the GEMINI gimbal system. 

It can be shown that 

w xb = ^ sin '!' + <!>, 

<jj yb = 9 cos <f> cos ^ ^ sin <fy , 

and 

w zb = -6 cos ^ sin <f> +*jt cos <f>. 


Also, 


and 


9 sin ^ + <f>, 

6 cos \f/ sin </> - cos <f >, 


w z = 9 cos <f> cos \f/ + ^ sin <f >. 
Multiplying by r, we have 

APm = A0 sin \jr + A<£ , 

Aq m = A9cos\j/ sin<f> - A \f/ cos <f>, 


and 


Ar m = A 9 cos <f> cos $ + A \f/ sin<f>. 

/\p m = Launch vehicle roll change in time r, 
Aq m = Launch vehicle pitch change in time r , 


and 


Ar m = Launch vehicle yaw change in time r . 
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Also 

A 9 = Change in 9 in time r , 

A \J/ = Change in in time r, 

and 

A<£ = Change in <£ in time r. 

To implement a change in gimbal angle, the following notation will be 
assigned to present and desired gimbal angles: 

#b> 'J'b’ ^b = present gimbal angles 

^N> = desired gimbal angles 

Now A 9 = 0 N - 9 b) 

-*b> 

and 

^ = *N - K 

As previously mentioned, a negative polarity is required to cause a 
positive vehicle rotation. Therefore, Ap m , Aq m , and Ar m must be multi¬ 
plied by -1 before being converted to analog signals. This can be accom¬ 
plished by the following operations: 

A0 C =9 b - 9^ 

A *c = ^ - *N 

A< £c = ^b " 

Subscript c denotes commanded change. 

When a specific body angular rate is desired such as a pitch only maneu¬ 
ver during Stage 1, the following relationships are required: 

ro> x = A<£ + A0 sin \ff 

rijy = A 9 cos ^ sin<f> - A if/ cos<f> 

ru> z = A 9 cos cos ^ + A\f/ sin <f> 
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or 


Setting 


and 


A B = 

A f = 
A <f> = 


sin <f> 
cos \j/ 
sin <f> 


cos <f> 

COB\ff 


sin<f> sin ^ 
cos *{/ 


cos <j> sin \j/ 
cos \f/ 


w x = u> z = 0, and « y = P, we have 

. . sin 

6 = P - , 

cos \f/ 

$ = -Pcos <f>, 

$ = - 9 sin \f/. 


B. UPDATE OF COMPUTED ANGLES AND VELOCITY 

Consider a platform which has been erected to the vertical precisely, 
but is misaligned by an angle A 17 about the platform Y axis. 



CLav 


I X) I Y ,Iz = Desired platform orientation. 

IXp, lYp,IZp = Actual platform frame. 

Atj = Misalignment in azimuth, positive 
rotation about platform Y axis 


l Y> J Y p Assume that the ground tracking device can perfectly measure velocity 
and transform the velocity into the desired platform frame. This frame is 
known explicity. 
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V X = Vox + f a X dt > 

t ° 

Vy = jay dt, 

V Z = V 0z + ya z dt. 


V X , Vy, and Vz = actual velocities in desired frame as well 
as measured velocities by ground station. 

V ox and V 0Z = calculated velocity due to earth's rotation 
in the desired frame. 

With a misoriented platform, 

V V 

V Xp = Vq x + / a x cos A t? dt - /az sin At? dt, 


v Yp = J 


a v dt, 


v Zp = V 07 + J a z cos A17 dt + f a x sin A17 dt. 
* o J or 

Here, V 0x and Vq^ are values inserted during initialization. 

If Arj is small and constant, 
cos A17 = 1, 
sin A77 = A?7 , 


V Z p 


°z 


J az dt + A77 J a x dt. 


V Z p - V z = Atj Tax dt = At? (V X - V 0x ), 

q j 

v Z p - V Z 


At? 


v x - v 0x 
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It is now possible to establish the relationship between measurements 
in the desired and actual frames. 


Ixp 


1 

0 -A17 

lx 

I Y p 

- 

0 

1 0 

I Y 

Izp 



0 1 

*z 


This matrix must be used to transform ground velocities into the platform 
frame. 

Vxp' = Vx-AtjVz 
V Y p’ = V Y 
V Zp M = A-r; V X + V z 

Future updates now require that the primed values denoting ground 
measured velocity be transformed to calculated platform axes. Since 1 ?x is 
originally assumed equal to zero, the following relationships are valid for 
the general case of angle update: 


V’Zp 

All 

= %v x + v z 

Vzp - vz'p 

Vxp - Vox 

^Xi 

= %i.i + A*7 


The effect of platform misalignment on initial calculated velocity errors 
(V 0x ' and Vo z ') will now be investigated. 



Initial velocity - Co Ie 
= C Q COS y 


= C 0 cos (y - A 7j) 
= -Co Sin (y- A17) 
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By expanding and making small angle approximations, we have 

V °Xp = Co cos y + Co sin yArj = V ox - V 0z A 17 , 

and 


Vo Zp . = C 0 sin y + Co cos y A 17 = V 0z + y Qx Ar t . 

Vxp and Vzp, however, must be updated to reflect the error in initial values. 
Therefore, 

V XPc = v Xp - v oz Atj , 

and 

Vz Pc = v Zp + V 0x A77 . 

Positions X and Z must also be updated. So, 

X c = X - V 0z A17 T, 

and 


Z c - Z + V 0x A 17 T 

T is the elapsed time from platform release until update. 

C. COMPUTATION OF VEHICLE AZIMUTH AND 
PITCH ANGLE ORIENTATION 


During Stage 2 flight, the actual vehicle attitude is required to imple¬ 
ment the steering commands. The required angles are shown in Figure E-2. 

Rb = I* X + r y Y + r z Z defines the vehicle position in the 
platform frame 


The X-axis of the launch vehicle can be defined with respect to the 
platform by using the following standard GEMINI platform-to-body relation¬ 
ships: 


IXb 


r Yb 

= 

r Zb 



b ij 

Platform to 
Body 


r X 

r Y 

r Z 


and 


I^b = ^ cos 'f'b cos ®b " fy cos '/'b sin ®b " Iz sin 'f'b 


E-23 



Also 


R b Ix b = ll^-ll cos 6 = R sin/3 M , 


and 


X Y Z 

sin/^M = — cos ^ b cos 0 b -cos^ b sin0 b -sin<// b 

R R R 


Similarity, a unit vector perpendicular to the orbit plane has been defined in 
terms of the platform axis: 

r ± = Vx^X + ‘n y *Y +V z Tz 

By using a similar equation, 

sin AY = t) x cos V b cos 8 b - Vy cos ^ b sin 6 b - t) z sin ^ b 



AY is the azimuth angle 
with respect to the or¬ 
bit plane. 


Figure E-2. Vehicle Azimuth and Pitch Angle Orientation 
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D. SUMMARY 


The guidance equation portion of the simulation consists of the follow¬ 
ing subroutines: (1) DI (discrete inputs, (2) FAST (IGS fast-loop attitude error 
computations), (3) GUID 1 (Stage 1 open-loop guidance equations) and (4) 

GUID 2 (Stage 2 closed-loop guidance equations). Fortran listings of these 
programs are included in Table E-II. 

1. Differences Between Math Flow and Equations Used in Simulation 


To date, simulation effort has been primarily devoted to the verifica¬ 
tion of the compatibility of the IGS equations and the vehicle model. The fol¬ 
lowing portions of the IGS math flow remain to be programmed and/or exer¬ 
cised: 


• Test for Lift-off 

• Gimbal Angle Update Following Platform Release 

• Launch Azimuth Offset 

• Platform Velocity Update 

• Orbit Velocity Adjust 

• Stage 2 SECO 

All of the above items will be incorporated into the simulation prior to 
the final simulation report. 

2. Operation of IGS Guidance Simulation (See Figures II-1 and D-2) 

The guidance equations accept ephemeris data from the main program. 
Guid 1 then updates this data while waiting on platform release. Upon re¬ 
lease, Guid 1 begins navigation; inputs to the navigation equations (the measured 
accelerations) are calculated by the main environment. After a nominal 
elapsed time increment, lift-off is simulated and fast-loop attitude error com¬ 
putations are begun. The fast loop accepts the gimbal angles from the main 
environment which is keeping track of vehicle inertial attitude. Commanded 
gimbal angles are generated in either Guid 1 or Guid 2. The attitude errors 
generated in the fast loop are then used as inputs to the simulated autopilot 
and vehicle model which, in turn, send vehicle angular rates back to the main 
program. The main program uses these rates to update the vehicle gimbal 
angles. Subsequently, another fast loop computation cycle may be performed. 

With the exception of a few discretes, no other data need be transferred 
between the environment and IGS equations. 

The fundamental computation cycle time used in the simulation of the 
Guid 1 as well as Guid 1 and 2 equations is 0. 5 sec. This increment may be¬ 
come slightly smaller for the Guid 1 operations and slightly larger for the 
Guid 1 and 2 operations when the computation cycle time of the IGS equations 
as programmed in the GEMINI computer is defined more accurately. 
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Table E-2. Symbols Used In 
Ascent Guidance Simulation (Stage 2) (cont) 


B. VARIABLES 


Symbol 

Env. 

Fortran 

Symbol 

* 

From 

t 

TIME 

E 

8B 

THEB 

E 

'f'B 

PSIB 

E 

N20 

A20 

1 

8* p 

DELEP 

2 

8« y 

DELEY 

2 


ER 

2 

A8 lv 

EP 

2 

A*lv 

EY 

2 

V x 

VP (1) 

1 

v y 

VP (2) 

1 

v z 

VP (3), 

1 

X 

RP (l) 

1 

y 

RP (2) 

1 

z 

RP (3) 

1 

^x 

ETX 

1 


*E - Environment 

1 - 1st Stage Guidance 

2 - 2 n( * Stage Guidance 
K - Constant 


Symbol 

Env. 

Fortran 

Symbol 

* 

From 

Vy 

ETY 

1 

Vz 

ETZ 

1 

g 

GT 

1 

N25 

A25 

1 

N26 

A25 

1 


DTS 

2 

LC29 

LC29 

2 

AtsECO 

— 

2 

P 

PI 

2 


PHIF 

1 

LCA232 

LCA232 

2 

R 

RM 

1 

*1 

- 

- 

PS2D 

PS2D 

K 

LC24 

— 

- 
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c 

c 


BOB COFER 
JAN. ,1963 
SUBROUTINE GUID1 




COMMON 


VG, 


THEB, 


DV1,Dll,012,013, 

XINCL,DFL,OT,OMEGAE,OMEGAP,CO,C5,YO,TU,A1,A 2 .GK.TBFfn 
TRF»TR,TP 1,TP2,TGC,CP1»CP2,C1*C2»G1P»FP, * DVS^JUpfF^* 


«Jf«y: 4 B ^^^:;ii FL,PHiFUDTHEN - o,,siN ->' >HiN - 

v FORMAT (6E18.8) 

IF (LC1) 114,113,113 

113 RRP(2) * —.20909749E8 

RP<2) * 20909749E8 

SINI * SINF(XINCL) 

7 COSI » COSF(XINCL) 

LC1 = -2 

114 N * N+ 1 


r LCA232, 


IF (LCA232) 118,119,119 
119 IF (LC29) 112,117,117 

117 IF (N - 11) 111,112,112 

111 IF (012) 116,115,115 

115 CALL 01 

116 CALL FAST 
RETURN 

118 DT = DTS 

112 N * 0 

71 T » T ♦ OT 

IF (LC18) 73,72,72 

72 T = 0. 

73 IF (LC29) 51,84,84 

84 DFL * DFL + DT * (OMEGAE - OMEGAP) 
B1 * SINF(OFL) 

B2 * COSF(OFL) 

SGM = +SINI*B2 

GM = SGM*(1. - SGM*SGM/6.) 

85 IF (LC26) 13,11,11 

11 IF (013 113,12,12 

12 SOL * —A 1*COSI + A2*SINI*B1 
CGM * 1.-SGM*SGM/2. 
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VPOC 11 = CO*CGM 
VPO(3)=—CO*SGM 

vp( i> » vpom 

VP(3) - VPO{3) 

ETY = SDL 
ETZ = 1. 

PHIN = C5 - GM 
THEN = 1.5707963 
PHIN 1 * 0. 

THEN 1 * 0. 

PSIN1 = 0. 

PHIFL * PHIN 
THEFL * THEN 
PSIFL ■ PSIN 
PHIB » PHIN 
THEB = THEN 
PSIB = PSIN 
GO TO 53 

13 R2 ■ RP( 1)**2 ♦ RP(2)**2 ♦ RP(3)**2 
RM « SQRTFIR2) 

GT =GK/R2 
131 XKG = GT/RM 
DO 18 I = 1,3 
G(I) = XKG *RPCl) 

18 VPKI1 = VP(I) 

FSTORE * F(2) 

F(2) » -F(3) 

F(3) * FSTORE 
152 DO 16 I =1*3 

VP C I > = VP(I) + F(I) - G(I)*DT 

Fpom = fid 

16 F ( I ) a 0 

161 DO 17 I = lf3 

R0FE(I)=R0FE(I)+(VP(I1 ♦ VPI(I))*QT/2. 

17 RP(I) =RRP(I) + ROFE(I) 

QT = DT 

171 IF (LC21) 25,19,19 

19 IF (TU - T) 23,21,21 

21 T1 = T 

DO 22 I = 1,3 

22 VKI) = VP(I) 

GO TO 31 

23 DELT » (TV - T11/DT 
DO 24 I = 1,3 
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24 VUP(I) * VIII) ♦ DEIT*CVP(I)- VI(II) 

TU = TU ♦ 50. 

GO TO 31 

25 LC21 - 2 

IF (308.- TU) 28,26,26 

26 ZDS * VUP(1)*ETX ♦ VUP(2)*ETY ♦ VUP(3)*ETZ 
OPSI *(VG(3) - ZDS)/VUP(1) 

DVPC 1 )* +0PSI*VP0(3) 

DVP(3 )* -DPSI*VP0(1) 

VP(1) * VP(U ♦ DVP(1) 

VP(3) - VP(3) ♦ DVP(3) 

RP(1) = RP(1) ♦ DVP(1) 

RP(3) - RPC 3) ♦ DVP(3) 

ETX - ETX - 0PSI*ETZ 
ETZ = ETZ ♦ DPSI*ETX 
GO TO 31 

28 TU * TU - 40. 

DO 29 I * 1,3 

29 VP(I) * VP C11 ♦ VG(I) - VUP(I) 

31 Z1 - SINF(PHIB) 

Z2 * COSF(PHIB) 

Z4 - COSF(PSIB) 

32 A20 * SINF(PSIB) 

A21 - Z4 * Z2 
A22 ■ Z1 

A23 - Z4 * Z1 
A24 ■ Z2 

A25 = Z4 *SINFITHEB) 

A26 - Z4 *C0SF(THEBI 
IF (LCA232) 322,323,323 

322 LCA232 » 2 
RETURN 

323 IF (TBECO - T) 52,33,33 

33 IF (LC18) 36,34,34 

34 IF (012 ) 35,53,53 

35 LC18 « -2. 

TR = TRF - 45.836624*(C5 -GN -I.570Z963 ♦C1*SDL ♦ C2»CDU 
PHIF * PHIN - (TRF - TR) *.021816616 

36 PHIFL * PHIN 
THEFL * THEN 
PSIFL * PSIN 

361 IF (TR - T) 37,49,49 

37 IF (TRF - T) 39,38,38 

38 PHIN1 =-.021816616 
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PHIN = PHIN + PHIN 1*0T 
GO TO 49 

39 IF C TP 1 -T) 42,41,41 

41 PHIN = PHIF 
PH INI = 0. 

GO TO 49 

42 IF (TP2 -T) 44,43,43 

43 PI = CPI 
GO TO 48 

44 IF (TGC - T) 45,46,46 

45 D057 =» -2 

46 PI = CP2 

48 PS IN 1 =-Pl*COSF(PHIN) 

THEN 1 - Pl*SINFtPHIN)/COSF(PSIN) 

PHIN1 =-THENl*SINF(PSIN) 

PHIN = PHIN + PHIN1*DT 
THEN = THEN + THEN 1*DT 
PSIN » PSIN + PSIN1*DT 

49 DPSIN * PSIN1*DT/10. 

DTHEN = THEN 1*0T/10. 

DPHIN = PHIN1*DT/10. 

GO TO 53 

51 CALL GUID2 
RETURN 

52 CALL GUID2 
GO TO 48 

53 GO TO 114 
END 

* WADDING 6059 GEMINI ASCENT GUIDANCE 48 

C JAN. ,1963 

C 6059 GEMINI ASCENT GUIDANCE-2ND STAGE DICK WADDING 

SUBROUTINE GUID2 

DIMENSION D1(24),D2(5),D3(49),D4(6),D5(4),06(2),D7{17),D1A(25), 

C D3A(31) 

COMMON D1,T,D1A, 

C A3,A4,A5,A6,A7,A8,A9,WN,D2,QN,VF,VL,VG,AFB,CS,CK,VA,TE,V 

1LF,C42,C112,TCC,RF,TK,TS1,GEPF,D3,TAU,D3A, 

C TB»PB,PHB,D4»EN20»05»TPO,PPO,06,OPLV 

20,DTLV0,DSLV0,D7,X1,Y1,Zl,X,Y,Z,ETX,ETY,ETZ,C3,EN25,EN26,DTS,LC29, 
3PHN ,P1 ,PF,LCA232,R,DEL1,MS2D,DTC2 
IF (LC29) 500,40,40 
40 VN2=X1**2+Y1**2+Z1**2 

LCA232=1 
VN=SQRTF(VN2) 

VYP=VY 
VPP=VP 

VY=—Xl*ETX+Y1*ETY+Z1 
9 VP=(X*X1+Y#Y1+Z*Z1)/R 

R1=-ETX*X+ETY*Y+Z 
H1 = R 1 

IF {LC24) 90,99,99 

90 IF (TSO-T) 92,91,91 

91 RETURN 

92 0064=1 

DELV=VF-VN 
DELR=RF-R 

IF t DELV-A3) 94,93,93 
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DVE=DV-AFB*DTU 

QNP=QN 

QN=DV/CS 

DQ=QN-QNP 

IF (ABSF(DQ)-CK) 100,100,103 

100 IF (ABSF(VF-VN)-VA) 102,101,101 

101 SENSE LIGHT 4 

WRITE OUTPUT TAPE 9,1010 
1010 FORMAT (20H ERROR AT BLOCK A206) 
RETURN 

102 LC24=-1 
GO TO 521 

103 WNP=WN 

WN=WNP*(1.+DQ*(l.+DQ«(.5+l./6.*DQ)1) 

C4=C3*SGM 

VG=VG-TAU*C4 

IF (DELI-1.5) 104,104,12 

104 DEL1=DEL1+1. 

GO TO 200 

12 VGF=VG—(TE-TAU)/2.*C4 
DTU={VGF+VLF)/AFB 
DEL3=T—C42 

IF (DEL3) 14,13,13 

13 VLF=VL-BLM*DVE 

14 WS2=VN2*CGM**2/R**2 
WS=VN*CGM/R 
GEP=C3-WS2*R 

Ul® (WNP-WN)/TAU 
AU®U1*CS 

TE=(WN—1.1/U1-0TU 
AT=AU/WN 

AF=AU/(1.+U1*0TU) 

AFB=.5*(AU+AF) 

ATT=AT*{1.-U1/2.*TAU/WN) 

SBG=GEP/AT 

BGNP=BGN 

BGN=SBG*(l.+SBG**2/6.) 

SBB=—VP/DVE 

BB=SBB*(l.+SBB**2/6.) 

93 DELY=VY/DELV 

94 IF (ABSF(OELY)-A4) 98,95,95 

95 IF (DELY) 96,97,97 

96 DELY=-A4 
GO TO 1800 

97 DELY=A4 

98 GO TO 1800 

99 VI* 1./VN 
SG=VP*VI 
CG2=1.-SG**2 
CG®.5*CG2+.5 
GM=SG*(l.+SG**2/6.) 

SGM=SG 

CGM=CG 

10 DV=VF-VN+VL+VG 
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141 Q5=1.-DVE/AF/TE 

Q6=Q5/DVE*CS 
IF IT-C112) 15,15,16 

15 C6=RF-R-TE*VP*(1.-Q6) 
WPR=C6/CS/TE**2/(Q6-.5) 

WYY={-R1+TE*(Q6-1.1*VY)/(CS*TE**2*(Q6-.5)1 

16 WG=2.*WS+GEP/CS 
BYB=VY/DVE 
SS=TE*WPR*IQ6-1.) 

161 BNPP=BGN+BB+SS 

8Y=BYB-TEMQ6-1.)*WYY 

P1=WPR-WS-WG 

BGF=GEPF/AFB 

BPP=B8+TE*Q6*WPR+BGF 

FLAM=.5*(BPP**2+BY*»2) 

XL=.5*{(BNPP-GM)**2 +(BY-VI*Z11**2) 

BLM = . 5* t XL+FLAM) 

VL=VL-XL*TAU*ATT 

1800 BM=ARSINF(X/R*EN26—Y/R*EN25—Z/R*EN20) 
SDT=ETX*EN26-ETY*EN25-EN20 
DT=SDT*l1. + SDT**2/6. ) 

IF (LC24) 171,170,170 

170 PP0=BY-DT 
TP0=BNPP—BM 
PFL=PB 
TFL=TB 
PN=PB 
TN=TB 
PHFL=PHN 

IF ITE-TK-TCC) 50 ,50,501 
501 LCA232= 1 

GO TO 200 
500 0TS=DTS+.05 

50 DTC2=DTS-(TE-TK) 

51 IF (0TC2) 52,53,53 

52 LC29=-1 

521 TS0=T+3.5 

TSl=T+20. 

520 OTLV0=O. 

DPLV0=0. 

DSLV0=0. 

P1=0 . 

GO TO 200 

53 0064=-1. 

LCA232=—1 
LC29=1 
LC24=—1 
GO TO 520 

171 IF (LC4F) 173,172,172 
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172 


LC1* 1 
LC3= 1 
LC8= 1 
LC4F=—1 

173 IF (LC1) 175,174,174 

174 GO TO 200 

175 IF (TS1-T) 176,174,174 

176 PHN=PF-1.5708 
DSSC*DT-DELY 
OTSC=BM 
OPSC*PHB-PHN 
C0C0=.349066 

IF (ABSF(DSSC J-COCO) 178,178,177 

177 DSSC=SIGNFlCOCO,DSSC) 

178 IF (ABSF(OTSC)-COCO) 180,180,179 

179 DTSC=SIGNF(COCO.OTSC) 

180 IF (ABSF(OPSC)—COCOl 182,182,181 

181 DPSC=SIGNF(COCO,OPSC) 

182 DRA=A5*VP**2+A6*DELV*A7*DELR 
ORP=-A5*VP**2-OELR 
DVXB=A8*0RA 

VGA=DVXB 

DVYB=A9*0RP 

VGP=DVYB 

0VZB=0. 

200 IF (MS2D) 201,202,202 

201 RETURN 

202 WRITE OUTPUT TAPE 9,203,LCA232,LC29,LC24,LC1,LC3,LC8,LC4F,VN, 
1,R1,0064,DELV,0ELR,0ELY,GM,DV,DVE,QN,DQ,WN,C4,VG,VGF,DTU,VLF, 
2P,U1,AU,TE,AT,AF,AFB,ATT,BGN,8B,Q5,Q6,C6,WPR,WYY,WG,BYB,BNPP 
3,BGF,BPP,FLAM,XL,BLM,VL,BM,SDT,DT,PPO,TPO,PFL,TFL,PN,TN,PHFL, 
4TC2,TS0,TS1,OTLVO,DPLVO,OSLVO.PHN,DSSC,DTSC,OPSC,DRA.DRP,VGA, 
5VZ8 

203 FORMAT (715/(8E15.7)) 

CALL PDUMP(D1(11,01(500),D 
RETURN 
ENO 


,VY,VP 
, WS,GE 
,BY,P1 
,DTS,D 
VGP»D 
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C BOB COFER 

SUBROUTINE FAST 
DI MENS IONOV1(125},DV2(35) 

COMMON DV1,DI2,DV2, 

C THEB,PS IB» PHIB»THEFL» PSIFL,PHIFL,DTHEN,DPS IN,DPHIN,A20 

C » A21» A22,A23,A24» DELEP,DELEY,DEP»DEY»ER»EP»EY 

IF (D12) 3,4,4 

3 PHIFL * PHIFL ♦ DPHIN 

PS IFL * PSIFL + OPSIN 

THEFL = THEFL + DTHEN 

DELEP = DELEP + DEP 

DELEY = DELEY ♦ DEY 

DPSI =-PSIFL + PSIB 
DPHI =—PHIFL + PHIB 
DTHE =-THEFL + THEB 

ER=+ ( A20*DTHE ♦ DPHI) 

EP=»+( (+A23*DTHE - A24*DPSI) - DELEP) 

£Y=+({A21*DTHE + A22*DPSI) - DELEY) 

2 FORMAT (6E18.8) 

QQ =0Q +1. 

IF (QQ - FP) 4,5,5 

5 WRITE OUTPUT TAPE 9,2, ER,EY,EP, PSIB,PHIB,THEB 

QQ = 0 

4 RETURN 
END 

C BOB COFER 

SUBROUTINE DI 
DIMENSION CT(3),DV1(118) 

COMMON DV1,CT,TFSPR,TFPR,DTIME,DI1,DI2,DI3 
T A= TA+ DTIME 
TFPR = TA -TFSPR 
*F (CT(3) - TA)3,3,8 

3 013 = -2. 

IF (LC1) 4,31,31 
31 TFSPR = TA 
LC1 = -2 

4 IF (CT(1) - TA) 5,5,8 

5 Dll = -2. 

6 IF (CT(2) - TA) 7,7,8 

7 DI2 = -2. 

8 RETURN 
END 
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Appendix F 


DATA CONSTANTS FOR 
GEMINI ASCENT GUIDANCE SIMULATION 


Table F-I is a list of constants used for the GEMINI Ascent Guidance 
Simulation; Table F-II lists the initiatization constants for Stage 2 Guidance. 

Table F-I 


DATA CONSTANTS 


Parameter 

Value 

Footnote 

Launch latitude 

28.46 deg. N 

_ 

Launch longitude 

83.55 deg.W 

1 

Gravitational constant 

32.174 ft/sec? 

2 

Launch point radius magnitude 

20. 909749 x 10 6 ft 

3,7 

Earth's equatorial radius 

20,925,640 xlO 6 ft 

4 

Oblateness, 2nd harmonic 

0.108228 x lO- 2 

4 

Orbital precession rate 

0.145 x 10 -5 radians/sec. 

— 

Earth's rotation rate 

0.7292115 x 10-4 radians/ 
sec. 


Vehicle initial weight, Stage 1 

329,000 lbs 

— 

Vehicle initial weight, Stage 2 

71,500 lbs 

— 

Fuel flow rate, Stage 1 

1662. 1622 lbs/sec. 

— 

Fuel flow rate, Stage 2 

325.41899 lbs/sec. 

— 

Burning time, Stage 1 

148 sec. 

5 

Engine thrust, Stage 1 

215,000 lbs 

— 

Engine thrust, Stage 2 

100,000 lbs 


Vehicle frontal area 

78. 5 ft2 


Initial orientation of vehicle - East to 

94 deg. 


-y 

Orbital inclination 

28.76 deg. 

— 

Central angle between launch longitude 
and longitude of ascending node at 
t = t r ef 

94 deg. 


Time from t ref to platform release 

2 sec. 


Angle between platform Z-axis and 
normal to orbit plane 

0. 2407797 deg. 

6 
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Table F-I. Data Constants (cont) 


Parameter 

Value 

Footnote 

Angle between platform X-axis and 

-1.93737 deg. 

6 

east | 

Time from platform release to launch 

6 sec. 

_ 

Time of end of roll program 

20. 48 sec. (20. 5) 

7 

Time of initiation of pitch program 

23. 04 sec. (23.5) 

7 

Time of pitch rate change 

89. 6 sec. (90.0) 

7 

Time of autopilot gain change 

105.0 sec. 

7 

Pitch rate No. 1 

-0.7727 deg. /sec. 

7 

Pitch rate No. 2 

-0.32 deg. /sec. 

7 

Gravitational constant 

0.14083914 x 1017 

7,3 

Time of change from GUID1 to 

156.04 sec. (156.5) 

7 

GUID2 

GEMINI Computer Comp, cycle time 

0. 5 sec. 

8 

GEMINI Computer Fast-loop Comp. 

0.05 sec. 

7,9 

cycle time 


Notes: 


1. True value 80. 55 deg. W 

2. Value used to convert fuel weight to mass — not used in gravi¬ 

tational acceleration determination 

3. Derived from Fisher Ellipsoid evaluated at latitude of launch 

point 

4. Defined in subroutine oblate 

5. Specified in lieu of Stage 1 final weight for determination of 

Stage 1 cutoff 

6. Computed result from initial conditions 

7. Value in GEMINI Computer — time values go to next 0. 5 sec. 

value 

8. Estimated — chosen for ease of programming 

9. Estimated — equal to 0. 1 of comp, cycle time 
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Appendix G 
VEHICLE DATA FOR 

GEMINI ASCENT GUIDANCE SIMULATION 
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Appendix G 


VEHICLE DATA FOR 

GEMINI ASCENT GUIDANCE SIMULATION 


The data in this appendix represents the composite vehicle model used 
for the Ascent Guidance Simulation. This data is not identical to the latest 
compatibility information because of the preliminary nature of the data avail¬ 
able during the initial stages of the simulation and the difficulty of incorpo¬ 
rating the latest data into the program. 

Table G-I summarizes the vehicle dimensions used in the program; 
Figures G-l through G-6 present aerodynamic and weight-related parameters. 


Table G-l 

VEHICLE DIMENSIONS 


Station of engines - Stage 1 . 1274.121 inches 

Station of engines - Stage 2 .500.00 inches 

Distance from vehicle longitudinal 

axis to engines - Stage 1 .29.07 inches 

Distance from vehicle longitudinal 
axis to roll engine - Stage 2.44.00 inches 
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Figure G-l. Center of Pressure Location vs. Mach Number 
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Mach Number 


Figure G-2. Side Force Coefficient vs. Mach Number 






















































Figure G-3. Normal Force Coefficient vs. Mach Number 
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Mach Number 


Figure G-4. Axial Force Coefficient vs. Mach Number 
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Figure G-5. Moment of Inertia and Center of Gravity vs. Vehicle Weight - Stage 1 
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